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ABSTRACT 

(Distribution  Limitation  Statement  A) 


A  one-dimensional  Lagranglan  hydrodynamics  computer  code  (SAP)  and  a  two- 
dimensional  Eulerlan  hydrodynamics  computer  code  (SHELL)  have  been  successfully 
written  in  the  GLYPNIR  language  for  the  ILLIAC  IV.  Timing  simulations  suggest 
a  speed  50  times  that  of  a  CDC  6600  for  the  GLYPNIR  SHELL  code. 
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SECTION  I 
INTRODUCTION 

The  Air  Force  Weapons  Laboratory  has  been  engaged  in  a  project  with  the 
Advanced  Research  Projects  Agency  (ARPA)  to  provide  applications  support  to  the 
ILLIAC  IV  project.  The  Air  Force  Weapons  Laboratory  (AFWL) ,  as  a  representative 
of  a  wide  variety  of  users,  has  undertaken  a  project  to  look  at  the  ILLIAC  IV 
from  the  point  of  view  of  the  user  and  assist  the  project  in  making  the  connec¬ 
tion  to  the  real  world.  Large-scale  production  codes  from  AFWL  are  being 
adapted  to  the  logic,  languages,  and  operating  system  of  the  ILLIAC  IV.  The 
reprogramming  of  these  codes  provides  an  insight  into  the  usefulness  of  languages 
and  operating  system  of  the  ILLIAC  IV,  and  a  basis  for  timing  comparisons  on 
real  problems.  It  will  also  give  an  indication  of  the  difficulty  in  logically 
modifying  real  programs  to  operate  on  the  ILLIAC  IV  system,  and  provide  an 
interaction  with  user  programmers  to  indicate  the  acceptability  of  the  total 
ILLIAC  IV  system  to  the  user  in  a  real  world  rather  than  an  academic  atmosphere. 
To  this  purpose  AFWL  has  sent  a  number  of  personnel  to  the  University  of  Illinois 
to  work  with  the  staff  there  and  adapt  these  programs. 

This  task  has  been  nontrivial  and  has  involved  considerable  system  and 
program  debugging.  The  purpose  of  the  preserir;  document  is  to  illustrate  the 
method  we  used  to  program  the  SAP  and  SHELL  hydrocodes  for  execution  on  the 
ILLIAC  IV.  Some  basic  knowledge  of  the  ILLIAC  IV,  and  associated  languages 
based  on  the  enormous  amount  of  documentation  already  available  on  the  machine, 
is  assumed.  Machine  details  will  appear  sketchily  where  required. 

The  programming  of  hydrocodes  for  the  ILLIAC  IV  has  involved  several  steps. 
The  first  of  these  was  to  become  acquainted  with  the  ILLIAC  IV  itself,  its 
languages,  and  the  "system"  at  the  University  of  Illinois.  This  latter  item 
included  such  things  as  learning  to  use  the  operating  system  on  the  B5500,  which 
was  at  the  University  during  our  first  visit,  and  the  B6500  which  replaced  it. 
With  this  in  hand  it  was  possible  to  program  a  simple  code  (SAP)  just  to  get 
practice  in  the  use  of  the  languages  and  as  an  exercise  in  testing  the  GLYPNIR- 
ASK-COLLECTOR-LOADER-SIMULATOR  route  that  our  later,  larger  codes  would  have  to 
traverse.  It  also  meant  writing  conversion  routines  to  communicate  with  the 
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simulator  and  its  Simulated  1LL1AC  IV  Disk  files  because  no  one  had  tried  this 
before  and  no  software  existed  for  our  use  in  this  communication. 

After  these  steps  were  accomplished  with  SAP,  the  initial  programming  of 
SHELL  in  GLYPNIR  was  relatively  easy.  The  initial  effort  was  to  write  a 
stripped-down  version  of  SHELL  (SHELL62)  that  would  contain  the  minimum  essen¬ 
tial  elements  for  an  airblast  calculation.  Extensions  could  then  be  made  to 
this  basic  code  to  include  larger  griding  and  other  frills  that  would  not  be 
essential  in  this  first  checkout  version. 

The  GLYPNIR  SHELL  code  was  tested  by  using  the  ILLIAC  IV  simulator  SSK/SSKX 
on  the  B5500  and  then  comparing  the  results  of  this  run  to  that  of  an  identical 
problem  run  using  the  standard  FORTRAN  SHELL  on  the  CDC  6600. 

Extension  of  SHELL62  from  this  minimal  start  to  the  more  complete  versions 
was  quickly  accomplished;  however,  this  job  was  extremely  time  consuming  on  the 
B5500.  The  code  itself  could  take  half  an  hour  of  clock  time  to  compile  and 
another  1-1/2  to  2  hours  to  assemble.  Simulation  of  a  cycle  for  a  grid  5  x  64 
would  then  take  another  hour  or  so  on  the  simulator.  This  then  was  not  an 
atmosphere  in  which  a  real  code  could  be  debugged.  To  reduce  the  number  of 
debug  shots  to  an  absolute  minimum,  a  great  deal  of  extra  labor  was  used. 

When  the  B6500  was  installed  at  the  University  of  Illinois,  the  expansion 
of  SHELL  was  resumed.  The  memory-contained  version  (called  SHELLN)  was  written 
and  tested.  A  version  was  then  written  to  use  the  disk  for  expanded  storage  of 
the  hydro  variables.  This  code,  SHELL/OF/THE/FUTURE,  can  handle  up  to  about 
two  million  cells. 
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SECTION  II 
SAP 

SAP  is  a  simple,  one-dimensional  Lagrangian  hydrocode  used  for  the  calcula¬ 
tion  of  nuclear  and  high-explosive  blast  waves  in  air.  It  is  the  most  elemen¬ 
tary  of  hydrodynamic  calculations,  but  it  does  exist  as  a  full-production  code 
from  which  useful  data  are  currently  being  derived.  We  have  used  this  code  to 
check  out  the  present  state  of  the  ILLIAC  IV  system  and  its  usefulness  for 
scientific  computation.  At  present,  of  course,  the  system  is  simply  a  compiler, 
assembler,  and  simulator  on  the  University  of  Illinois  B5500  (now  B6500) . 
Nevertheless,  this  was  the  beginning  of  the  adaptation  of  much  larger  and  more 
extensive  programs  in  expectation  of  the  ILLIAC  being  available. 

SAP  is  particularly  useful  for  this  purpose  because  it  has  been  the  standard 
first-check  code  for  AFWL  for  a  number  of  different  new  machines,  operating 
systems,  and  equipment.  In  fact,  the  particular  deck  we  started  with  for  this 
project  was  the  one  used  to  check  out  the  extended  core  storage  facility  for 
the  CDC  6600.  For  that  exercise  a  program  which  simulated  the  extended  core  in 
normal  core  was  prepared.  This  was  run  for  logic  verification.  Then  the 
simulated  reads  and  writes  were  replaced  by  system  reads  and  writes  to  extended 
core  and  the  results  were  checked  to  ensure  that  they  were  identical.  Although 
all  the  extended  core  portions  of  this  program  were  removed,  we  mention  this 
here  because  the  general  technique  was  the  same  as  we  followed  in  checking  out 
for  the  ILLIAC  IV.  The  program,  as  it  existed  in  the  particular  card  deck  that 
we  started  with,  had  ample  extended  core  reads  and  writes,  plus  such  features 
as  an  isothermal  sound— speed  calculation  that  is  included  for  purely  historic 
reasons  and,  in  fact,  is  only  used  to  calculate  the  time  step.  SAP  involves  a 
fair  amount  of  coding  which  has  in  the  past  proven  useful  for  checking  out  new 
systems.  We  emphasize  again  that  SAP  is  the  simplest  of  all  possible  hydrocodes, 
yet  historically  it  has  turned  out  to  be  extremely  useful  for  checking  out  new 
machines.  This  simplicity,  of  course,  is  important  because  it  is  easy  to 
discover  where  any  errors  might  have  occurred.  For  all  its  simplicity  SAP  has 
a  sufficient  amount  of  logic  and  a  wide  enough  distribution  of  types  of  arith¬ 
metic  statements  to  exercise  both  the  compiler  and  a  fair  number  of  the  instruc¬ 
tions  that  are  frequently  used  in  any  hydrocode. 
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The  general  technique  used  in  this  and  the  following  code  developments  has 
been  to  first  reduce  the  code  to  the  simplest  possible  configuration  that  will 
run  and  produce  results,  and  from  this  build  up  in  steps,  to  a  full-blown  code. 
SAP  has  an  extensive  amount  of  dump  and  restart  capability  in  it;  it  calculates 
standard  times,  that  is,  certain  selected  times  are  forced  to  the  output  so  that 
comparisons  may  be  made  between  different  runs  at  exactly  the  same  time;  and 
there  is  a  fairly  complex  analytic  equation  of  state,  various  options  for  output 
and  plotting,  including  tapes  and  a  printer  plot,  and  a  fairly  elaborate  auto¬ 
matic  rezone  capability.  All  these  were  removed,  leaving  simply  a  card-read 
input  generation  of  the  problem,  raw  execution,  and  printout  of  the  entire  out¬ 
put  at  every  cycle.  Further,  since  we  were  setting  up  the  code  to  execute  on 
the  ILLIAC  IV  (hopefully  in  parallel)  in  the  simplest  possible  configuration, 
we  reduced  the  number  of  zones  to  64  (actually  62).  This,  of  course,  makes 
absolutely  no  difference  in  the  FORTRAN  version,  but  conceptually  it  is  satisfy¬ 
ing  as  a  starting  point  for  working  with  the  ILLIAC,  therefore  we  actually  did 
make  that  change.  Appendix  I  gives  the  SAP  code  at  that  point— the  so-called 
stripped,  basic  version. 

It  is  important  to  realize  that  the  ILLIAC  IV  has  a  number  of  other  features 
besides  simply  parallel  processing.  Perhaps  an  even  more  important  point  is  the 
fact  that  it  is  really  just  a  number  cruncher  as  far  as  some  code  adaptation  is 
concerned.  It  does  not  have  the  features  normally  associated  with  the  input/ 
output  and  data  handling  of  serial  computers.  One  cannot  reasonably  expect  the 
ILLIAC  IV  to  print  output  or  read  cards.  Its  only  source  of  information  from 
the  outside  world  is  through  the  ILLIAC  IV  disk  (which,  we  understand,  may  be 
changed  in  the  future).  This  must  be  initially  loaded  from  the  driving  computer, 
that  is,  a  core  load  and  binary  data  prepared,  put  on  the  ILLIAC  IV  disk,  read 
into  the  ILLIAC  IV  core,  executed  using  whatever  additional  information  may  be 
available  on  the  disk,  then  written  entirely  to  the  disk.  This  disk  information 
can  then  be  postprocessed  by  the  driving  computer  to  provide  various  visible 
forms  of  output.  The  structure  of  a  code  must  be  made  to  conform  to  these 
limitations. . 

As  it  turns  out,  SAP  and  a  number  of  other  hydrocodes  are  particularly  well 
suited  for  this  adaptation.  Looking  at  the  basic  SAP  code  in  appendix  I,  we 
see  that  it  starts  with  an  input  section  that  reads  cards  and  sets  up  the 
initial  mesh.  It  then  goes  to  a  main  computation  loop  and  ends  in  an  output  or 
print  section.  In  the  more  complicated  code,  of  course,  the  input  section 
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involves  the  reading  of  tapes  and  various  restart  features  and  the  output 
section  includes  writing  of  tapes,  printer  plots,  etc.  Nevertheless,  the 
general  structure  is  still  the  same.  Therefore,  we  chose  to  break  the  code 
into  three  parts  for  convenience  and  labeled  these  program  I  (the  input  section), 
program  IV  (the  center  computational  loop  which  would  be  put  on  the  ILLIAC  IV), 
and  program  II  (the  output  section).  By  input  and  output  we  mean  those  portions 
which  connect  with  the  outside  >vorld,  rather  than  simply  to  the  ILLIAC  IV  disk. 
This  nomenclature  and  this  way  of  thinking  has  been  extremely  useful  in  the 
early  adaptation  of  codes.  Thus  we  have  broken  up  a  code  into  these  three 
constituent  parts  for  running  it  on  the  complete  ILLIAC  system.  The  first 
portion,  program  I  (the  input  section),  and  the  last  portion,  program  II  (the 
output  section),  would  be  run  on  a  serial  machine,  either  the  driving  computer 
or  perhaps  some  other  computer  in  the  ARPA  net,  while  the  main  computation  loop, 
program  IV,  would  take  place  on  the  ILLIAC  IV.  Communications  between  these 
programs  would  be  through  the  ILLIAC  IV  disk.  This  structure  is  not  an  unusual 
procedure  because  a  number  of  standard  hydrocodes  already  use  this  method. 

SHELL,  for  instance.,  which  is  a  two-dimensional  Eulerian  hydrocode,  has  a  very 
extensive  setup  section  known  as  CLAM  and  runs  as  a  separate  program.  Further, 
the  output  processing  is  also  very  involved  because  it  is  contour  plotting  for 
the  most  part,  and  is  also  a  separate  program,  SHPLOT.  The  normal  procedure 
with  SHELL  is  to  prepare  a  tape  with  CLAM;  SHELL  then  reads  this  tape  and 
executes  with  that  data,  writing  out  its  resulJis  on  tape  at  appropriate  times. 
This  output  is  then  postprocessed  with  SHPLOT  and  other  programs. 

After  running  the  basic  SAP  program  in  its  present  form,  which  is  structur¬ 
ally  the  form  we  wish  to  use  on  the  ILLIAC  IV,  we  can  then  adapt  the  central 
section  to  a  language  appropriate  for  the  ILLIAC  IV  and  provide  the  links 
between  this  section  and  the  two  ends  lying  in  the  driving  computer. 

These  two  efforts  proceed  fairly  independently  and  will  be  described  in  no 
particular  historical  order.  Therefore,  let  us  consider  the  problem  of  getting 
the  connection  between  the  exterior  programs  (I  and  II)  and  program  IV. 

There  were  a  number  of  nontrivial  problems  associated  with  this  communica¬ 
tions  task,  and  apparently  very  few  of  these  had  been  faced  to  any  extent  by 
those  previously  working  on  the  ILLIAC  simulator.  We  were  not  adapting  a 
program  for  the  ILLIAC  IV,  but  rather  for  the  B5500  simulator  (SSK/SSKX)  at  the 
University  of  Illinois,  and  to  a  certain  extent  we  had  to  work  within  this 
limitation. 
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The  fact  that  the  B5500  employs  a  48-bit  word  while  the  ILLIAC  IV  and  its 
simulator  both  use  a  64-bit  word  created  one  of  the  main  problems.  Moreover, 
the  general  format  of  the  words  is  entirely  different  between  the  two  machines. 
(Note  that  this  will  continue  to  be  a  problem  in  the  final  machine  configuration 
with  the  B6500  driving  the  ILLIAC  IV.)  The  first  problem,  therefore,  was  to 
provide  a  program  that  would  take  the  initial  setup  data  generated  in  program  I, 
store  it  in  some  common  file,  then  rewrite  this  data  in  a  form  acceptable  to  the 
simulator.  Later,  the  reciprocal  operation  would  have  to  be  performed — the  out¬ 
put  from  the  simulated  ILLIAC  IV  disk  file  would  have  to  be  transformed  back 
into  the  B5500  word  format  for  processing  by  the  postprocessor,  program  II. 

Appendix  VIII  gives  a  short  explanation  of  the  formats  of  both  word  types. 

It  is  noted  that  because  of  the  vastly  different  format  of  the  words,  it  is 
important  to  differentiate  between  fixed-point  and  floating-point  words  when 
making  these  translations.  In  particular,  the  ILLIAC  IV  makes  the  distinction 
between  word  formats  quite  rigorously,  while  in  the  B5500  the  fixed-point  format 
is  just  a  special  case  of  the  floating-point  format;  a  number  carried  in  the 
FORTRAN  program  as  an  integer  may  appear  in  floating-point  format  internal  to 
the  machine.  For  this  reason  it  is  recommended  that  in  the  future  the  records 
being  written  for  the  ILLIAC  IV  core  load  be  entirely  in  one  format— for 
instance,  floating  point.  Minimally,  the  fixed-  and  the  floating-point  variables 
should  be  rigorously  separated  in  the  records  written  or  even  be  put  in  separate 
records.  Presumably,  the  same  caution  holds  for  alphanumeric  or  other  coded 
data.* 

Finally,  an  interesting  problem  related  to  the  number  conversion  arose  that 
was  not  recognized  until  after  all  of  the  other  features  of  the  code  were  well 
checked  out  and  we  were  able  to  compare  and  expect  identical  (at  least  to  near 
the  39-bit  precision  of  the  B5500)  bit  patterns  from  the  FORTRAN  and  GLYPNIR 
results.  When  we  compared  the  results  in  that  fashion,  they  were  not  identical 
and  after  running  a  number  of  cycles,  actual  physical  differences  were  noted. 


*Although  at  the  time  this  was  written  we  could  not  know  what  the  as-yet-to-be- 
written  operating  system  (OSK)  and  the  ILLIAC  control  language  (ICL)  would 
support  in  the  real  ILLIAC  system,  it  is  highly  desirable  that  all  numbers  in  a 
given  file  should  be  of  the  same  form.  The  ILLIAC  control  language  will  support 
the  number  conversion  to  and  from  the  ILLIAC  disk  converting  from  integer,  real, 
or  double  presision  in  the  B6500  to  integer  or  real  in  the  ILLIAC;  however,  all 
of  the  numbers  in  a  given  file  will  receive  the  same  conversion. 
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Now  we  come  to  the  main  part  of  this  section  of  the  report  where  we  will 
describe  the  actual  programming  of  the  ILLIAC  IV  part  of  the  code.  This  was 
our  first  effort,  so  the  procedure  at  this  point  was  extremely  cautious.  We 
took  the  block  that  we  called  program  IV,  the  central  section  of  the  FORTRAN 
code  with  which  we  started,  and  reprogrammed  it  in  FORTRAN  to  adhere  to  the 
logic  of  the  ILLIAC  IV  system.  That  is,  we  restructured  DO  loops  where  neces¬ 
sary  to  range  from  1  to  64  and  isolated  those  portions  of  the  program  which  were 
adaptable  for  parallel  processing.  This  was  essentially  a  trivial  job  in  the 
case  of  this  particular  code,  but  the  technique  is  worth  pointing  out  (and  it 
is  more  interesting  in  the  case  of  a  more  involved  code). 

At  this  point  we  were  in  the  position  to  make  an  almost  card  for  card 
translation  from  the  FORTRAN  into  an  ILLIAC  higher  level  language.  There  were 
two  possible  higher  lever  languages  for  the  ILLIAC  available  at  the  time.  The 
first  of  these,  a  language  having  an  operating  compiler,  is  called  GLYPNIR.  The 
syntax  of  this  language  is  based  on  ALGOL,  therefore  it  looks  a  little  strange 
to  FORTRAN  programmers.  The  logic  is  not  all  that  different  from  FORTRAN  and, 
after  the  initial  shock  wears  off,  one  is  able  to  think  in  this  language  very 
conveniently.  The  GLYPNIR  compiler  produced  an  assembly  language  program  from 
an  initial  set  of  input  cards.  This  program  was  then  put  together  by  the  ASK 
assembler.  The  other  higher  language  available  for  the  ILLIAC  was  a  specially 
extended  FORTRAN  designed  for  the  ILLIAC  IV.  We  conveniently  (but  Improperly) 
called  this  *FORTRAN . 

This  language  was  designed  by  the  University  of  Illinois  but  not  thoroughly 
Implemented.  The  idea  was  to  have  short-order  implementation  of  the  language 
(locally  known  as  COCKROACH)  that  was  essentially  a  *FORTRAN  to  GLYPNIR  trans¬ 
lator  to  gain  initial  experience  with  the  language.  This  would  eventually  be 
replaced  by  full  FORTRAN  compiler  which  would  produce  ILLIAC  code  directly  and 
the  GLYPNIR  compiler  and  the  assembler  would  not  be  used.  Thus,  being  FORTRAN 
users  from  time  immemorial,  we  started  by  transforming  the  simple  FORTRAN  code 
into  what  we  call  the  *FORTRAN  version.  The  reader  attempting  to  use  this 
procedure  will,  of  course,  be  familiar  in  some  considerable  detail  with  the 
syntax  of  *FORTRAN,  so  we  will  not  reproduce  an  explanation  of  the  program  at 
any  length.  The  casual  user  may  wish  to  glance  through  this  program  (shown  in 
appendix  II)  and  quickly  peruse  the  following  short  explanation.  An  asterisk 
in  parentheses  (*)  in  the  position  of  a  subscript  for  a  FORTRAN  variable  implies 
that  an  entire  row  across  the  PEs  is  associated  with  this  particular  variable 
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and  the  arithmetic  statement  there  associated  should  be  done  simultaneously 
across  the  PEs.  This  is,  of  course,  modified  by  the  mode  statement,  the  mode 
being  a  64-bit  word  whose  bits  are  set  to  1  for  those  PEs  you  wish  enabled  and 
0  for  those  PEs  you  wish  disabled.  The  *+l  or  *-l  constructs  simply  refer  to 
the  quantity  in  the  adjacent  PE  to  the  right  or  left,  a  quantity  derived 
eventually  by  routing  in  the  machine.  The  system  is  presently  limited  to  64  in 
a  single  variable- implied  dimension.  It  will  be  extended  to  take  into  account 
variables  being  dimensioned  greater  than  64;  however,  we  shall  hold  this 
possibility  in  abeyance  for  a  while  because  it  appears  that  the  overhead  in 
handling  that  sort  of  thing  in  a  very  general  manner  might  be  excessive  (or 
TRANQUIL-like) .  While  handling  it  in  the  program  explicitly,  it  does  not 
appear  too  much  of  a  problem.  The  program  in  appendix  II  gives  some  idea  of 
the  resemblance  between  this  technique  and  ordinary  FORTRAN.  Even  with  all  the 
asterisks  the  program  does  look  vaguely  familiar.  However,  as  we  shall  later 
see,  this  small  comfort  is  not  entirely  necessary. 

At  the  time  this  program  was  written,  it  was  impossible  to  actually  compile 
the  program  by  machine.  In  some  very  vague  sense  we  translated  it  into  GLYPNIR 
by  hand  and  the  result  of  this  translation  is  shown  in  appendix  III.  Two  or 
three  things  should  be  pointed  out  here.  First,  the  overall  ALGOL  syntax  of 
the  program  is  a  xittle  strange  to  the  FORTRAN-adapted  eye,  but  most  of  the 
arithmetic  statements  and  the  bulk  of  the  logic  does  transform  almost  one  for 
one  from  the  FORTRAN.  There  are  some  major  differences  which  have  to  do  with 
the  declaration  of  variables  at  the  beginning  of  the  code,  i.e.,  the  FORTRAN 
COMMON  structure.  In  this  language,  all  of  the  subroutines  that  we  call  must 
be  placed  previous  to  their  first  use.  Thus  the  lowest  order  subroutines 
appear  first  and  the  main  program  last. 

We  have  here  added  to  the  program  those  portions  which  couple  it  to  the 
ILLIAC  IV  disk,  the  simulated  read  and  write  statements.  One  of  the  features 
that  we  should  bring  up  at  this  point,  which  has  since  been  changed  somewhat 
but  was  necessary  in  the  original  version  of  this  code,  was  that  all  input  had 
to  be  done  with  a  single  long-read  statement  from  the  ILLIAC  IV  disk.  Thus  all 
information  came  in  in  the  form  of  words  written  across  PE  memory.  This  was 
perfectly  appropriate  for  things  like  pressure,  velocities,  and  other  zone 
quantities,  but  it  is  less  useful  for  those  quantities  in  an  ordinary  code  that 
represent  single  variables  rather  than  dimensioned  variables.  These  variables 
will  be  identified  as  CU  variables  in  the  ILLIAC  and  must  be  read  in  as  elements 
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of  a  PE  variable  and  then  set  into  their  respective  CU  variables.  For  this 
purpose  we  used  what  is  known  as  the  Z  block.  A  variable  Z  is  dimensioned  as 
necessary  for  the  total  number  of  variables  we  wish  to  bring  in.  All  the 
appropriate  CU  variables  in  program  I  are  equivalenced  into  the  Z  block  that 
is  normally  used  as  the  first  dimensioned  variable  in  the  program.  One  must 
then  perform  an  action  in  the  GLYPNIR  program  equivalent  to  the  equivalence 
statement  in  program  I.  Thus,  at  the  beginning  of  program  IV  the  individual 
numbers  are  pulled  out  of  the  Z  block  and  put  into  CU  variables  with  a  GRABONE 
statement.  Of  course,  before  every  write  on  the  disk  the  appropriate  CU 
variables  must  be  placed  back  into  the  PE  variable  Z.  (NOTE;  This  is  another 
case  where  floating-  and  fixed-point  numbers  must  be  handled  somewhat  differently 
and  reemphasizes  the  usefulness  of  either  separating  the  variables  or  keeping 
them  all  floating  point.  Our  codes  have  tended  to  use  the  latter  technique.) 
Strangely  enough,  the  Z-block  technique  and  terminology  are  not  unique  to  this 
code  system  or  machine.  It  is  taken  directly  from  the  standard  SHELL  procedure 
to  arrive  at  its  present  form.  Other  than  these  items,  the  GLYPNIR  program  is 
quite  straightforward  and  very  closely  resembles  the  FORTRAN  program.  It  is 
this  that  makes  us  say  that  the  general  structure  of  the  GLYPNIR  is  in  fact  very 
FORTRAN-like,  although  disguised  behind  the  ALGOL  syntax. 

A  number  of  individual,  very  special  problems  turned  up  during  the  checkout 
of  this  code.  Some  misunderstandings  and  difficulties  in  the  compiler  were 
straightened  out  and  these  need  be  of  no  concern  to  the  reader  because  they  have 
been  changed.  There  are,  however,  a  couple  of  points  which  might  be  mentioned 
as  still  being  of  some  use.  The  first  and  foremost  of  these  is  the  use  of  mode 
control  with  the  ILLIAC  because  this  is  something  that  the  programmer  will  have 
to  get  used  to.  It  is  somewhat  tricky  and,  in  general,  we  have  tended  to  play 
safe  in  certain  instances  by  just  setting  the  mode  TRUE  (i.e.,  all  Is)  on  all 
occasions  for  which  we  do  not  have  a  very  specific  reason  for  setting  it  to 
something  else.  Routes,  for  instance,  can  occur  under  node  control,  but  a 
statement  routing  a  row  of  numbers  under  the  current  mode  may  leave  several  PEs 
out  in  the  cold  when,  in  fact,  you  wanted  valuer  brought  up  from  them  and  routed. 
Secondly,  we  note  that  the  ALGOL  IF  ...  THEN  ...  ELSE  statement  is  perhaps 
somewhat  superflous  in  the  context  of  the  parallel  machine  in  that  its  main 
advantage  in  ALGOL  is  to  provide  simply  two  different  routes,  only  one  of  which 
is  executed  during  single  flow  through  the  procedure.  In  a  serial  machine  this 
obviously  is  the  desired  technique.  In  a  parallel  machine,  however,  we  normally 
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will  have  a  number  of  processors  operating  for  one  branch  and  others  operating 
on  the  other  branch.  Moreover,  it  is  found  that  this  particular  GLYPNIR 
construct  generates  a  great  deal  more  actual  code  than  just  setting  the  ELSE 
clause  initially  for  all  PEs  and  then  following  this  with  an  IF  ...  THEN  state¬ 
ment.  This  is  a  small  point,  but  it  explains  why  the  more  obvious  uses  of  the 
IF  ...  THEN  ...  ELSE  statement  was  not  employed  in  this  particular  code. 

This  deck  then,  combined  with  programs  I  and  II  shown  in  appendix  IV, 
represents  the  complete  package  that  was  simulated.  The  code  produced  the 
results  expected  on  comparison  with  the  original  FORTRAN  program. 

One  should  say  something  about  the  actual  physical  operation.  These  remarks 
are  transitory,  of  course,  and  will  change  as  the  operation  changes.  By  the 
time  this  report  is  published,  the  B6500  at  the  University  of  Illinois  will  have 
an  entirely  different  operating  system.  Nevertheless,  the  following  difficulties 
may  still  be  encountered.  First,  the  B5500  expects  a  BCL  character  set  (the 
B6500  accepts  EBCDIC  or  BCL)  and  this,  of  course,  will  make  certain  transforma¬ 
tions  a  little  more  difficult  to  those  of  us  with  only  026  key  punches.  This 
can  be  programmed  ahead  of  time  at  one's  local  installation  and  transformed 
character  by  character  into  an  acceptable  deck.  There  is  one  rather  unfortunate 
feature  in  the  Burroughs  system,  however — the  system  will  stop  the  card  reader 
if  an  illegal  character  is  encountered  during  a  read.  This  is  particularly 
upsetting  in  situations  where  the  error  occurs  in  columns  that  the  program  or 
compiler  will  later  ignore.  We  arrived  with  several  decks  that  had  illegal 
multiple  punches  in  columns  73  through  80.  Card  by  card  by  hand  with  a  key 
puncn  we  had  to  correct  every  one.  It  would  be  a  happier  situation  if  the 
system  ignored  illegal  punches  altogether  and  treated  them  as  blanks  other  than 
in  column  1  (which  signifies  a  control  card),  but  this  is  not  the  case. 

The  system  does  not  accept  serial  batch  runs,  so  the  three  Individual 
programs  had  to  be  loaded  separately,  one  following  the  completion  of  the  other, 
otherwise  the  system  would  attempt  to  execute  a  program  before  its  predecessor 
had  been  completed  and  find  no  data  to  work  on.  Further,  the  printer  unloaded 
regularly,  so  the  compilation  and  the  execution  listing  need  not  be  together  or 
even  in  order.  Certain  limited  types  of  machine  failures  were  common  during 
the  operation  and  it  was  found  advisable  to  be  able  to  restart  the  simulator 
where  the  machine  dropped  off  because  the  simulations  themselves  often  ran  for 
an  hour  or  so.  The  program  that  performed  this  task  is  shown  in  appendix  V. 

Here,  we  pick  up  the  last  completed  cycle  of  our  calculation  from  our  output 
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file,  have  the  capability  to  make  changes  as  desired  to  the  'Z  block,  and  write 
a  new  input  file  to  pick  up  the  simulation  with.  It  is  noted  that  the  more 
straightforward  use  of  the  restart  feature  of  SSK/SSKX  was  not  possible  at  the 
time  we  ran  these  simulations  because  it  lost  track  of  its  disk  files  during  a 
restart. 

The  extension  of  SAP  from  a  straight  64  (or,  in  fact,  62  zone)  mesh  to  one 
of  arbitrary  width  was  straightforward  and  perhaps  fairly  simpleminded.  Never¬ 
theless,  it  works  well.  We  view  the  mesh  as  a  number  of  64-word  blocks,  over¬ 
lapping  one  zone  at  each  end  with  the  next  block,  as  shown  in  appendix  VI.  Due 
to  the  first  order  differencing  scheme  used  in  SAP,  we  must  have  information 
from  the  previous  cycle  about  the  neighbors,  left  and  right,  of  any  cell  before 
we  can  update  it  in  the  mesh.  Thus,  to  update  zone  62  in  this  first  order  code, 
one  must  have  information  from  zones  61  and  63,  both  of  which  are  available  in 
the  first  row.  One  cannot  update  zone  63  with  the  information  available  in  the 
first  row.  We  have  information  about  zone  62,  but  we  do  not  have  information 
about  zone  64.  Zone  63  (the  last  zone  of  the  first  row  since  the  present 
nomenclature  has  been  fixed  with  PEs  running  from  0  to  63),  the  zone  associated 
with  PE  63,  must  therefore  be  repeated  in  the  second  row  so  that  it  can  be 
updated,  but  as  noted  we  will  need  one  zone  to  the  left  of  it.  Thus  the  infor¬ 
mation  from  PE  62  of  the  first  row  must  be  loaded  into  PE  0  of  the  second  row, 

PE  63  into  PE  1,  and  the  second  row  continues  until  the  last  two  Pts  are  reached. 
These  last  two  PEs  must  overlap  the  first  two  PEs  in  the  third  row,  and  this 
continues  for  as  many  rows  as  are  needed.  One  must  point  out  that  while  this 
overlap  must  occur  during  the  run,  the  actual  movement  of  the  numbers  does  not 
occur  until  after  a  cycle  has  been  completed.  That  is,  we  must  use  the  old 
value  of  the  quantity  in  the  62nd  zone  in  PE  0  of  the  second  row  to  update  the 
63rd  zone  in  PE  1  of  the  second  row.  Of  course,  one  can  apply  the  same  sort  of 
thing  or  essentially  create  a  zone  0  for  the  center  boundary  condition  in  this 
spherically  symmetric  code.  These  arrangements  and  scheduling  of  zones  must  be 
made  in  program  I  and  taken  out  in  program  II  before  the  final  upward  processing. 

This  initial  storage  scheme  with  the  data  in  the  PEs  just  described  made 
the  loading  of  these  end  PEs  run  at  extremely  low  efficiency  because  the  data 
from  PE  62  of  each  row  had  to  be  moved  to  PE  0  of  the  next  row.  Similarly, 
data  from  PE  1  of  each  row  had  to  be  moved  back  to  PE  63  of  the  previous  row. 
Thus,  one,  or  two  at  the  most,  PEs  could  be  active  during  this  adjustment  phase, 
making  it  much  less  efficient  than  the  actual  calculation  where  62  or  64  PEs 
were  active  at  a  time. 
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This  was  a  very  crude  way  of  handling  data  and  a  much  morn  sophisticated 
method  immediately  came  to  mind.  For  instance,  one  can  rotate  the  second  row 
two  PEs  to  the  left  with  respect  to  the  first,  the  third  with  respect  to  the 
second,  by  an  additional  two  PEs,  etc.  In  this  way  one  can  make  the  entire 
readjustment  process  in  parallel  with  all  PEs  active,  rather  than  making  a 
separate  adjustment  for  each  row.  This  is  perhaps  a  trivial  additional  compli¬ 
cation  but  one  which  we  did  not  wish  to  get  into  the  first  time  through.  SAP, 
although  a  simple  code,  has  many  of  the  features  of  the  more  complex  codes,  and 
because  of  its  simplicity,  it  allowed  major  charges  in  programming  and  storage 
with  only  a  minimal  rewriting  of  code.  The  coding  of  SAP  brought  out  many  of 
the  difficulties  of  the  present  operating  system  and  allowed  us  to  exercise 
solutions  thereto.  It  also  provided  a  great  deal  of  experience  in  actual 
physical  manipulation  of  codes  in  the  machine. 

This  version  of  SAP  is  certainly  not  the  west  general  code  even  for  its 
particular  limitations  to  one-dimensional  spherical  hydro.  The  main  production 
codes  are  normally  run  with,  for  instance,  complex  equations  of  state.  The 
equations  of  state  normally  used  in  this  code  are  interchangeable  with  other 
codes  (like  SHELL)  and  were  written  as  a  part  of  that  effort.  Some  of  the  more 
elaborate  rezone  techniques  have  not  been  included  in  this  version.  Because 
we  believe  these  are  really  details  and  do  not  further  illustrate  our  main 
purpose  in  this  initial  run,  we  did  not  include  any  of  the  fancier  additions 
such  as  radiation,  complex  boundary  conditions,  etc.  We  did  not  face  some  of 
the  restart  problems  because  they  will  be  very  specifically  related  to  the 
driving  machine  and  necessary  details  were  not  available  for  our  investigation. 
Many  of  these  problems  will  be  taken  care  of  in  later  versions  of  the  code  and 
perhaps  one  of  these  later  versions  of  SAP  will  be  brought  up  to  full  strength 
immediately  before  the  full  availability  of  the  ILLIAC  IV  itself. 

Finally,  for  comparison  purposes  the  ASK  code  assembly  language  which  was 
produced  by  GLYPNIR  was  compared  with  an  assembly  language  program  that  was 
written  from  the  same  groundwork,  but  hand-coded,  rather  than  using  the  GLYPNIR 
compiler.  This  might  conceivably  be  looked  upon  as  what  a  better  compiler  should 
be  able  to  produce.  It  was  noted  that  this  hand-compilation  gave  a  fair  reduc¬ 
tion  in  number  of  instructions  generated  (a  factor  of  3),  but  that  the  improve¬ 
ment  was  not  especially  startling  in  execution  time.  One  can  say  that  for  a 
first  cut  the  GLYPNIR  compiler  generates  quite  respectable  code. 
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The  pure  assembly  language  program  Is  listed  in  appendix  Vll 
compiled  and  executed  as  a  part  of  the  three-program  operation. 
Identical  results  to  the  other  procedures. 


and  was 
It  gives 
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1.  THE  METHOD 

With  the  learning  experience  of  SAP  behind  us,  we  initiated  an  effort  to 
convert  the  much  more  extensive  SHELL  hydrocode  to  the  ILLIAC. 

This  effort  proceeded  along  two  lines.  On  the  one  hand,  Major  Whitaker  ar.d 
Mr.  Needham  took  a  FORTRAN  SHELL  listing  and  in  the  spirit  of  the  SAP  effort, 
wrote  a  SHELL  version  from  it.  This  then  could  be  held  until  the  COCKROACH 
translator  was  available.  This  also  could  be  used  as  a  stepping  stone  in  the 
translation  of  the  serial  FORTRAN  logic  to  parallel  GLYPNIR  logic.  This  effort 
produced  a  nearly  complete  * SHELL  deck  for  rows  up  to  62  zones  wide.  This  code 
was  later  run  through  the  *FORTRAN  translator  but  was  not  compiled  or  executed. 

2.  FORTRAN  SHELL 

One  is  appalled  by  the  large  amount  of  code  in  a  working  FORTRAN  version  of 
SHELL.  This  is  traceable  to  the  limitations  of  the  various  machines  and  com¬ 
pilers  that  SHELL  has  run  on.  Once  added  to  SHELL,  this  code  has  never  been 
removed.  Indeed,  SHELL  works  as  it  stands.  In  any  removal  process  the  most 
careful  attention  to  detail  would  be  required  to  ensure  that  necessary  and 
currently  useful  data  are  not  deleted  inadvertently. 

Leaving  everything  that  has  accumulated  over  the  years  in  the  code  is  quite 
acceptable  when  moving  the  code  from  one  FORTRAN  compiler  to  another  because 
extensive  time  and  effort  at  a  keypunch  would  be  needed  to  remove  those  parts 
of  the  code  which  are  of  questionable  use.  Only  a  few  minor  changes  are  usually 
necessary  to  get  the  entire  code  running  on  a  new  compiler.  This  is  hardly  the 
case  in  reprogramming  the  code  for  a  new  language  (GLYPNIR)  and  a  new  machine 
architecture  (the  ILLIAC  IV).  Since  the  code  had  to  be  rewritten  from  scratch 
and  carefully  debugged,  this  was  the  most  reasonable  time  to  delete  much  of  the 
overhead  from  SHELL  and  make  the  programming  task  in  GLYPNIR  at  least  smaller. 

In  rewriting  the  code,  it  was  decided  to  remove  the  capability  (such  as  it 
was)  to  handle  vacuum  zones  because  this  only  gives  stopgap  answers  to  the 
questions  one  tries  to  solve  by  using  this  technique.  The  artificial  viscosities 
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were  also  removed  from  PHI  of  the  code  because  these  have  not  been  used  in  the 
FORTRAN  version  of  SHELL  in  many  years.  The  viscosity  calculations  almost 
double  the  number  of  lines  of  code  in  PHI. 

The  INPUT  and  EDIT  routines  were  moved  to  the  driving  computer  programs  I 
and  II  and  with  the  deletions  the  length  of  SHELL  was  reduced  from  about  2500 
to  4000  cards  in  its  various  FORTRAN  versions  to  approximately  500  to  700  cards 
in  GLYPNIR.  Thus  substantial  deletions  to  PHI  and  PH2  (artificial  viscosities, 
vacuum  zones)  dictated  that  we  start  from  scratch.  If  we  had  tried  to  adapt 
the  code  as  it  stood,  the  known  inconsistencies  in  the  code  that  handles 
variable  zoning  in  the  mesh  could  have  caused  trouble.  In  this  stripped-down 
version,  these  inconsistencies  could  be  corrected  easily. 

Our  approach  was  first  to  write  a  stripped-down  version  of  the  code  and 
work  up  from  there.  We  decided  to  rederive  the  basic  difference  equations  on 
which  SHELL  is  based.  This  rederivation  crystallized  our  thoughts  and  helped 
to  bridge  the  gap  between  the  serial  logic  in  the  FORTRAN  version  of  the  code 
and  the  actual  equations  we  were  trying  to  solve.  This  approach  seems  highly 
desirable  for  many  large  codes  being  transferred  to  the  ILLIAC  IV. 

At  this  stage  the  problems  to  be  encountered  by  a  computer  translation  of 
a  serial  code  into  a  parallel  code  (or  indeed  the  translation  of  a  large  code 
by  someone  who  is  not  familiar  with  it)  became  only  too  apparent.  The  programmer 
who  is  familiar  with  the  physics  of  the  code  could  jnake  a  number  of  changes, 
deletions,  etc.,  but  the  programmer  trying  to  make  a  strict  one-for-one  trans¬ 
lation  of  the  code  would  be  forced  into  a  number  of  costly  inefficiencies  if 
he  wished  to  preserve  the  code.  What  is  basic  to  SHELL  is  perhaps  a  half  dozen 
lines  of  physics  presented  in  the  next  section  and  not  the  thousands  of  cards 
of  superstructure  into  which  they  are  embedded  so  as  to  apply  them  to  a  problem. 

Now  we  will  briefly  sketch  the  rederivation  of  the  SHELL  difference  equa¬ 
tions  as  made  by  Captain  Durrett.  Except  for  the  notation,  which  is  now  much 
clearer,  there  are  no  changes  in  the  basic  equations  from  those  that  were 
derived  by  others  at  earlier  times.  There  are  differences  in  the  boundary 
conditions  derived  and  in  the  update  to  the  total  energy  in  the  system  (ETH) 
that  were  discovered  in  this  rederivation.  Earlier  authors  had  made  those 
cosmetic  changes  necessary  to  convert  the  constant  zoning  results  to  variable 
zoning  results,  but  had  missed  several  dimensionless  multipliers  which  would 
have  been  included  if  a  consistent  rederivation  had  been  done  with  variable 
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zoning  in  it  from  the  start.  Thus  it  may  be  that  the  poor  relative  error  that 
SHELL  has  reported  over  the  years  while  using  a  variably  zoned  mesh  is  really 
an  error  in  the  computation  and  reporting  rather  than  an  actual  error  in  the 
calculation. 

3.  THE  SHELL  DIFFERENCE  EQUATIONS 

The  basic  hydrodynamic  equations  integrated  by  SHELL  are 


/  ^ 

\ 

i  4 

2- 

+  u.V 

p  +  p  V.u  =  0 

(1) 

\  3t  / 

-v  ) 

k  4 

P  /I- 

+  u.V, 

u  +  Vp  =  0 

(2) 

\  3t  f 

P  (f 

\3t 

+  u.vj 

|  E  +  V*  pu  -  0 

(3) 

P  =  P  (E,p) 

(4) 

Taking  equations  (2)  and  (3)  and  ignoring  the  (u*V)  operator,  we  get 
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Equation  (5)  gives  us  one  of  our  difference  equations  directly.  Equation 
(6)  must  be  manipulated  somewhat  before  it  can  be  used.  The  E  in  equation  (6) 
is 

E  “  +  0.5  (u*u) 

i.e.,  specific  material  plus  specific  kinetic  energy.  Therefore, 
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We  now  expand  the  divergence  in  equation  (6)  to  get 

9E  +  1^  (pV.u  +  u  .  7p)  =  0 

9t  p 

and  then  use  equation  (7)  to  obtain 

m  -  /V£\.  U  +.  pV  .u  +  u.Vp  B  0 
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or  canceling  terms 
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This  can  be  expanded  in  two-dimensional  cylindrical  coordinates  to 
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This  then  is  the  second  equation  we  want  to  solve. 

The  finite  differencing  of  equations  (5)  and  (8a)  is  straightforward. 
Equation  (5)  becomes,  in  cylindrical  coordinates  and  the  nomenclature  of 
figure  1, 


Equation  (8a)  becomes 
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where 


r 


Figure  1.  The  SHELL  grid 

The  u,  v,  and  En  are  the  new  velocities  and  internal  energy  at  time  (n+1)  based 
on  the  pressure  effects  of  equations  (2)  and  (3)  only  before  accounting  for  any 
mass  motion. 

The  distance  r1  refers  to  the  right  edge  of  the  cell  k  while  rk  refers  «_o 
the  center  of  the  cell,  i.e., 

r,  =  ri  +  ri-l 
k  - 2 - 

Similarly,  the  pi,  u^  and  are  defined  at  cell  boundaries  while  u^,  v^,  pk, 

p, ,  and  E  are  defined  at  cell  centers.  Thus,  for  constant  zone  size  p . ,  for 
k  m  i 
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example,  is  defined  by 


pk  +  pk+l 


In  the  case  of  nonconstant  zone  size  this  becomes 


4rk+i  +  Vi  4rk 

4rk+l  +  4rk 


where 


Now  this  choice  of  u  and  v  Is  used  ir  the  energy  equation  to  achieve  energy 
conservation;  however,  in  the  case  of  variable  zone  size  the  conservation  is 
no  longer  assured.  SHELL  for  the  ILLIAC  was  programmed  to  allow  for  variable 
zone  size. 

Equations  (9) ,  (10) ,  and  (11)  are  solved  in  and  constitute  what  is  called 
phase  1  of  the  calculation.  Equation  (1)  Is  solved  In  what  is  called  phase  2, 
as  follows. 

Rewriting  equation  (1) 


+  V*  pu  =  0 
3t 

or  in  two-dimensional  cylindrical  coordinates 
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the  finite  differencing  of  this  equation  gives 
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SHELL  uses  this  foraulatlon  except  that  the  p^  and  p^  are  taken  not  on  the 
boundary  as  Indicated  in  equation  (12),  but  rather  the  density  chosen  is  that 
of  the  center  of  the  cell  donating  the  mass.  Also,  the  velocity  chosen  is  not 
precisely  u^  or  as  indicated,  but  rather  a  second  order  choice  is  made.  This 
is  done  to  improve  the  behavior  of  the  method  in  areas  behind  a  strong  shock. 

The  velocity  chosen  is  called  the  "rear  of  slug  velocity"  and  is  shown  in  figure 
2.  The  volocity  u|  is  chosen  as  a  linear  Interpolation  between  u^  and 
to  the  point  r  such  that  with  velocity  u^  the  rear  of  the  slug  will,  in  time 
At,  move  exactly  to  position  r^,  the  boundary.  The  amount  of  mass  between 
positions  r^  and  r^  will  therefore  move  into  cell  k+1  in  the  current  time  step 
(at  a  density  of  p^) .  With  these  modifications  equation  (12)  is  solved  in 
phase  2. 

With  an  appropriate  equation  of  state  this  form  of  SHELL  was  coded  for  the 
ILL I AC  IV. 


r 


Figure  2.  Rear  of  Slug  Velocity 


4.  SHELL62 

To  program  the  equations  above,  it  was  first  necessary  to  decide  on  a 
method  of  data  storage  as  the  entire  algorithm  for  solution  is  tied  to  this 
choice.  If  one  choses  a  computational  technique,  it  is  necessary  that  the 
correct  data  be  in  the  appropriate  FEM  (processing  element  memory)  or  nearby 
for  efficient  use  of  II.LIAC  IV.  Alternately,  given  a  choice  of  storage  scheme, 
some  algorithms  will  be  much  more  efficient  than  others. 

The  scheme  that  first  comes  to  mind  for  SHELL  is  a  contiguous  mapping  of 
PEs  across  a  row  of  the  SHELL  mesh  as  shown  in  figure  3.  With  this  scheme  one 
can  then  obviously  store  all  the  required  information  for  the  first  column  in 
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PEO  PEI  PE2  PE 3  PE61  PE62  PE63 

PEM  MAP 


Figure  3.  SHELL62  Mapping  of  Cells  and  PEs 

PE  1,  for  the  second  column  in  PE  2,  etc.  Because  of  the  nature  of  the  SHELL 

equations  and  the  explicit  time  differencing  scheme  used,  the  future  values  of 
hydrodynamic  variables  for  any  cell  are  determined  from  values  in  the  cell 
and  its  four  neighbors  to  the  right  and  left  and  above  and  below  only.  These 
data  are  conveniently  available  from  the  right  and  left  by  routing  instructions 
in  GLYPNIR  and  above  and  below  from  different  sections  of  the  corresponding 
cell  PEM.  For  the  initial  implementation  of  SHELL  we  chose  to  limit  the  width 
of  the  mesh  such  that  one  row  of  the  problem  could  be  contained  across  the  set 

of  PEs.  Our  experience  with  SAP  showed  that  it  would  then  be  quite  easy  to 

extend  the  code  to  arbitrary  width.  However,  as  became  apparent  in  coding  SAP, 
including  this  extension  in  the  original  code  would  have  added  an  unnecessary 
load  of  detail  (mostly  in  keeping  track  of  the  mode)  to  the  coding  effort  at  a 
time  when  attention  should  have  been  paid  to  more  significant  problems. 

As  it  was,  writing  SHELL  first  to  only  62  columns  wide  (SHELL62),  then 
extending  it  to  arbitrary  width  (SHELLN) ,  we  discovered  only  two  or  three  errors 
(all  in  handling  the  mode)  in  the  first  code  (SHELL62) .  No  additional  errors 
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were  found  in  the  extension  to  arbitrary  width.  To  handle  the  boundary  condi¬ 
tions  conveniently  (and  simplistically) ,  we  required  that  the  row  width  be 
limited  to  62  cells  or  less  in  this  version.  PE  0  and  PE  63  (or  the  last 
active  cell  +1)  are  then  used  to  store  values  which  will  give  the  appropriate 
boundary  conditions  when  the  row  is  calculated.  To  this  end  we  chose  to  insert 
precalculated  values  in  these  PEs. 

The  implementation  of  the  SHELL  code  in  GLYPNIR  for  62  zone  width  is  shown 
in  appendix  IX.  A  description  of  the  general  features  of  the  code  follows. 

This  was  intended  to  be  a  code  complete  enough  to  execute  on  the  ILLIAC  IV 
simulator  on  the  B5500  computer,  and  this  was  in  fact  accomplished  successfully 
in  January  1971. 

We  have  no  intention  of  expanding  the  code  on  a  line  by  line  bases  here 
because  we  feel  that  little  could  be  gained  by  this  exercise.  But  we  do  feel 
that  a  certain  number  of  comments  is  in  order. 

The  FORTRAN  version  of  SHELL  proceeds  through  the  calculational  mesh  by 
columns — bottom  to  top,  then  left  to  right.  The  GLYPNIR  version  processes  the 
mesh  by  rows— bottom  to  top  (left  to  right  in  SHELLN) .  This  is  a  minor  change, 
but  it  makes  the  code  much  easier  to  handle  in  parellel.  Because  FORTRAN  SHELL 
updates  its  arrays  in  on  top  of  themselves,  some  of  the  unupdated  quantities 
must  be  saved  for  use  in  updating  the  next  cell.  This  logic  in  the  code  has 
been  extensively  changed  by  processing  a  row  of  cells  at  a  time;  in  some  cases 
a  whole  row  must  be  saved  and  in  others  none  at  all.  Further,  certain  geomet¬ 
rical  constructs  that  involve  X  and  DX  as  well  as  the  mode  are  calculated  once 
in  SETDXDYETC  rather  than  continually  calculating  them  throughout  the  code. 

Another  interesting  point  can  be  noted  in  ES,  the  air  equation  of  state  used 
in  this  code.  In  the  FORTRAN  version  of  the  code  the  arguments  of  the  exponen¬ 
tials  are  extensively  checked,  basically  for  speed.  By  checking  the  argument, 
we  can  often  set  the  result  to  0  or  1  and  avoid  the  loss  of  time  spent  in  calcu¬ 
lating  the  exponential.  Since  the  equation  of  state  is  called  at  least  once 
per  zone  per  cycle,  even  a  small  saving  in  time  is  worthwhile.  In  the 
GLYPNIR  version,  however,  the  shoe  seems  to  be  on  the  other  foot.  Here  we  can 
expect  that  of  the  64  zones  being  calculated  at  once,  one  will  require  us  to 
calculate  an  exponential.  Since  this  is  the  case,  it  takes  no  longer  to  calcu¬ 
late  64  exponentials  than  it  does  to  calculate  one,  so  we  have  eliminated  the 
argument  checking  as  excess  overhead. 
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The  subroutine  ERROR  references  the  file  "LINE."  This  is  a  construct  of 
the  simulator,  not  of  the  ILLIAC  itself  and  is  our  only  deviation  from  the  use 
of  constructs  that  will  be  available  in  the  real  ILLIAC  system.  (Of  course, 
with  the  coming  of  the  real  operating  system  for  the  ILLIAC,  the  SIMREADs  and 
SIMWRITEs  will  become  simply  READs  and  WRITES.)  This  is  used  because  communi¬ 
cation  with  0S4,  the  ILLIAC  operating  system,  had  not  been  defined  or  imple¬ 
mented  at  the  time  this  code  was  written.  When  this  implementation  is  complete, 
the  SIMULATE  (LINE,  -  )  will  be  replaced  by  an  appropriate  call  to  the  system  to 
return  our  error  code.  (See  SHELL/OF/THE/FUTURE  listing  in  appendix  XI.)  As 
noted  before,  certain  constructs  that  are  familiar  to  the  FORTRAN  programmer 
either  have  not  been  implemented  in  GLYPNIR  or  they  exist  in  a  rather  primitive 
form.  Various  things  can  be  cited  that  fall  into  this  category,  such  as  the 
partial  implementation  of  the  FORTRAN  EQUIVALENCE  facility  and  the  rather  primi¬ 
tive  I/O  facility. 

None  of  these  lacks  is  really  serious  because  each  can  be  programmed  around. 
The  I/O  Is  probably  least  serious  because  we  will  probably  want  to  do  unbuffered 
I/O  with  our  large  data  arrays.  The  ability  to  equivalence  CU  variables  into 
PE  variables  would  greatly  simplify  our  INPUT  and  OUTPUT  routines.  As  it  is, 
we  merely  use  replacement  statements  and  the  GRABONE  and  PUTCPR  functions  to 
pack  and  unpack  our  PREAL  variable  Z.  Our  only  use  of  the  equivalence  facility 
as  it  currently  exists  is  found  as  the  seventh  line  of  the  subroutine  ES,  where 
we  equivalence  RHO  and  CS. 

One  ray  now  look  at  the  code  as  a  whole  and  question  the  interaction  between 
the  storage  scheme  and  the  code  of  the  routine.  Surely  the  storage  that  we  are 
using,  storing  the  grid  by  rows,  is  the  simplest  conceptually  that  we  can  hope 
to  achieve  (figure  3),  so  it  is  of  interest  to  discuss  the  reasons  that  we,  or 
others,  would  go  to  some  other  method  of  storage. 

In  other  types  of  calculations  one  is  given  an  array  of  some  predetermined 
size  to  operate  on,  and  it  is  desirable  to  be  able  to  come  close  to  this  size 
in  one  storage  scheme.  For  instance,  in  the  code  just  described  with  storage 
by  rows,  a  grid  62  zones  wide  (or  n  x  64  -  2,  In  general)  would  be  most  effi¬ 
cient  in  terms  of  storage  and  PE  usage.  If  we  were  given  a  mesh  instead  of 
being  able  to  choose  one,  we  could  run  into  serious  efficiency  problems.  We 
would  want  to  choose  a  more  compact  storage  scheme  (say  8x8  rows)  and  reduce 
this  waste.  We  do  not  have  this  problem  with  hydrocodes  because  we  will  be 
working  on  meshes  considerably  larger  than  64  in  a  given  direction  and  we  can 
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always  use  the  extra  rows  In  a  given  direction  to  secure  more  resolution.  The 
reasons  for  considering-  such  storage  would  be  with  two  ends  In  view.  First, 
the  least  efficient  part  of  the  SHGLL62  code  comes  In  presetting  the  boundary 
conditions  at  the  axis  and  the  right  of  the  grid  in  preparation  for  the  actual 
hydro  calculation.  This  Is  done  with  essentially  only  one  PG  turned  on.  A 
more  complicated  storage  scheme  (say  8  x  8  or  C-skew)  would  allow  us  to  process 
more  than  one  boundary  condition  at  a  time  because  the  boundary  cells  would  not 
all  be  In  the  same  PG.  This  would  Increase  the  efficiency  of  that  part  of  the 
code.  For  instance,  using  8  x  8s,  we  could  process  eight  cells  at  a  time.  With 
a  7  x  9  C-skew  we  could  effect  100  percent  PG  efficiency  by  updating  64  boundary 
conditions  at  a  time. 

This  more  compact,  blocked  storage  would  also  help  In  calculations  of  GS 
and  opacities  of  materials  because  one  could  expect  a  compact  region  would  have 
essentially  the  same  conditions  throughout.  Hence,  the  number  of  Iterations 
for  each  cell  would  be  approximately  the  same  or  the  part  of  a  table  searched 
would  be  small. 

For  these  reasons  we  Investigated  reprogramming  part  of  SHGLL  using  7x9 
C-skew  storage.  Although  the  boundary  conditions  would  be  performed  at  higher 
efficiency,  the  overhead  In  doing  just  the  normal  hydro  calculations  would  more 
than  outweigh  the  advantages. 

Time  was  not  available  to  Investigate  the  less  than  optimal  8x8  case.  Nor 
was  there  time  to  Investigate  the  possibility  of  transferring  the  arrays  before 
and  after  certain  operations.  These  are  not  currently  thought  to  hold  any 
probability  of  a  big  payoff,  but  certainly  should  be  considered  In  the  future. 

5.  SHGLLN 

SHGLLN  Is  an  expansion  of  SHGLL62  to  handle  rows  of  any  width.  The  principle 
difference  in  the  codes  Is  In  the  handling  of  the  row  storage.  Figure  4  shows 
the  mapping  of  the  SHGLL  grid  onto  the  processing  element  memories.  As  in 
SHGLL62,  one  dummy  cell  Is  required  at  each  end  of  the  row  to  handle  the  pre¬ 
stored  boundary  conditions.  Because  we  are  free  to  choose  the  number  of  cells 
In  a  row  in  a  typical  hydro  calculation,  we  will  arrange  to  have  only  two  cells 
per  problem  row  that  do  not  represent  hydro  cells. 

SHGLLN  was  written  to  handle  a  problem  of  any  size  provided  it  could  be 
memory  contained.  A  listing  of  the  code  is  given  In  appendix  X.  The  code 
shown  compiled  on  the  B6500  at  the  University  of  Illinois  In  approximately  12 
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minutes  of  processor  time  (30  minutes  real  time  but  with  other  users  in  the 
mix) .  One  cycle  was  then  executed  on  the  ILLIAC  IV  simulator  for  a  problem 
mesh  5  rows  tall  and  70  columns  wide  in  30  minutes  of  processor  time  (60  minutes 
real  time  with  no  one  else  in  the  mix) .  The  answers  agree  precisely  with  those 
obtained  from  the  FORTRAN  version  of  SHELL.  One  may  notice  chat  PH3  and  REZONE 
are  included  in  the  code;  these  sections  compiled  but  were  not  executed  because 
of  time  limitations. 

Figure  5  shows  a  sample  line  of  code  from  SHELL62  and  SHELLN  demonstrating 
the  principle  difference  between  them  and  the  solution  of  a  problem  unique  to 
parallel  processing.  The  line  shown  in  figure  5  is  quite  relevant  to  the  code 
and  constitutes  the  computation  of  the  acceleration  of  the  fluid  in  the  radial 
direction  for  an  entire  row  of  the  problem.  This  is  the  finite  difference  form 
of  the  radial  component  of  one  of  the  three  partial  differential  equations  that 
SHELL  seeks  to  solve.  A  detailed  explanation  of  the  code  line  follows. 

MODE  is  a  boolean  variable  (single  word)  that  specifies  bit  by  bit  which 
PEs  are  on  and  which  are  off  for  most  computational  purposes.  For  a  SHELL62 
mesh  as  shown  in  figure  3,  the  appropriate  MODE  pattern  would  be  Oil  ...  110; 
that  is,  the  first  (zeroth)  and  last  (63rd)  PEs  are  off  and  all  the  rest  are  on. 
This  boolean  value  is  stored  in  MI  and  used  throughout  the  code. 

The  variable  U  is  a  PREAD  vector  dimensioned  to  100;  that  is,  101  storage 
locations  are  reserved  for  U  in  each  PEM  (processing  element  memory) ,  and  it  is 
type  REAL  (as  opposed  to  INTEGER  or  BOOLEAN).  The  variable  K  is  type  CINT, 
which  means  it  has  a  single  integer  value.  The  construct  U[K]  then  refers  to 
the  (K+l)th  value  (because  we  start  counting  from  zero,  naturally)  of  U,  which 
in  general  will  be  different  in  each  PEM.  TAU  is  an  area  term  and  has  a  differ¬ 
ent  value  in  each  PEM.  DY  is  single  valued  (for  each  J  which  is  also  single 
valued)  and  represents  cell  height. 

We  now  come  to  an  interesting  construct  of  GLYPNIR,  that  of  a  route.  As 
shown  in  figure  5,  it  is  a  route  to  the  right,  an  amount  1,  of  variable  PR,  the 
pickup  being  under  a  mode  pattern  shifted  1  to  the  left  from  MI.  To  see  what 
happens,  first  shift  the  mode  one  bit  to  the  left  (end  around).  The  pattern 
will  then  be  11  ...  100.  Under  this  mode  pick  up  variable  PR.  This  pickup 
will  be  accomplished  in  PE  0  through  PE  61  and  nothing  will  be  done  in  PE  62 
and  PE  63.  Next,  route  (routing  registers  are  always  on  regardless  of  the  mode) 
the  value  picked  up  one  element  to  the  right.  The  mode  is  returned  to  its 
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SHELL62 : 

MODE=MI 

U[K]=U[K]+TAU*DY[J] * (RTR(1 , ,PR) -PR) *DT/ (DX*AMX[K] ) ; 

WITH  TYPES: 

PREAL  VECTOR  U,  AMXflOO]; 

CINT  K,  J; 

PREAL  TAU,  PR,  DX; 

CREAL  VECTOR  DY[ 100] ; 

CRRAL  DT; 


SHELLN: 

LOOP  1-0,1, IBLK  DO  BEGIN 
MODE* IMODE [ I ] ; 

U[K]=U[K]+TAU[I]*DY[j]*(RTR(l,,PR[I+IR])-PR[I])*DT/ 
(DX[ I ]*AMX[K] ) ; 

END; 

Where  IR  is  type  PINT  and  *  0,0,0,...,-! 


Figure  5.  Sample  Line  of  GLYPNIR  Code 
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original  value  (MI)  and  che  routed  values  are  available  for  computation  in  the 
PEs  to  which  they  have  been  moved.  Now  since  PR  is  type  PREAL,  there  can  be  a 
different  value  in  each  PEM.  As  can  be  seen  in  appendix  X,  a  pressure  term  is 
actually  stored  in  PR  so  that  when  PR  is  subtracted  from  RTR  (1,,PR)  we  actually 
take  the  difference  in  pressure  between  each  cell  and  its  neighbor  to  the  left. 
Of  particular  concern  is  the  value  of  PR  in  PE  0.  Although  this  PEM  does  not 
represent  an  actual  hydro  cell,  a  value  of  PR  must  be  available  for  the  route. 

An  appropriate  boundary  value  is  stored  in  PEO  before  this  line  of  code  is 
reached.  The  remaining  variables  in  the  line  represent  the  time  step,  cell 
width,  and  cell  mass. 

Figure  5  also  shows  the  same  computation  for  SHELLN  and  the  changes  required 
to  accomplish  it  for  any  row  width.  The  single  statement  is  replaced  with  a 
loop  over  the  number  of  PEM  rows  required  to  store  one  row  of  the  actual  hydro 
mesh.  We  now  have  an  appropriate  mode,  stored  in  IMODE,  for  each  PEM  row.  The 
value  for  IMODE[0]  would  be  Oil  ...  1  for  the  mesh  of  figure  4,  while  IM0DE[1] 
would  be  all  on  and  IM0DE[2]  would  be  11100  ...  0.  For  1*0,  the  line  of 
SHELLN  code  executes  similar  to  that  for  SHELL62  above.  Now  since  the  area  and 
cell  width  are  different,  in  general,  for  each  cell,  TAU  and  DX  have  been 
expanded  to  vectors. 

These  usages  are  simple  and  straightforward.  The  significant  development 
comes  in  the  handling  of  PR.  First  of  all,  it  is  expanded  to  a  vector  like  TAU 
and  DX  and  then,  to  do  the  end  around  route  correctly  (that  is,  to  connect  cell 
64,1  with  63,1  (and  not  127,1)),  it  is  necessary  to  subtract  1  from  the  subscript 
on  PR  in  PE  63  only.  This  is  accomplished  by  adding  IR,  which  is  equal  to 
0,0  ...,  0,  -1,  to  1  after  the  mode  is  shifted  end  around  one  to  the  left, 
similar  to  that  described  for  SHELL62  above.  This  construct  simply  and  effec¬ 
tively  connects  the  three  PE  rows  to  form  one  problem  row  for  all  routes  to  the 

right.  For  routes  to  the  left,  a  similar  construct,  employing  IL  ■  1,0,0 . 0 

in  place  of  IR,  is  used. 

6.  SHELL/OF/THE/ FUTURE 

Although  SHELLN  completed  the  development  of  a  code  similar  to  the  regular 
production  SHELL  code  used  at  AFWL,  the  fact  that  it  is  memory  contained 
restricts  the  prob.t  .m  size  to  about  20,000  cells,  the  same  number  that  can  be 
run  on  our  CDC  6600.  We  therefore  expanded  the  code  once  more  to  use  the 
1LL1AC  IV  disk  for  problem  storage.  The  result  was  SHELL/OF/THE /FUTURE,  which 
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solves  the  problem  by  blocks, reading  and  writing  them  to  disk  In  a  fashion 
similar  to  the  use  of  extended  core  at  AFWL.  A  listing  of  this  code  Is  given 
In  appendix  XI.  Programs  to  write  and  read  the  ILLIAC  IV  disk  files  are  given 
In  appendixes  XII  and  XIII.  These  correspond  to  the  programs  I  and  II  for  SAP 
discussed  earlier. 

SHELL/OF/THE/FUTURE  will  be  able  to  use  the  full  capabilities  of  the  ILLIAC 
IV  In  solving  a  two-dimensional  hydrodynamics  problem  with  a  single  material. 
However,  since  this  code  was  developed,  we  have  written  a  new  hydrodynamics 
code  at  AFWL  that  we  will  probably  want  to  code  in  GLYPNIR  and  use  on  the 
ILLIAC  IV  in  place  of  SHELL/OF/THE/FUTURE.  This  new  code  will  include  the 
ability  to  handle  any  number  of  materials  as  well  as  containing  other  features, 
and  it  will  employ  new  differencing  techniques. 
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SECTION  IV 
CONCLUSIONS 

The  basic  purpose  of  this  exercise  was  to  translate  working  one-  and  two- 
dimensional  hydrodynamic  computer  programs  from  serial  logic  applicable  to  a 
machine  such  as  the  CDC  6600  to  parallel  logic  appropriate  for  the  ILLIAC  IV. 

To  accomplish  this  we  found  it  most  convenient  to  reduce  the  problem  to  its 
basic  elements  and  start  from  the  beginning.  Even  with  the  more  adventuresome 
effort,  the  SHELL  code,  we  were  able  to  make  rapid  progress  after  rederiving 
the  basic  difference  equations. 

We  expected  to  encounter  difficulties  in  using  the  ALGOL-like  GLYPNIR 
language  but  actually  found  after  working  with  it  for  several  weeks  that  we 
were  able  to  think  in  the  language  with  relative  ease.  The  real  difficulties 
that  we  found  in  the  code  translation  effort  were  all  related  to  the  parallel 
structure  of  the  ILLIAC  IV. 

The  SHELL  code  as  it  exists  at  AFWL  consists  of  a  few  simple  routines  which 
solve  the  hydrodynamic  equations  coupled  with  a  fairly  large  number  of  support 
routines  which  make  a  working  and  useful  code.  Included  in  this  group  are  the 
massless  (for  some  applications  massive)  particle  transport  routine  generally 
used  to  follow  Lagrangian  interfaces  and  an  automatic  rezone  routine.  For 
problems  with  several  materials  a  diffusion  limiter  is  available. 

The  coding  in  GLYPNIR  of  the  basic  hydrodynamic  equations  was  relatively 
straightforward.  It  was  not  as  simple  as  coding  for  a  serial  machine,  but  we 
were  never  completely  baffled  in  our  search  for  efficient  parallel  algorithms. 

The  support  routines  were  more  challenging  for  certain  of  the  authors.  We 
found  that  it  was  not  a  trivial  exercise  to  generate  an  efficient  particle 
transport  routine.  Likewise,  the  automatic  rezone  routine  presented  certain 
difficulties  in  handling  the  rezone  in  the  radial  direction.  All  of  these 
problems  were  eventually  overcome  and  satisfactory  algorithms  were  developed. 

We  observed  a  B6500  in  action  at  the  University  of  Illinois  as  it  compiled 
GLYPNIR  programs  into  ILLIAC  machine  language.  For  our  SHELL  code  the  overall 
time  turned  out  to  be  a  little  better  than  100  cards  per  minute  of  processor 
time.  This  is  for  the  complete  cycle  from  GLYPNIR  source  code  through  ASK  and 
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the  collector  and  loader.  The  SHELL  code  consisted  of  about  1200  GLYPNIR  cards 
and  took  approximately  12  minutes  of  B6500  processor  time.  The  speed  of  a 
B6500,  we  understand,  can  be  greatly  Increased  by  enlarging  its  memory  size. 

We  hope  that  the  B6500  or  similar  computer  that  drives  the  ILLIAC  IV  will  be 
equipped  with  a  large  memory. 

We  anticipate  that  a  running  SHELL  code  will  produce  too  much  data  to  be 
stored  on  magnetic  tape,  but  the  laser  store  as  advertised  should  handle  the 
problem  quite  nicely.  We  estimate  that  approximately  8  x  109  bits  of  data  will 
be  sufficient  to  describe  the  complete  time  history  of  a  500,000-zone  problem. 
To  process  and  analyze  this  volume  of  data,  analysis  routines  will  be  written 
in  GLYPNIR  for  execution  on  ILLIAC. 
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APPENDIX  I 

COMPUTER  LISTING  OF  THE  BASIC  STRIPPED  VERSION  OF  SAP  IN  FORTRAN 
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1  IeXa  1  IX  A0*tFbAi3XAlbHAl  lX,3HhKlA)  3XAlEbAl3XAlhTAl2XAiE./FiA3XAlhjA  ‘Je  >1 

2  /  1  h  )  i .  £  1 1. 

EbFEATllH  .I3AbE14,bAl«> 

LE«b  bj£i/ 


SE  UM.  NT 


E 

E 

E 

Ei 

F 

E 

F 

E 

E 

E 

Ei 

h 

h 

h 

h 

h 

♦■ 

h 

h 

h 

h 

h 

h 

h 

h 

h 

h 

K 

r 

h 


APPENDIX  II 

COMPUTER  LISTING  OF  ‘FORTRAN  SAP 


CQMh0N/MAlN/Z(64),U(6«),P(6<l)»Q(6«)»RHD(64)»X(64)>7H(64)» 

1  EPG(64),CS(64),T(64),PnRT(64) 

C0MM0N/TEMP/RH0N(64)#nVf 64)#CRNT(64)#GM0NE(64)#P2(64)#E1(64), 
lPI(64),RHDf(64)»DTZj(<>4) 

EMJIVAIENCF  (PHOB#Z(l  ))#(CN,Z(2))#(TIME#Z(3))#(DFLT#Z(4))# 

1(U VC»Z(5))#(JSTAR#Z(6))#(JSTP1»Z(7 ) ) » ( MJSTAR# Z( 0) ) # (MJSTPl » Z( 9) )» 
2(JF IN»7(l0))»(PMAX»Z(ll ) ) H JPMAX »Z(1 2 ) ) # ( XPMAX »Z ( 1 3  )  )  » (C02#Z(1 4 ))# 
3(Cl,Z(15))HTCl#Z(l6))#(FCD2#Zd7)),(0TZJM#Z<  18)) 

FILE  I NPUT ( 0#  704 )  SERI AL 
FILE  OUTPUT ( 0# 704 )  SERIAL 
BINARY  NJST ARC64 ) »  MJSTP1 (64) 

23  DELT«OTC 
DTC«DT?JM 
MJSTAR.(0#63)»0 
MJSTAR.(1# JSTAR)«1 
DU  341  POOE-MJSTAR 

341  U(*)«U(*)*DELT*((P(0*0(*)"0(**l  >"P(**1  )))/  ( RHO(  *  )*(X(  * )-X(  *-l ) )  ♦ 
lRH0(**l)*(X(**l)“X(O)) 

00  351  MOOE-MjSTAR 
IF (ABS(U(O).LT.UMIN)  U(*)«0. 

351  X(*)«X(*)*0TC*U(O 

IF (U( JSTARJ.GT.O. )  JSTAR«JSTAR*1 

IF (JSTAH.GE. JFIN)  STOP  351 

JSTP1bJSTAR*1 

MJSTPl. C0#63)«0 

MJSTP1.(1,JSTP1)«1 

DL  471  M0DE«MJSTP1 

RF*0N(O*ZM(*)/(X(O**3"X(*“l )**3) 

OV(*)«J ,/RHON(*)-1./RhO(*) 

0U(*)«U(*)-U(*-1  ) 

0(*)»0. 

IF (OV(*).LE.O,)  0(*)«RHnN(*)*(AMINl(0U(*)#0,)*C02-Cl*CS(*))*AMINl( 
1DOCO.O.) 

IF (0(*).LF.OMIN)0(*)B0. 

CRNT (*)«0. 

IF  (QCO.NE.O. )  CRNT(*)«TC1*CS(*)-FC02*DU(*) 

GHDNEf  *  )«0 .4 

P2(*)bGM0NE(*)*HH0N(*)*FPG(*) 

El(*)«AMAXlt<EPG(*)-0V(*)*P(*))»0.) 

PI ( * )bGMOnE( * )*RH0N(*)*E1 (* ) 

EPG(*)bAMAX1((EPG(*)-(P2(*)*P(*)*0(*)*0(*))*0V(*)/(2.*(f2(*)-P1(*) 

1)/P(*)))#0.) 

P(*)«GMONEC  *)*RHON(*)*EPG(*) 

P2( *)b«U.0C1 

IF  (DV(*  )  ,f,E  .0, )  P?(  *  )b*P2( * ) 

RH0F(*)«RH0N(*)*( 1 ,*P?(*)) 

PFU0GF(*)»GM0NE(*)*RH0F(*)*FP6(*) 
CS(*)bSQRT((P<O-PFU0GE(*))/(RHDN(*)*RHDF(*))) 
0TZj(O«NW*(X(*)-X(*-l))/  (CRNT(*)*1 .E"2«CS(*)) 

IF (DTZJ(J).LT.DTZJM  )  DTZJMbOTZJT J ) 

RF*0(*)«RH0N(*) 

PMAX«AMAX1(P(J)#PMAX) 

471  CONTINUE 

TIME«TIME*DTC 
C  EXIT  CONTROL 

WRITE  (OUTPUT)  Z 
GO  TO  23 
STOP  7777 
END 
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APPENDIX  III 

COMPUTER  LISTING  GLYPNIR  SAP62 


SCPDF  LIST  TIP  pRHGA  SliMPV 
SprOfPT  PROnllCTTON  SAVF 
RFGTN 

PRFAl  SUPPOl'TINF  SORT  AS  RGAfPREAl  y  AS  RGA ) ) 

PTC  TP 

SI  PC  IMS  110  ASKRATCW  14 
$$TAPF?«  GL vPMr/GSORT  SFRlAll 

fppj 

SUBROUT  TPF  PI)TCPR(CRFAL  X.PRFAL  OUT  ?#flNT  T)| 

BFGTN 

FOR  ALL  PfP.I  00  Z«-XJ 
FPPJ 

SUPROUTIPp  PUTf P I ( C I  NT  X»PRF AL  OUT  7,CTPT  TtI 
rfgtn 

cn  RFAL  A I 
A* X  I 

FOR  ALL  PFN«T  00  Z«-A| 

FkDI 

CPF AL  PRPP.TIMF»PFLT»OTC#PMAy»yPPAy»CO?»Cl#TC 1 »FC0?»PTZJM, GPIP#UPTP»WNJ 
CTRT  P. JSTAP* JSTPl , JFTN»7PRPR»7N»ZTImF#7DF|  T»7DTC»7JSTAR»ZJSTP1  I 
TINT  7JFTm,7PMAX»ZJPMax,ZXPPAX,7C0?,ZC1 ,ZTC1 #ZFCO?.70T7JM»7oMTNJ 
CIPT  JPPAV»7I'PTp»7WW»ZNSTPP,NST0PJ 

PRfAL  7»U.P»0»RH0»X»ZM»FPfi»CS»T,RPRT,y S»RH0P»0V»CRPT#GMONf»P?#F1 f Pi  I 
PPFAI  PDU,RHPF.OT7J»DU.BfU0GF» 

LAPFL  START»AGAlP»FlP#STOPl 
*  VARTABLFS  MAVF  BFEP  OFC! ApfP 
»SFT  UP  IPPiiT  APO  OUTPUT  FlIFS 
FIIF  I40lSK?«"SAP"/"BTnl«(!l  ROWS  FULL )  I 
r II  F  T4DISK1»"<AP"/"IT0P«(!1  ROWS  FULL  )  I 
OPFPf  I4DlSK1#0»m 
PPFPf  I40TSKP»0#m 
STARTl 

SIPPFAOCI40TSKP»7)I 

SlPhRlTt(I4PlSPl»7)| 

7PP0B«-0j 
7P«-1 1 
7TIME*? J 
7PFI T*3l 
7DTC«-4I 
7 JSTAP*5l 
7JSTP1+AI 
70PIN«-7| 

7UMIN*P| 

’ JF I I 
7PPAX«-10I 
7JPPAXM1 i 
7XPMAX«-|?| 

7C0ZO3J 
7C1«-14I 
7TC1+15J 
7FrO?«-16J 
7PTZJP  +  17ii 
7PWM8J 
7PST0P*! 9 l 
DT7J4-100000.0I 

OMIP«T,PARPNF(Z*ZOMIW)I 

NW«-GRaP0PF(7#ZNN)I 

0TZJH«-GRAPPPF(Z#Z0TZJM)I 

JSTAR«.f,RApPNF(Z#ZJSTAR)l 

PROP*r.PAROPF(Z»ZPROR)» 
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N*GPAPPNrr7»ZN)l 
TIHF*GPABpNF(Z#ZTIMF)l 
DTC*6PABPNFf Z#ZDTC)I 
O^TN«'nRAPn^FCZ»ZllMTN)> 

JF  IM*r.RARpNF(Z#ZJFIN)  I 
C0?*GPAPPpF(Z#ZCP2)l 
Cl*GR4PtlMF(Z#ZCni 
TC1*GRABP*Ff Z#ZTC1  )l 
NSTPP*GRAppNF(7#7NSTOP)! 

FCn2«-fiP»RPAT(Z»ZFCn?)i 
AGAIN!  OFl?*DTC  I 
DTODT7JM» 

MPPF*fPFN<jSTAP)ANPfNPT  BOOLE ANf -0) ) I 

U*U*0FI  T*f  (P*0)-RTLf  1  #TBtlE#f  P*0)))/(RHPx(yPTL(-l  #TP||F»  v  )  )*RTI  (  1 
#TPUE  fPHOl*fRTLM#TPUF»X)-XJ)l 
IFF  APSf|i><HMIN)  THFF  |i«-0  •  I 

x*x*oVr*iii 

IFfGRAR0NFfl'#jSTARl>0.l  THFN  JSTAR* JST AP*1  I 
IF  f  JSTAR>,|FIN)  THFN  Gp  TO  STOP! 

JSTP1*JST»P*1  I 

MCDF*SHIFTPf  1#M0PF  PR  RPOLFANf -0) )  I 
X3*X*y*XI 

RH0N*7M/fy3"RTLf"l#TRUF»X3))l 
DV  +  (Rpn-RHnN)/(RNON«Rt'n)  J 
PU*fl'-PTl  f  *1  #TRllF#U)  )  I 
0*0.1 

PCU*Otil 

IFfDU>P.Ol  THFN  POlM-0 .0 « 

lFfDV<0.)  THFN  0*RM0NxfPDUxCPP*fU  ('^JNPPUI 

lF(0«PTNi  THEN  0*0.1 

cpnt*o.o» 

If (0>p. 01  THFN  CRNT*Tf 1*CS-FCP2xPDH| 

GPPNF*0.#i 
P2*GMPNFXRHPN*EPGI 
F1*(FPG*PV*P)I 
IFfFlcO.)  THFN  F1*0.I 
P1*GMPNF*RHPNXF1 I 

FPG*FPG»fP?+P*0*0)*0V/f ?. *(P?-P1 )/P)l 

lF(FPr.«0*.  )  THFN  FPG*0*.  I 

P*GNOnF*PhPN*EPGI 

RHOF  *RHPNxl .0011 

IF (DV>O.Oi  THEN  RHPF*PH0N*.999» 

PFllPGF*GHpNFxRHOF*FPGl 
C5*S0pTffP»PFl!)DGF  )  /(RHPN»RHPF)  )  I 
DTZJ**W*f X-RTl(-l*TRUF»X) ) /( CRNT*.01 *CS ! J 
DT7jM«.HlNfPTZJ)l 

rhp*rhpm 

PMX*MAXf  p)l 
TIHE*TlHF*PTCl 

mopf*trufi 

PUTCPRf PRPB#7#ZPR0RJ| 

PUTCPR(DT7JH#Z#ZPT7JMil 
PUTCPP(PMAX#Z#ZPHAX)I 
PUTCPPf TTMF#Z#ZTTMF)I 
PUTCPBfDTf#Z#ZDTC)l 
PUTCPT(N#7#7N)J 
PUTCPTf JSTAR#Z#ZJSTAR)» 

PUTCPPf  DFI  T#Z#ZDFlT)l 
putcpr(xpmax#z#zxpmaxii 
pi'TCPrf  jpmax#z#zjpmaxti 
PUTCPTf J*TP1#Z#ZJSTP1  II 
PUTCPTf NSTPP#Z#ZNSTOP)  I 
SIMWRITEflAniSNI #7)1 
N*N*1 i 

IF (N<wFT0P )  THEN  GP  TP  AGAIN! 

FlNlSTOPlFNn. 
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n  r>  r> 


APPENDIX  IV 

COMPUTER  LISTINGS  OF  PROGRAMS  I,  SAP62  INPUT,  AND  II,  SAP62  OUTPUT 


scarp  list  singi  r 

FILF  lBSAP/PTOT,ltNlTBPlSK,SAVF«1»l PCK#SFRT4l  , ARFAba, 

1BLPCKTNG«1 1  RFC0RDb540.RUFFFRb? 

C0mmPn/MTI'V7C64),IK6a  >»P(64)#0t64),RMn(64),X(6a),7M(6A  )» 

1  FPGr64>,C$(64)»T(64>»P0PTC64> 

peal  pfsct(12) 

FCMIVAI FNpF  C PNUM, 7(1)1,  (N,7t 2) ) # (TIMF » Z( 3) ) , (PFLT , 7f  A )  ) # 
1(DTC,7(5))#f JSTaR,7(6))*C JSTPl,7f 7) ),(<H'IN  ,7(P)>» CUMIN  ,7f 0)), 

2(JFlN.7f  10))#CPMAX.7(11  ))»CJP**A»#7f  l?)i.(XPMAX#7f  n>l#CfP?»7(1*))» 
3(C1»7f 1 5  >  >  *  f TCI »7C 16) ).tFCP2»7(17))#(DT7JM,7( l«)),f  WW»7C19)), 
4(NST0P.7(?O) ),(OFSCT(1  1,7(52)) 


r  READ  sTATfmFNTS  0019 

pfap  Ri#rrr$cTm#T«i#i?)  oo?o 

r  ooi6 

JF  INb*P 

PFAP  P3»  PNl'M, C0»C1,WN,TS,CKF',IIMIN,0MlN,TIMF,B0MAS»Hf)B»ANGlF»PxMUl  00?2 
TF  (IIM’tP.FO.O,  J  (JMTNb  1.0F-10 

PRINT  P«#(prSCT(T)»T*1 .12)  0029 

print  pf.nc .ntfdit, jfin, jcycs»n»itimf»i  tapf,pnum,co,ci ,hw»ts  0031 

PRINT  P7,rKP»UMIN,0MlN.TIMF,BPMAS»H0R»  AK'Gl  F  0032 

XNTNUSsO.  0041 


umeo.o 

X(1  )«0.0 
JSTARb.JF  IN 
I«  0 
J«  2 

DO  16  KbI , JFTN 

C  OX  IS  THF  PFLTA  X  FDR  TMF  70NFSCCR1 

r  f77  IS  TMF  FNEPGY  FDR  THF  7PNFS(fRcS/Gm> 

P77  Is  THF  DENSITY  FOP  THF  7PNFS(GV/CM3^ 

U77  Is  THF  VELOCITY  F PR  THF  7PNFS(PP/SFr ) 

I  IS  THF  MAY  70NF  NIIMPFR  FOR  THFSF  CPNPTTinNS 
I F  C  K  .1  F.li  GO  TO  11 
RFAD  PO,Py»F77,R77,il77»I 
IF (U77.NF.0)  JSTARb  I 
IF(  JSTAP.GT.  JFIN)  .ISTARb  JFIN 
11  X(  J)ByMiN|iS+PX 

IF (F77.GT, 0)  GO  TO  701 
PH0(J)«l".??5F*3 
FPG( Jl«2.F9 
GO  TP  70? 

701  RHP(J>bR77 
FPGI JibF77 
70?  GHONF  »0.4 

P(J)«C,HOMF  *PHO(  J)*FPGC  J) 

CSC J)bSObT  CPC J)*1 ,47PH0( J)  ) 

0(  J)B(1. 

IK  J)b|i77 

7N(J)bPHO(J)*(X{J)**3»YHTNU5**3) 

XMTNUSByC.'l) 

TF(MOrKK»50).F0.1  )  PRINT  2004,  OFSCT.PniiM 

PRINT  ?00B»  K,I)X*X(  Jl.FPGt  J),P(J),RHO(.l),'K  J),T(  Jl,ZM(  J)»K 

JB  j*i 

C 

16  CONTINI'F 
C02«CP*CO 
TC1b?.*C1 
rco?Bi.*ro? 

DTZJMbWV»*XC?)/(S0RT(P(7)*1 . 4/PH0(2)  )  )  /?. 

DTCbDt7JM 
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2cn*(e»iiM 

Z(6)«viSTAR 

7(10)bJFIN 

Z(lA)«f07 

Z(15)sCl 

Z(16J«TC1 

znr)«Fcn? 

znB)«nT7.iM 

Z ( 1 9  )*WW 

NSTOPbIO 

CALL  wPTIAn(7#70A) 

STOP 

61  FORMAT  ( 1 ?A6  ) 

63  FORMAT  (6f10.0> 

6A  FORMAT  (1h1»1?A6)  _ 

66  FORMAT  ( 7rM0SYMRnL.3».l0HPFF!MITl0N#66v,5HWALUF/3HOMC#FX»36MNl'MBFP 
1  OF  DlFFFRFNT  JPLTR  AMD  NOUMPS #40 X# I 6/7M  NTFDTT,3X#30HNUMBFR  Of  DI 
2FF  FREnT  FnlT  tlMFS#46y.I6/5w  JFlN»6y# 26MNUMRER  OF  7PNFS  TN  PROBLEM 
3#50X#T6/6m  JCVCS#  Ay#3PMMAXTM(lM  WllMRER  OF  PROBl  FM  CVOLFS ' 44X, I6/?H 
AN#  6y  »  ?6HC YCL  E  MlMBFR  OF  INPUT  DATA#50X» T6/SH  ITIME.AX.76HPOWFP  OF 
5TEN  FOP  MRST  AUinT»46X»I«/6H  ITAPF»4X»41hmSED  TO  OFsTGNATF  TMPUT 
6T APES  FPR  RFStART#35X»T6/5H  PNUM,5v#14HPR0Pt  Em  MUMRFR,aOY#F6.3/3H 
7CO»7x.PAHol'AORATTC  VISCOSITY  TERM»A6x»1PE1P.3/3m  (M  »7y»  ?1  HI  TNFaR  V 
6 1 SCOST TV  tFRM,A9X#F12.V3H  WW#  7V#  3SHMUI  TIPI  TFS  ThF  F  Al  f  liLATFO  TIMF 
9  STFP.35v.F12*.3/3H  TS.7v#20HTIMF  Or  PROrLFM  STOP»50y»f1  7 . 3  Y 

67  FORMAT  f  Ap  CkP# 6X # A 1  HOT  ST  A wf F  OF  MAX.  PRESSURF  AT  PRORlEM  ST0P.29X 
1 # 1PE17.3/6M  UMIM,5y,?7HMfNIMUN  VELOCITY  IN  PROPI  FM,a3y#F17.3/5h  OM 
2IN#5X.70HmINIMUN  0  TN  pRUPLFM#50y # F 1 2. 3/ 

3  6M  TIMF**iX,2AuTTMF  OF  START  OF  PROP!  EM, A6X»E1 

42.3/  6H 

6  P0MrtS#AY#l?HMASS  OF  mPF#56X#F1 7. 3/ 

74H  H0p#6V#77HHFI6HT  OF  7FRO  POINT  IN  CM . * 43y »F 1 2 ,3/fH  ANOLF# 4 V# 36H 
BDlRECTlON  OF  RUN  (OFGPFFS»POS.  IIPMARO>»32x»F12.3> 

69  FORMAT  (Ar!2.4M5)  .  „  . 

2004  F0RMATT1M1»1?A6/1H  /5y.7HPNUM  s  #F10.2/1H  73X# 1 M J» 1 ?X#?h0X» 1 3 V# 

1  IhX.11X#3HFP6»13v#1mP*11X#3mRmO»13X#1hU»13X#1MT#12X»2m7M#3X#1HJ# 

2  /1H  ) 

2005  FORMAtMM  » I 3#6E1 4 .5*  T  A 1 
FNO 

SURROllT INF  WRTI40(7»NY 
RE  Al.  aT5A0)»7(N) 

PRINT  99 

99  FORMATMMll 
LI  *  0 

NN*  N 
5  K  *  MM 

TF(K.fiT.2S6)  K«?56 
00  11  T»1.M 
L«  U+t 

IF  (L  ,r.F,5?.  AND.L  .LF.63T  60  TO  50 
CALL  PMAFf7Cn»A(2*I-n»A(?*T>> 

60  TO  10 

50  CALL  RI4Tf7fL)*AC7*I*l)»Af?*m 

10  IF (7(1  l.ro.O.)  00  TO  11 

PRINT  100#  l #Z(L  1»7(Ll»A(?*I*l T#A(?*I) 

11  CONTINUE 

100  FORMAT(I10.F30.6.E?0.6.?020) 

write ( i >  ( a( n# i«i »54(ii 
nn«  nn-k 

if(nn'.fc'.o)  return 

LL«  LI +256 


027A 
0276 
0  777 
0760 
0761 
0767 
0763 
076A 
0765 
0766 
0767 
0766 
0769 
0790 
0791 
0792 
0793 
079A 
0795 
0796 
0797 
0799 
0770 
0771 
0772 
0773 
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60  TO  * 

r  kd 

SUPROliTINF  R5I *F { A  #P#  C  J 
FCIITV»I  FNcF  ( XX# 1 1 ) 

C 

XE»COnrATrO# A#«1#?#7J 
IFOF'.Fc'.OT  60  TO  10 
CALL  PT4FfA»B#C) 

RFTURN 

r 

10  xx>  A 
Xe  TT 

x«  x+n.i 
x*  x-o.i 

CALL  PT«Ffy#R»Ci 

RETURN 

FKH 

SUPROliTINF  R!4F(A#P»C1 
FCl'TVilFNrF  CXF»IF8) 

DATA  TPTAs/OAOOOO/ 

IFIA.NF.OT  60  TO  10 
R«  0 
C«  0 
RETUPN 
C 

10  XS»CONPATfO#A#47»2#l ) 

Xf ■C0NrATf0#A#42#3»6) 

IF(XS'.NF'.O)  IFB  —  TFR 

ieb«mfb 

B9«CnurATfO»A#A7#9#n 
B10«CPNCAT(0# A»A7»10» 1  1 
TFBF«? 

IF(PIO.NF.O)  1FBF»1 
lFfP9*.NF  .0)  TEBFbO 
ieomfbmfbf 

IF»  IPIAS-IF09 

BB«C0wfATr0»A»16#l#n 

BB«C0NCATfPB»IE»17#33#l5) 

B«C0NrATfRR»A»32#9*lEflF#16) 

C«CONfATfO> A#16#25*IE8F#23-IEBF) 

RETURN 

FKO 

SUPROliTINF  Bl«l(A#P#Ci 
BB«COnC AT (  0#A#1A#1#1) 

P*COKTATf RP» A»A1#9#7) 

C»CnNrATfO#A»16#l6#32i 

RETURN 

ENO 

SAP  STRIPPFO  ANO  READY  FOR  ACTION 

1*.  1  ,8  0.5  0,2  30.  1.00  0.  0.0 

•0000  0,  0,  0*  0, 

r.  2, FI  3  ,  1  ,F*3  6.F05  30 

l‘.  o'  0*.  0.  63 


m 


ACARO  LIST  SU’Gir 

r  ILF  2«SAP/IT0P#IINITbDISK»SFRIAL»IIL0CKInGb1,pECORD«540 

C0MM0N/MATN/Z(64).|l(64).Pf64).0(64)«RHnf 64).Xf64).7Mf 64 ). 

1  FPGf64),CS(64)#Tf64>,R0PT(64) 

REAL  pFSC  T(12) 

FOUIVlIFNCF  (PNI)M»7fl  })»  fN»Zf2))#fTIMF»Zf3))»(PFl.T»Zf4))» 

If DTC#7f5))#f  JSTAR.7f6)).(  JSTP1 # 7f 7) ) # ( OMIN  #7(6)),fUMIN  ,7f9))» 

2f  JFIN,7M0))#(PMAX,Z(11 D.fjPMAY.Zf 12)),(XPMAX#?fl1)),fro2,Zf1 4)), 
3fCl#Zf 15))#fTCl.Zfl6)).fFC02#7fl7))#(0T7JM#?(l6)>,f  WW»7f 19)1, 

C 

4frrscrf  n,7f52)) 

IPO  CALL  PFDT4Pf7#704) 

PRUT  094 
999  FORMATflHl) 

DO  199  T  *  1 #704 
IFf7(T).Fo.O.)  GO  TO  199 
196  PRUT  200,  1*7(11*7(11. 7(1) 

199  continuf 

200  FORMATfI10»F20.6,F30,6»n40) 

N«Zf?i 

JSTAR.7f6) 

DO  601  JJ.l.JSTAR 
J«  JJ*1 

IFfMOPf  JJ  ,50),NE‘.1  )  GO  TO  601 
PRINT  M»fOF$CT(I)»I»i»l2) 


PRTNT  96.PNIIM.N.TIMF.0TC  0179 

601  PRINT  99. jj.X(J)»U(J).0(J)»P(J)#FPG(J).fSf J).RHO(J),  T(J)»JJ 
GO  TP  100 

f  0259 

C  FORMATS  FOR  SAP  0?60 

64  FORMAT  (1H1.12A6) 

98  FORMAT  doHOPROB.  NO ,F« . 3» 3X» 2HNbi4. 3X.5HT tMF«F1 4 .5, 3X.4HDTC «rl 4 .5  0316 

1  /4 HO  J.7y.lHy.1?X.lHl'.12X.lHCl.l?X.lHP.l?X.lHF.10x.4H  CS  »10X»3HRH  0317 

20.9X.6HP/PP-1.6X.1HJ/)  0318 

99  FORMAT  (1H  14# 1 P8F1 3.5. 1 4 )  0S!9 

FNO  0320 


SUPROllTINF  RF0I4D(7#N) 

RFAL  77Z( 1 6?0)#7(N) 

LOGICAI  FIRST 
DATA  FIRST  /.TRUF./ 

C 

IF(FIRRT)  PFWIND  2 
FlRSTi  .FALSF, 

REAP  f  2»FnP»200)  (77ZfTl.lBl.540) 

REAOf ?)  f7ZZf 14540). 1.1 .540) 

RFADf?)  f777f I41060)#Tb1,540) 

00  10  IbI ,256 

IF  f  I.f,F.5?.AND.I.LF»63)  GO  TO  5P 
CALL  T4BFf77Z(2M"l ),777f2*T),7(I )) 

GO  TO  10 

50  CALI  I«BTf  777f2*T-l  )#7Z7(2*I)#7f  ID 

10  CONTINIIF 

DO  11  Ial.256 

11  CALL  T4BFf7ZZ(?*T4539),Z7Zf2*T4540)#7(T4256)) 

DO  12  1.1. 192 

12  CALL  T4RFf7ZZ( 2* 1410791.7 ZZf 2* I*1060)»7f 14512)) 
RETURN 

C 

200  PRINT  201 

201  FORMATf 24H1END  OF  FILF  ON  14  DISK.) 


STOP 

ENO 

SURROilTINF  1 ABF  v  *y Pj»C  ) 

INTEGFR  OlOO 

data  ipi a$#pioo/oaoooo»oioo/ 

EOlItVAl  FNfFfXE*  IE) 

C 

C 

XE  ■  CnK'C»T{0#A#33#1T#15) 

IE*  IftIAS*39-IE 

c  iep  is  power  nr  8 

10  TE8*  TF/3 

c  ifpe  ts  rt  shift  of  MANTISSA  AT  IF6 

IE8E  •  TF-1ME8 
If  (IERF.Of.O)  r.o  TO  20 
!E8F»  TE0F+3 
IEB*  tF8»l 
20  IF8*  -JF8, 

IF(IFS.IT.O)  IE8*  M00-IE8 
C  PlrK  HP  SIGN 

OUT*  r0NCAT(0#A#l»16#n 
C  PICK  UP  FXPONFNT 

OUT*  rnNCAT(0UT#ie«#2.«!»7) 
t  PICK  UP  TOP  OF  MANTISSA 

OUT*  rONCATIOUT# A»9*IF«F»3?» 18) 
f  PIf K  IIP  POTTOM  OF  MANTISSA 

C  *  r0NfAT(0UT»fl#?5+IF8F»18»?3-TFPF) 

PFTllPN 

ENO 

SURROilTINF  T88I ( A#P#C ) 

CC-COt^f  ATfO»A»l#l6#l  ) 
CC*C0NCATfCC#A#9#Al#T) 

C*CONt ATI CC#B» 16*16*3?) 

RETURN 

FNP 
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APPENDIX  V 

COMPUTER  LISTING  OF  THE  SAP62  CHANGEZ  CODE 


sciro  LIST 

FILF  l«5AP/!Tnp»IINIT*OISK»sr*IAL»SLOCKlNG*l»PECnPO«5#0 
FILF  2bSAP/PT0T»UNITb0ISK.SFRI AL.BLOCKInGbI ,pECDR0b540 
COMMON  7(704) 

COMMON  /Fl  AGS/  IIT.IOT.CY.DMY(1«?0) 
r 
f 

I  it*  i 

IOTb  9 
CY*  S 

r 

r 

CALL  PFDT4D(7.704) 

?0  REID  1 00. T  »  F 

100  FORMAT! IS.ri5.0l 

1 F ( X • I F.0>  GO  TO  *0 
Z(  I  )b  F 
GO  TO  20 
C 

AO  CALL  WBTI 40(7.704) 

STOP 

END 

SUBROhTINF  RFDI40(7»N) 

RFIL  7(  N) 

COMMON  /FI  AGS/  I  IT. IOT.C Y. 77Z( 1 620* 

C 

REWTMr  ITT 

101  RE  AO  f TTT.FN0»2OO)  (77*( I)# !■! #540) 

RFAO(ITT)  (727(1*540). T»l»5401 
REAO(TTT)  (7Z7( 1*1000 i»I«l .540) 

00  10  T*1.256 

IF(I.GF.5?.AN0.I.LF.63l  GO  Tn  50 
CALL  T4pFf77Z(2*I-l)»777(2*I)»Z(T)) 

GO  TO. 10 

50  CALL  T4BI(77Z(2*T-1)#777(2*I)#7(I>) 

10  CONTINUE 

IF  ( 7( 2 )  -CY)  101.102.103 

102  00  11  TbI .256 

11  CALL  tABF(77Z(2*T*539). 777(2*1*540). 7(1*256)) 

00  12  TbI. 192 

12  CALL  t4RF(Z7Z(2*T*1079)»Z77(2*I*10R0)»7(T*^12)) 
return 

c 

103  PRINT. 104,7(2)»CY 

104  fORM»T(1»H1HaVE  FOUND  CYCLE  F5.0.19N  LOOKING  FOP  CYPlF  F5.0) 
STOP 

200  PRINT  201 

201  FORMAT! 24H1FN0  OF  FILF  ON  14  OISK.) 

STOP 

END 

SUBROUTINE  I4RF( A.R.C) 

INTEGFR  0100 

DATA  TR I As# 01 00/0*0000. 01 00/ 

EOUIVAI  ENCE(XE.IF) 

C 

C 

XE»  CONCAT(0»A»33. 17.15) 

IEb  ip!AS*39-IE 
C  IEP  15  PONFR  OF  0 

10  IE0b  TF/3 


c 


20 

C 

C 

C 

C 


99 

5 


50 

10 

100 


iepf  is  st  shitt  or  mantissa  at  ies 

IE8E*  IE-3*IE8 

iF(iepf,gf.o)  go  to  20 
IE6E-  IE4f*3 


IE8*  TEft-1 
IF8*  -TE8 

IF(  lES.LT’.O)  IE8*  OI0P-IE8 
PIrK  UP  SIGN 
(IUT»  rnNC»T(0» A»1 ,16, 1 ) 

PICK  UP  EXPONFNT 
OUT*  rONC*T(OUT» IEG»2,41  #7) 

PlfK  IIP  TOP  OF  MANTISSA 
OUT*  rONr»T(OUT#A#9*IF8E#32#l6) 

PICK  UP  BOTTOM  OF  MANTISSA 
c  «  rPNt AT(01IT»B»25*TF8F»  16»23-!FPE  ) 
RETURN 
END 

SUPROuTINf  f4RI(A»R»C) 
CC»C0nC»T(0#A#1#16#1  > 

CC*CONCATf CC» A»9»41 ,7) 

C*CONr AT! CC»B,I6»16»32) 


RETURN 


end 

SUBROuTINF  WRTI*0(7#N) 

SiSoI1!}.  AGS/  Ilf#I0TtCY*A(54Ol»DMYClOSOl 
REWINr  IOT 
PRINT  99 
FORMAT! 1  HI ) 

LL*  0 
NA  *  N 


K*  AN 

IF(K.BT.758)  K*25G 
DO  10  T*1.K 


lF(UfiF.52.AND.L.LF.831  GO  Tfl  50 
CALL  p5I4f(7(L)»A(?*I*1 l#Af2*IJ) 

GO  TO  10  .... 

CALL  PIAT(7(L),A(2*I-1)#A(2*T))  , 

PRINT  IPO,  Lf2(L)»7(Ll.A(2*I-n»AtP*TJ 

FORMAT(I1P»F30,6,E20,A,P020) 

WRITE! IOT)  CACI 3*1*1 ,540) 


NN*  N  N  •  W  ^ 

IF  (nn’.fo'.o)  60  TO  1000 
LL*  LI  *25S 
GO  TO  5 


C 

1000  REWIND  IPT 
RETURN 
FNP 

SURROIITINF  P5I4F! A,B»C) 

ECUlVil  FNCF  (XX, ID 
C 

XE*CONfATf 0, A, 41 ,2,7) 

IF(XE.FO.O)  GO  TO  10 
CALL  BT4F ! A,B,C) 

RETURN 

C 

10  XX*  A 
X"  II 
X*  X*P.l 
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x»  x-n.i 

CALL  Rl4FfX#B#C) 

RETURN 

EM» 

SUBROiiTINF  PJ4FCA#R»0 
FfilllVilFNcF  CXEMER) 

DATA  TRI Ar/OAOOOO/ 

IF(A.WF.O)  GO  TO  10 
B*  0 
c«  o 
RETURN 

c 

10  XS»C0NrATf0»A#47»2»l> 

Xf«C0wf‘*TfO*A#A2*3»6) 

IF(XS'.NE.O)  IE8--TE8 
IF8«-TF8 

B9«CCNfAT{0#A»A7»$#l> 

B10cCnNCATC0#A#A7#10#n 

IE8F*R 

IF(BIO.NF.O)  IFflF*l 
IF(B9>E.0)  IE8F«0 
IF  «3*TF6*TF8F 
IE*  !RTAS-IFt39 
EB«C0^CATf0#A#16#l#n 
BB»C0nCATCBB#!F#17,33#15) 

P«CONf ATfRP#A#32#9flE8F#16) 
C*C0NrATf0#A»16#2S*lE8F»23*IF8F) 
RE THRn 
ENO 

SURKOliTTNF  BI4I(A»8»C) 

BB«C0nC AT (  0#A#16#l#ll 
i*C»:iNr  ATCpR#  A#A1#9#7) 
C*CONrATfO#A#16»16#32) 

RETURN 

FM 

20  3  o’. 
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APPENDIX  VI 

COMPUTER  LISTING  OF  THE  SAPN/SKEWED  CODE 


sconr  list  7iP  dbuga  sijmry 

S  PROMOTION  SAVE 

BEGIN 

CREAL  PR0B.TIME»nFLT»DTC»PMAX*XPMAy#C02»Cl*TCl*FC02»DTZjM*OMTN#UMTN,WMl 
CRFAL  OTZJMR#PmAXRJ 

CTNT  N# JSTAR# JSTPl » jFTN*ZPROB»  ZN*ZTIME»70E| T*70TC*7JSTAR*ZJSTP1 I 

PINT  7JElN»7PMAX*ZJPMAX»ZXPMAX*Ze0?*ZCl*ZTCl*Zrcn2»7DT7JM#7OMlNl 
CTNT  JPMix*ZnMIN*ZNM*ZNSTOP,NSTOP*JSX»J*JJ,JMAX| 

PREAL  X3, RHDn»OV*CRNT»P2#e1 #P1 *GMONE*POU»RHOF »DT7J»DU»PFHDGE*Z| 
pr  REAL  wFCTpR  Ul41*Pt4l»Ot4 1 *RHOl 4)*Xl4I*7Mf  4 1# PPG T  AT  * 
CSr«1#TtAJ#RDRTmj 
LABEL  STAR?*  AGAIN*  STOPl 

«  FOLLOWING  CA(?D  FOR  SURROUTINF  AOJUSTSl FMFD. 

PINT  FROMpOW*  TOROMIROOL  EAn  TABLESFVI 
X 

«  VAPIABLFS  HAVF  BEEN  DECLARFO 
* 

ISET  UP  INPUT  InD  OUTPUT  FILES 

* 

FILF  I4DlSK1»"SAP"/"ITnB"<  51  ROWS  FULL)! 

FIIE  I4DISK?»»'SAP"/"BTni"C  51  RONS  FULL  1 1 

« 

PRFAI  SUBROUTINE  SORT  AS  RGAIPREAL  X  AS  RGA)I 
BEGIN 

SINCTNS  110  ASKRATCH  14 
SSTAPE2*  01 VPNIR/GSORT  SERIAL! 

FNDI 

SURROUTINF  PUTCPRICRFAL  X»PREAL  OUT  Z»rlNT  1)1 
BEGIN 

FOR  ALL  PfMrI  DO  Z*X| 

FNOl 

SURROUTINF  PUTCPICCINT  X*PRFAL  OUT  7*CINT  III 
BEGIN 
Cll  REAL  Al 
A*X  I 

FOR  ALL  PEW»T  DO  Z*A| 

ENDI 

SUBROUTINE  SLtWlPCPOINT  X*  CINT  JMAX1I 

BFGIN 

CINT  Jl 

LOOP  J»1*1.JMAX  00  XIJURTI  cZKj«MOOF*XtJl)l 
END! 

SUBROUTINE  UNSl FHCPCPOINT  X*  CINT  JHAy) I 

BEGIN 

CINT  Jl 

LOOP  J+l  *  1  •  JMAX  00  XI  JURTRf2Kj*M0DE»Xl  Jl)l 
ENOI 

SUBROUTINE  AOJIlSTSLEHFOf  PCPOINT  X#  CTNT  JMAXl! 

BFGIN 

« 

c  OUF  TO  THE  LACK  OF  DNN  VARIABLES  AND  OATA  STATEMFNTS  IN  GLVPNIR 

c  thf  fGiiowIng  Two  caros  must  appear  in  the  outermost  block 

S  IN  WHICH  ADJUSTSlEwrn  IS  USED, 

S 

I  PINT  fromrom*toromiboolean  tablesfti 
itarlesft^falsei 

!  MOOE  MUST  bf  true  on  entry  to  aojustslemed  to  pass  ali  args* 

I  IT  IS  TRuF  ON  FXIT. 

s 

CINT  JILABEL  MOVE! 


H7 


IF ( TABLESET )  THEN  60  TO  MOVFI 
t a 8i  fsft*tbufj 
TORON*r6A-PFN)  DIV  21 
FROMRON*RTl  ( 1 »MODF#TOROM) J 
T0R0W*RTR(1»M00E#T0R0N)I 

MODE«-BPOLFaN(5555555555555555(16)))FROMBOM«-TORONIMOOF«'TB|iF) 
TOROM-RTK  t»MODF»TORON)J 
MOVE  I  LOOP  J*0»S’»(JMAX-1>  OlV  32  00 
BEGIN 

LABEL  HOVFNI 

lF(JMAy*J*3?)  THEN  60  TO  HOVFNI 
M00F«-SHlFTB(?x(JO2-JMAX>»M00F>l 
MOVEMl  XT TOBON*J)*XtFROMROW*j)J 
ENOI 

MOOF«-TRUFI 
endixadjustsi  fhfo 

* 

*  START  or  NAIM  PROGRAM. 

* 

* 

stabti 

t  FOLLOWING  CARO  FOR  SUBROUTINE  AD JUSTSl FWEO . 

TABI  ESET*  FAUFI 
0PFN(IADISK1»0»1)I 
OPFNf IADISK9*0*1 )l 
jMAy+A  I 

SIHRfA0(IA0TSK?,7)I 

SIMNRlTECIAnT8Hl»Z)l 

ZPRnpKIJ 

7NB1J 

7TIMF*2l 

Z0FLT*3) 

ZOTCtAj 

?JSTAR*5I 

7JSTPi*6J 

?0HIN<7| 

7UMN*BJ 

7.if;n*9i 

7PMAXM0) 

7JPMAXM1 1 

ZXPHAX«-12l 

7C02H3I 

ZCIMAl 

ZTf  1051 

ZFr02*16) 

ZDTZjMUri 

7NW«-16) 

7NST0P*19l 

QMIN+GRAR0NE(Z*Z0MIN)J 
NN*6R*B0NFf7#ZHN)l 
0TZJH«-6R*RONE(Z»  ZDTZJM)) 

JSTAR»GRARONE(Z#ZJSTAR),l 

PRPB*GRARPNF(Z»ZPR0P)I 

N*GRAR0NFf7#ZN)J 

TIMF«-gRARONF(Z»ZTIMF)| 

OTC*GRABOnF (Z»ZOTC)| 

UHIN«gRARONF(Z*ZUMIN)I 
JFIN«-GRAR0NE(Z*ZJFIN)I 
C02«-GRAB0NF(Z»ZC02)I 
CUGRaR0NF(?»ZC1)I 
TCI  «-6rAB0NF  (7»7TCl)l 
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N5T0P4GRABPNF<Z#ZNRTDP1I 
fcnz«-r,PAPnwf{Z»zrcn?ii 
AGAIM  Of LT*-DTC  I 


DTODT7JM> 

OTZJP«-10flOOO,Ol 

PMAX40.0I  „  - 

SlEw(U»JMAy)lR|  fW(P#JMAXllRI  F W< 0# JMAX) I SLEHIRHO# JMAX) I  SI FWf X* jMAy > » 
SLFWfZM#JMAX 11*1  fW(fPG# JHAXllSLfWtCS# JMAXIISI  fW( T» JMAX II SI  FN(RDRT,JMAX)I 
JJ*JSTAR  HIV  621 


inop  no 

RfOIR 

wnor^TRuri 

pi«-pt.ii40r  jii 

MPOf  «-Rf  Vi  (PEN*f  J*6?ll<J$TARlANDfN0T  PDDlFAN<-mil 

Uf  J1*0t  Jl  +  nfLTxCPl-RTtt  1  #TRllF»Pl  >  >/<RHOI  Jl*fXt  Jl-RTlf-l.TRliE#Xt  JD) 
♦  RTLC1#TR|lE»RH0t  Jl)y(RTL(l#TRl'f  #Xt  JD-XT  J]))l 
IF(ABR<Uljl><UMTN)  THFN  UtjUO.OI 
XtJ]4Xt Jl4PTC*UtJll 
ENDI 

Mn0f4TPUri 

ADjllSTSlFWrP(U.  JMAXil  AOjUSTSlrWFO(X#JMAX)| 

J4JSTAR*6Ak JJI  If ( J<0 1  THEM  J4J46AI 

If (GR(»BONf fU(JJ)«J)>0.1  THFN  JSTAR4 JSTAR41  I 
IF f JSTAR2 JFINI  THEN  Gfl  TO  STOP! 

JSTPUJSUR41I 

LOOP  J40.1.JSTP1  DIV  62  00 
BFGIN 

H0DF4TRUF I 
P14XC jll 
X34P1*P1*P1I 

MnDE*REVi  1 2*  j#f  <PFN4(.|>«6?))SJSTPt  JANDfNOT  PODI  FAN(*1  1)11 
RH0N47Ht<1|)/(X3-RTL(-l,TRllF#X3n» 

DV*( RmOI JI*PHON) /( RHONyPHOt  Jl ) I 
0U4(ur  Jl»PTL(-l#TRllF#lit  JDII 


Ot J]4O.0» 


P0U40IU  , 

IF(DU>O.Ol  THEN  P0y40.0» 

If(0V<0.)  THEN  OtJURHONxIPDlIPCPP-ClxCRt  JlIxPOUl 
IFIOt.lISOMlN)  THEN  0lj|40.0l 
CRNT4O.0I 

If fOt jl>0’.0)  THFN  CRNT4TC1XCSI Jl-Ff 0?*POUI 


GH0NE*0 • A  I 

P24GMnNE*RH0NKFPGt JII 
FUtFPGt  Jl-OVPPt  JJI) 

IF(E1<0.1  THEN  E1*0.I 

Pl4GHnME*RH0Nyf 1 1  .  . 

EPGIJT4FPf,I  Jl-fP?4P[Jl40tJl40tJl)KPV/(?.4(P2-PlT/PU1'»l 

iF(rPftrji<o.o)  them  epgijuo.oi 

Pt  JJ4r.MPNF*BH()NKFPGC  Jll 
RHOf 4RH0N* 1 .0011 

-  IF(0V>0.P1  THFN  RHOf 4RH0N* «999l 

PFUOGF*GHnNExRHOF*F PGf Jl I 
CStJ3  4RORT((PtJl-^rUOGF)/(RHON-RHDn)l 
RHOt JI*RHONI 

DTZj4WWx(Xt  Jl»RTK"l#ThUf#Xt  Jl))/(PRNT4.014fSt  Jill 
DT7JHR4M!N(0TZJ)I  * 

PHAXR4HAXfPt  JDI 

If(0TZJHR<0T7JH)  THFN  OT7JH40TZ JHR I 
IF(PHAXR>PHAX)  THEN  PMAX4PHAXRI 

ifvPI 

M00E4TRUEI 


AP  JliSTSLFWrr  C  P*  >  |  *0  jySTRt  r wr D« 0#  JMAW  >  I  AO  ;*)ST8LFWrO(fS#.lHAX>l 

*D  JllSTSl  FWFOf  PuP#  JNAX)>»DJHSTSlF*lFI»tFP6**>  LVlt 
N*N*li 

TlM»>TTMF*OTCf 

pUTcPprpBPB#z»zPROe> J 

PUTCPP<DT7JM#Z#Z0TZJH)I 

PUTCPP(PMAX#Z#ZPM4X)f 

PUTCPPCTTMr#Z#ZT?»«ni 

PUTCPPCDTr#7#Z0TC)l 

PUTCPt(N#7#ZN)| 

PUTCPIf JST»P»Z#ZJST»im 
PUTCPPCDFl  T#Z»ZDFLT)J 

PUTCPPf XPM4X#Z#ZXPP4Xll 

PUTCPtC jPMX#Z#ZJPP»X>t 
PUTCP T  t J*TPl #  Z»Z JSTP1 ) I 

UBSI  E!(u!jM»niUNSLr5(pIjM21)lUNSLPH(®#  JHAX>pONSirn(PMO#J»»*X)l 

II M SI  rx(X»Jp»X)lUN8Lrw(ZM»Jp»X)lUNStrw(EPC»JMAX)lliNSLrwCfX#.IM*X)l 

UNSI  FH(T# jP»X>illNStEW(POPT.  jP*X)l 
SlMPPlTf  (I4MSK1»Z)I 

1F(N<nRTPP)  THFN  60  TP  4GMNI 

STOPI 

FNP. 


ft 
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APPENDIX  VII 

COMPUTER  LISTING  OF  SAP  IN  ASK 


!****•••*••****•  ASK  PROGRAM  SAP  . . 

BFGTN 

CHNS  6ai  SSEl FCT  64  PIT  WORD  MODE 

(PREAl  SUBROUTINE  SORT  AS  RftA(PRFAL  X  AS  R6A) 
(BFGTN 

*  SORT  ROUTINE  FROM  GLVPNIR 

DEFINE  COMMONbA.SKPATCH##! 

ASKRATCMfENTRVIi  I 


BLk 

1 4 1 

RASBStENTRYl I t 

NDS  64 I 

SORT! 

F 1 1  l 

1 

FXCHL (1) 

SD53I 

TXFFM( 1 ) 

#2l 

sl'ttmi 

b5i 

MAI  T 

i 

SToRrd ) 

SCSI 

(rfgtn 

FXrHLf  u 

SD53l 

JUMP 

GSORTI 

FII  t 

16* 

GSOTl  I 

DATA 

(?.??l5?)l»(3.16?l4)il 

DATA 

(*?«03146)l»( *5, 86116)1 

data 

(0.fll?5)l»(4, 75)11 

DATA 

(4000100000000000000001 

FII  l 

1  6  * 

GSOTPI 

STI (?) 

snr>0i 

RlN(l) 

SD40I 

1  I Tf 1 5  ■! 

#47#40l 

GS0T3 I 

LDi  (?) 

0(1)1 

LDS 

SC?! 

RTS 

•  Ol 

XI 

B  1  1 

TXl TAMf 1) 

»6S0T31 

LOl  (?) 

SD5BI 

EXfMLI?) 

SICRl 

GSORTI 

FII  l 

I 

Clr(?) 

1 

SI  TT(?) 

■GS0T2 1 

CLf(l) 

1 

SLTT(1 ) 

■GSOTl  1 

LDV  bARKRATCH*6I 

FXfMl (?) 

SICRl 

IDFF1 

SCOi 

IP 

0| 

SFTT  I 

•AND. El 

SETf (?) 

II 

ZF PTA(?) 

,i »i 

SITT(H 

b 1 ? |  f TRIFO  TO  FIND  SORTi 

HAI  T 

l 

JM7 

1 

SFTf (1 ) 

Jl 

ONFST ( 1 ) 

»GS0T4I 

lpx 

bOi 

SETF  -J 

.AND.Fl 

IB 

151 

SFTf (1 ) 

El 

SFTFj 

E.OR.Fl 

STA  ArKPATCH+5 i 
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6S0TAI 

IENDI 


SFTF 

I  •  AND  .F 1 

SETE1 

F.OR.Fl 

LDv 

■  1 1 

adfx  ASKRATCHaI?! 

SHAHR 
STt  1 

1 1 

ISKRATCHa$I 

Lorn 

sell 

IB 

11 

SHIR 

1 1 

SFTF 

I. AND. FI 

SETF1 

F.OR.Fl 

SIR 

1 1 

RAP 

?  i 

LOFF 1 

sri  i 

SETF 

-t.and.fi 

SETF1 

E.DR.El 

SIR 

7\ 

LDFF1 

SC1I 

ST» 

ASKRATCHa2| 

LOA 

ASKRATCHA5I 

LEX  ASkRATCNaIOI 

NORM 

1 

LITM)  “3.01 

SAN 

1 

STA 

ASHRATCH+31 

HLRN 

♦  IISKRATCHAJ 01 

SBRN 

* ASKRATCH+8 I 

HLRN 

ASKRATCH+3I 

ADRN 

♦askratch+ai 

STA 

ARKRATCHaAI 

HLRN 

SAI 

mlrn 

ASHRATCH+31 

ADRN 

SCI  l 

HLRN 

ASKRATCHa4| 

SHAMR  11 

NDRH 

1 

STA 

askratchaai 

HLRN 

SAI 

HLRN 

askRATCh+Sj 

AORN 

sri  i 

ML  RN 

AskRATCM+AI 

SHAHR  11 

NPRH 

1 

STA 

ARKRATCHaA  1 

HI  RN 

SAI 

HLRN 

A$KRATCHA3l 

ADRN 

sei  i 

HLRN 

AsXRATCHaAI 

SHAHR  ll 

NDRH 

1 

SAn 

1 

HLRN 

ARKRATCHaSI 

ADFX 

ASKRATCHAZi 

norh 

1 

LDFF1  *CO| 

EXPHLfll  »D5S| 
LOADMT  »c3i 
ILTTtn  *777777771*1 
EXrHien  *0531 
STiOt  *ICR» 
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FIl  t 


RFGTN  SAP  FXECUTION 


SPFXl  LITFO)  W 
STLFO)  '  * 
LDFE1  * 
LlTFOl  ■ 
STL FO)  * 
LlTFOl  ■ 
STL  F  0 i  * 
STl FO) 
LITFO) 


*17777777777777777777771 61 

*0331  *STORE  MODE  TO  ENABLE  ALL  PEIS  Tn 

*COI  *FNARLF  ALL  PFlS 

■  O.OJ 

*0351  *STORF  0.0  TO  *035 

■  II 

*0361  *STORF  1  TO  *036 

i)  *036 |  *STORF  WRITE  RECORD  COUNT 

i)  1#RA5BS+6A#rA5BSI 


STl  (0) 


*053l 


SImRFADFIA0ISK?»7)I 

CLCF3).  I 

SLITF •)  »AAJ 

CR0TRF3)  1*1 

CLCF?)  » 

SLITF?)  Zl 

CSMRF?) 

C0PF3)  *C?I 

CR0TRC3J  33 | 


LlT<?>  ■  ?!  > 

CORF  3)  SC ?  I  1 

CP0TRF3)  171  1 

CSPF3)  II  1 

DISKIO  I 

SlMWRITFFI4DISK1»Z)l 

ClrF3)  i  1 

SLITF3)  *441 

CR0TRF3)  l«l 

CLfF?)  I 

SLJTF?)  »7l 

CSMRF?)  *1 

CORF  31  *C?I 

CR0TRF3)  331 


•CLEAR  AC AR3 

(WORD  COUNT  TO  ACAR3 

•ROTATF  ACAR3  RIfiHT  1«  BITS 

•clear  acaR? 

%7  ADDRESS  TO  ACAR? 

•SHIFT  ADDRFSS  RIGHT  A  BITS 
•ADDRESS  10  ACAR3 
•ROTATE  A(AR3  RTfiMT  33  BITS 
•15  FOR  ADDRESS 
•  ?  FOR  CONFIGURATION  CONTROl 
•16  FOR  RECORD  NUMBER 
•FILE  NUMBER  TO  ACAR? 

•FILE  NUMBER  TO  ACAR3 
•ROTATF  ACAR3  RIfiHT  17  BITS 
•RFAD«1 


BITS  FOIVISABIF  BY  16) 


LITF?) 
CORF  3) 
CROTRF  3) 
C  R  R  F  3  ) 

OISKin 


•CLFAR  ACAR3 
•WORD  COUNT  TO  ACAR3 
•ROTATF  ACAR3  RIfiHT  1«  BITS 
•CLFAR  ACAR? 

*Z  ADDRESS  TO  ACAR? 

•shift  addRfss  right  a  bits 

•ADORESS  TO  ACAR3 

•  ROTATE  AC AR3  RIfiHT  33  BITS 

•15  FOR  ADOPFSS 

•?  FDR  CONFIGURATION  CONTROl 

•16  FOR  RECORD  NUMBER 

•FILE  NUMBER  TO  ACAR? 

•Fit F  NUMBER  TO  ACAR3 
•ROTATE  ACAP3  RIfiHT  17  BITS 
•WRITF-0 


FDIVISABIE  BY  16) 


STORE  PF  7  CONSTANTS  TO  ADR  F*D0-*031)I 


LOA  7 1 

CLCF1)  t 
CSRFll  01 
LlTFOi  1.31*01 
CLCC3) 

LDFF1  Ifll 


•LOAD  7  TO  BOA 

*pf  fnablf  moof  for  index  loop 

•LOOP  INDEX 

I  XCLEAR  ACAR3  FOR  INDEX 
•FNABLF  INOFX  PF 


LDCC2)  *AI  *2  TN  ENABLED  PE  TO  ATAR? 

STL(2*  *00(3)1  SSTORE  Z 

CApD(3)  *0361  SINCREmENT  AODRESS  INDEX 

CSHR(i)  II  f ENABLE  NEXT  PE 

TXlTM(O)  »-6l  SLOOP  ON  INDEX 

L0l(03  *n33l 

L0EE1  *fO I  SENARLE  All  PElS 

s 

S  DT7J*1 00000.01 

s 

LlT(l)  ■100000.01  s 
LD»  sell  s 

STi  DT7JI  TSTORE  OTZJ 

S 

LITd)  «63|  *63  TO  ACAR1 

CSUB(l)  S05l  S63-JSTART  TO  ATARI 

STL (1  )  *0371  SSAVE  63-JSTAR  IN  *0*7 

SAGAINl  OEi TsOTCI 
t 

AGAINlLOL (0)  SD4I  SOTf  TO  ACARO 

STL(0i  *n3l  SSTORF  BELT 

t 

S  0Tr*0T7JM| 

S 

LDL(O)  *0171  S0T7.IN  TO  ACARO 

STL(0>  *D«I  SSTORE  PTC 

s 

*  MODE ♦( PEN SJSTAfc)  AND  (NOT  ROOLEAN(-0))I 

s 

LOLd)  *0331  SSET  ALL  BITS  IN  ACAPl 

LDL(?)  *05|  SLnAO  JSTAR 

CADD(?>  *0361  S JSTAR  ♦  1  TO  AT ARP 

CSHR(l)  0(?)|  SSHTET  ACARl  RIGHT.  fJSTAR+1)  BITS 

coHPccn  i  scomplement  acari 

CRBd )  01  SCI  EAR  SIGN  PIT 

STLd )  *0341  SSAVE  MODE  IN  *034 

LPEE1  SCll  SSET  PE  NODE 

S 

t  U+jl+DEl  Tx((P*0)-RTi  d. TRUE. (P*0)))/fPH0*fX-RTL(-1 .TRUE. X)) 

S  SRTl  d.TRUE.RHO)x(RTLd.TRIIF.X)*X))l 

S 


LDl  (03 

S033I 

s 

L0EE1 

SCOl 

SENABLE  ALL  PElS 

LOA 

X) 

SLOAD  X  TO  RGA 

LOS 

SAI 

SSAVE  X  IN  PGS 

RTl 

SA.ll 

SRTL (»1 .TRUE.X)  TO  PGR 

SBRN 

sKi 

SX"RTL("l»TPUE»y)  TO  RGA 

MlPN 

RHO| 

SRHOxf X"RTLf *1 .TRUE »X3)  TO  PGA 

LDP 

SAi 

SSAVE  ABOVE  IN  PGB 

LOA 

SSI 

SX  EROM  PGS  TO  RGA 

LOs 

SR| 

SRGB  TO  RGS 

RTl 

SA.63I 

SRTLd.TRUE.X3  TP  RGR 

CM*  A 

1 

S-X  TO  RGA 

AOPN 

SB| 

SRTL(1»TRuE.X)-X  TO  RGA 

LDP 

SAI 

SSAVE  ABOVE  IN  RGB 

LOA 

RHO| 

SLOAD  RHO  TO  RGA 

RTl 

SA.63I 

SRTLd.TRllE.RHO)  TO  RGR 

LOA 

SRI 

SRTLd. TRUE. RHO)  TO  RGA 

MLPN 

SPI 

SRTL(l»TRUE»RHO)x(RTLd»TRUE.X)-X)  TO  RGA 

A  CRN 

SSI 

SDENOMTNATOR  TO  RGA 

LOs 

SAI 

SSAVE  DENOMINATOR  TP  RGS 

p 

ID* 

PI  (LOAD  P  TO  RGA 

ADRN 

Ol  »P*0  TP  PGA 

RTl  SA*6i|  SRTL(l»TRUF>(Pl05)  TP  RG* 

SBRN 

SRI  *(P*0)-RGA  TO  RGA 

L Ol  (05 

SD3A|  * 

iDrri 

SC Ol  SRESTORE  EXECUTION  MODE 

LDR 

■  Ol 

DVRN 

SSI  *(.....)  TO  RGA 

LDS 

$A|  SSAVE  . . ) 

L Ol  (15 

$D3 1  * 

LD* 

SCll  SDEIT  TO  ACAR1 

MLR*1 

SSI  *PELTx( . 5 

ADRN 

Ul  SUPDATF  II 

ST* 

Ul  sstore  II 

¥ 

* 

IFCARS(llXUMIN)  THFN  U*O.I 

¥ 

*0 

IN  RGA  FR(1M  RRFVIOUS  STATEMENT 

¥ 

RAB  >0 I 

SAPS(II)  TO  RGA 

LDi  (?) 

SDfll  tUMlN  TO  ACAR2 

I  *L  *C?> 

SSFT  PE  MODF  RFGlSTrP  f  IF  ARSCUXUMJN 

5ETC (2)  fl 

*PF  I  MODE  RTT  PATTFrN  TP  AC AR2 

CRB(2)  ol 

fCl far  acaR?  sign  bit 

LOI  (35  . 

*D37| 

CSHP(’) 

0(3)1 

CSHl (?5 

00)1 

LDFF 1 

SC'.*I  SENABLE  PEIS  FOR  WHICH  APS(tl)<IIMIN 

LPI  (35 

SD3'j  1  SlOAD  ACAR3  WITH  0.0 

LD| 

SC  3 1  * 

ST* 

Ul  f STORE  |i»0.0 

LDI  (35 

SD34I  f 

LDFF1 

Sc 3 1  f STORE  CURRFNT  FvECUTION  MOnF 

¥ 

* 

XtX+DTC*U 

¥ 

LDL(l)  spai 

SDTC  TO  ACAR) 

LD*  ill 

fU  TO  RGA 

MLRN  *Cll 

SDTCxU  TO  PGA 

ADRN  XI 

*X+nTC*ll  TO  RGA 

STA  XI 

fSTOPE  RG*  TO  X 

¥ 

I 

IF(6RApONFfU* JSTAR)>0,5  THFN  JST*R«JSTAr+1  J 

¥ 

CLC(O)  i 

SCI  FAR  ACARO 

SL  IT(Oi  Bill 

SPF  ADDRESS  OF  U  TO  ACARO 

CADO(O)  >P5I 

SAPORFSS  OF  II  I*'  JSTAR  PF 

LOAD( 05  SC1I 

SLOAD  U( JSTAR)  TO  ACAR! 

STL ( 1 )  SC?I 

SSAVF  U( JSTAR)  IN  ACAR? 

L1T(35  b63| 

SSHIFT  COUNT  TO  ACAR3 

C  SHR ( 1 5 

0(3)1 

ZFRT ( 1 1  .11 

SSKIP  IF  U( JSTAR)  i  0 

SKIP  ,61 

SSKfP  IF  U(JStAR)  <  o 

ZFRT(2i  ,51 

SSKTP  IF  U(JST*R)  »  o 

LDL(O)  SP36I 

SU(.ISTAR)  >  0.0 

CADD(0i  $P5| 

STL(O)  *n5l 

sstore  jstap 

CStiP(  35 

SD5I  * 

STi  (3) 

SD3? 1  f  UPDATA  63-JSTAR 

¥ 

t 

IF  ( JSTAR> ,|FIN5 

thFn  go  to  stopi 
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B 


* 

* 

* 


B 

B 

* 


« 

« 

« 


« 

* 

* 


B 

B 

B 


LDL(I)  . 

SD9J 

SJFTN  TO  ACAR1 

csuecn 

SD5I 

BJFIN-JSTAR  TO  ACAR1 

ZERXTti) 

» 31 

BSKIP  TO  HALT 

STL(l) 

*C?I 

BSAVF  ACAR1  IN  ACAR? 

CSHRfH 

?3I 

BSHTFT  ?*  PIT  NPRD  SIGN  PIT 

ZERXTM 1 

•  11 

BSKIP  IF  JSTAR  <  JFTN 

MALT 

1 

BSTDP 

JSTP1  «■  JSTAPM 

LDLtlJ 

SD36I 

BSFT  ACAR1  TP  1 

CADDC1  i 

S05l 

BJSTARM  TP  ACAPl 

sucn 

BP  SI 

bstprf  JSTAR+1  TD  JSTP1 

MDDF  ♦  SHTFTRC 1 

•  MPDF  DR  BOOLEAN  (»0l)J 

LDL(O) 

|D3*I 

SLOAD  CURRFNT  EXECUTION  MODE  TO  ACAPO 

CSB(O) 

PI 

BSFT  AC ARO  SIGN  BIT 

CSHRCOi 

11 

BSHTFT  ACARO  1  PIT  RIGHT 

LDEF1 

BCOJ 

bstpre  MODE  TO  PE  S 

STL (0  ) 

*P3*I 

BSTPPE  MODE 

X3  ♦  X*XxX 

LOi  (0) 

SD33I  * 

LDffi 

JCbl  BFNABLE  ALL  PEIS 

LDA 

XI 

BX.TP  RGA 

LDS 

SAI 

BSAVF  X  IN  R6S 

MLBW 

SSI  BXKX  TP  PE  RGAIS 

mlrn 

SSI 

SXSfc  XxXmX  TP  PGA 

BMUN  ♦  7M/I x3»RTLf-1 # TRUE# X3)) 

RTL 

*A»1| 

BRPIlTF  X3  LEFT-1 

SBRN 

SRI 

BXS-RTLI-lf TRUE#X3) 1  TP  RGA 

IDS 

SAI 

BX3-RTL(-1#TPUE#X3))  TO  RGS 

LDi  (0) 

SD3AI  B  ' 

LDrri 

SCOl  BRFSTORE  t'XFCUTlPN  MODE 

LDA 

7MI 

BZM  TP  PGA 

LDR 

■01  BCLEAR  PE  RfiBtS 

DVRN 

SSI 

BZM/(X3-RTL(-!»TRUE#X3)1  TP  RGA 

STA 

RMDNI 

B(POA)  to  phpn 

DV  ♦  (PMO 

■RHDN)/(RHONxRHD ) 

LDS 

SAI 

BRmON  TP  RGS 

LDA 

RMDl 

brhp  TO  RGA 

LDR 

SAI 

bsave  rho  to  pgr 

MLRN 

SSI 

BRHP*RHDN  TO  RGA 

LDP 

SAI  BSAVE  RHOxRMON  TP  RGP 

LDA 

SR i  BRHD  TP  RftA 

LDP 

SPt  BSAVE  RHflxRHDN  TM  RGP 

SBRN 

SSI 

BRHP-RHDN  TD  PGA 

LDP 

■01  BCLFAR  PF  RGBtS 

DVRN 

SRI 

B(RHD-RHDN)/(RHPNxPmP)  TO  RGA 

STA 

•Si  vi 

BSTORF(RHD-RHON>/RHPkRHP)  to  dlv 

DU  ♦  (li»RTl (•1#TRUF#U)>I 


B 

B 

B 


LDl  m 


SD33I 
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LOrri 


SCO)  (FNABLF  ALL  PEIS 


LOA 

III 

(LOAD  RGA  WITH  U 

RTL 

SA> 

11  «RTI (-1#TRUE»U)  TO  RftR 

SBRN  SRI 

( 

LOI  (0) 

$034 1 

( 

LOrFi 

SCO) 

(RESTORE  FXFCUTlON  MODE 

STA 

nui 

(Oll«<U*RTL  (•!  #  TRUE# II )  )  TO  OU 

LDft 

SAI 

(SAVE  OU  IN  PE  RftS 1  ft 

0  ♦ 

0.01 

LIT(l)  »0.0J  fO'. 

0  TO  ACARI 

LOA 

sen 

(ACARI  TO  RGA 

STA 

01 

((PftAxO.O)  to  0 

IFf  OUPfl’.O 

THEN  P0U*0.0  ELSE  PDU»Dll 

LOA 

oui 

(OU  TO  RftA 

LOI  fl) 

SD35I 

(0.0  TO  ACARI 

I  AC, 

sell 

XSET  PF  I  MODE  RlT  IF  DU»0.ft 

SFtC(?) 

II 

(PE  i  Mode  pit  pattern  to  acar? 

CRPf2$ 

Ol 

(CLEAR  PF  0  MODE  BIT 

LOI  O) 

snsri 

CSHPf^i 

00)1 

CSMl (?) 

00)1 

LOrri . 

SC?I 

(FNABLF  PEIS  FOR  WHICH  DU»0'.0 

LOI  O) 

$0351 

(LOAD  ACAR3  WITH  0.0 

LDa  . 

SCSI 

(0.0  TP  FNAPLFO  PElft 

Loi  O) 

S034I 

( 

LOrri 

SC  3 1 

(STORE  CIJRRFNT  EXECUTION  MOOF 

STA 

POlll 

(STORE  POll 

I 

«  IF(PV(0‘.  )  THEN  0*RMPNX(PDUvCn?»ClvCS|vPPllJ 

( 

LOi  (I) 

LOA 
I  Aft 

sftcm  1 

CRPfl) 

LOI  (3) 

CSmRI 1) 

CShL ( 1 > 

compcm  ) 

CRP( 1  ) 

CSwRM  ) 

CSMl  (11 
LOrri 
LOA 

Loi  rn 

MLRN 
LOS 
LDA 
LOR 

loi  m 

MLRN 
SBRN 
MLRN 
MLRN 
ST  A 

LOI  (0) 

LOrri 


$035l  (0.0  TO  ACARI 

OLV|  (LOAD  RGA  WITH  01  V 

SCI  I  (SET  RGO  I  BIT  TF  0LV>0.0 

II  (PE  RGO  ft  TO  ACARi 

Ol  (CLEAR  ACAR1  SIGN  PIT 

SOSTl 
00)1 
00)1 

)  (CDMPLFMFNT  ACARi 

Ol  (CLEAR  ACARI  SIGN  PIT 

00)1 
00)1 

SCli  IENABLF  PEIS  FOR  WHICH  OV<0. 

CSI  (LOAD  RGA  WITH  CS 

SOI  4 1  SLOAD  Cl  TO  ACARI 

SCI)  XClMCS  TO  RftA 

SAI  (SAVE  CivCS  IN  RftS 

PDUi  (LOAD  RGA  WITH  POU 

SAI  (SAVE  POU  IN  RQR 

$01 3J  (LOAD  CO?  TO  ACARI 

SC1I  IPDUvCO?  TO  RGA 

SSI  *(POU»CO?-ClvCS)  TO  RGA 

SR |  ((PDUvCD?-C1 vCS)xPOU  TO  RGA 

rhdn;  (RHONv(Pou»eo?-ei vcs)vpou  to  rga 

01  XSTORF  0 

S03A)  ( 

SCOl  (RESTORE  CURRENT  EXECUTION  MOOr 
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IF  f OSOMTN) 

THFN  0*0 

.1 

LDl 

01 

(LOAD  R6A  WITH  o 

LDl  M) 

sort 

tLDAO  OMIN  TO  ACAR1 

IAr. 

SCll 

tSET  ROD  I  PIT  IF  0>6MIN 

SETCM  ) 

ii 

XPF  R6P  $  TO  ACARl 

CRR(l) 

Ol 

tCLFAR  ACARl  SIAN  PIT 

LOi  131 

S03?i 

CSHRM  ) 

0(3)1 

CSMl <1  ) 

0(3)1 

COMPCO  ) 

1 

xcomrlfmfnt  ACARl 

CRPfl ) 

0| 

XCLEAR  ACARl  SIAN  RTT 

csi«R(i) 

0(3)1 

cshi  m 

0(3)1 

LOrFl 

SCll 

tenable  PFIS  FOR  WHICH  QSOHtN 

loi  m 

sn3Si 

f LOAD  0.0  TO  ACARl 

LDl 

SCll 

STI 

Ol 

XSTORE  0 

Loi  m 

S034I 

X 

LDFP1 

seoi 

XRFSTORF  CURRENT  ExFCUTION  MODF 

crnt*o'oi 

IFf 0>0.0) 

THFN  CRNT«-TC1*CS-FC02*P0UI 

0  IN  RGA  FROM  PRFVIOUS  STATFMENT 


L0K1) 

SD35I 

X 

LD» 

sen 

X 

ST* 

CRNTl 

XSTORE  CRNT 

LOA 

Ol 

XLOAP  0 

IAr, 

SCI  1 

XSET  RGO  I  PIT  IF  0>0.0 

sFTrri ) 

II 

XPE  RAO  S  TO  ACARl 

CRRfl ) 

01 

XCLFAR  ACARl  SIAN  RTT 

LDl  (3) 

S03TI 

CSWRC1 ) 

0(3)1 

CSML (1 ) 

0(3)1 

LDFF1 

SCI  1 

XENARLF  PEIS  FOR  NHTCH  0>ft.ft 

LDA 

POUl 

XlOAD  RGA  NITH  PftU 

LDl  (1) 

S016I 

XFC02  TO  ACARl 

MLpN 

SCll 

XFC02KPDII  TO  RGA 

LDS 

SAl 

X SAVE  FC02XPDU  IN  RAS 

LD* 

csi 

XLOAD  RG*  NITH  CS 

LO!  (1  ) 

S015I 

XLOAD  TCI  TO  ACARl 

MLRN 

sci  i 

XTCIxCS  TO  RGA 

SPRN 

SSI 

XTClxCS-FCOPxPDll  TO  RGA 

ST* 

CRNTl 

XSTORE  CRNT 

LDl (0) 

Sf>34| 

LDFFl 

SCOl 

XRFSTORF  CURRENT  EXFCUTTPN  MODF 

6H0NE  «-  0.41 

LIT(l)  bO.«I 

XSFT  ACARl  TO  0.4 

LOA 

SCll 

XLOtD  PE  RG*  S  WITH  ft. 4 

STA 

AMPNFl  1 

(STORE  0.4  TO  GMONE 

P2  • 

GMONFkRHONkEPGI 

HLRN 

l  RHONI 

XGMONkRHON  TO  pf  rga  s 

LDS 

fpgi 

XL0A0  PE  RGS  S  NITH  FPft 

MLRN 

i  SSI 

XGMONEkRHONxFPG  to  pf  rga  s 
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ST*  p?|  ISTORF  GHPNFkRHPNkEPG  TP  P2 

FI  ♦  (FPG-PV*P)I 

LO*  DLV I  *DLV  TP  PE  PGAIS 

CHS*  I  It-PV  TO  PF  RG*  S 

HLRN  PI  *»PV*P  TO  PE  R6*  S 

*0RN  SSI  *(FPG-0VkP)  TO  PE  RG*  S 

IF f F1«0.  )  THFN  rUOJ 


LOj  (1) 

*0351 

*L0*D  0.0  TP  ACARl 

1*1 

*Cll 

*SET  RGO  I  PIT  IF  Fl<0. 

SFTP(1 > 

II 

*PE  RGO  S  TP  ACARl 

cRRf  n 

LOi  (3) 

01 

SP3*I 

9CLEAR  ACARl  SIGN  RTT 

CSmRO  I 

0(3)1 

CSMl  (1  > 

0(3)1 

kfnapif  pf«s  FOP  WHICH  FI <0.0 

LDrrl 

SCI  i 

loi  (3) 

*035| 

HOAD  ACAP3  WITH  0.0 

LOt 

SCSI  *0.0  TP  EN*PLFD  PEIS 

LDi  (3) 

SD34I 

* 

LOPFl 

Sc3| 

SRESTORE  CURRENT  ExFCUTIPN  mPOF 
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APPENDIX  VIII 

B5500/6500,  ILLIAC  IV,  AND  CDC  6600  WORD  FORMATS 

The  B5500  and  B6500  have  the  same  48-bit  word  formats  because  the  machines 
were  intended  to  be  compatible.  The  difference  in  the  description  of  these  two 
word  formats  is  that  the  high  order  bit  in  the  B5500  is  labeled  bit  zero  and 
the  low  order  bit  is  thus  bit  47;  in  the  B6500  the  low  order  bit  is  numbered 
zero  and  the  high  order  bit  is  numbered  47.  Bits  that  can  be  manipulated  by 
the  system  have  been  added  outside  the  left-hand  end  of  the  word,  so  these 
added  bits  become  48  through  50.  This  distinction  is  relevant  only  when  the 
user  employs  such  functions  as  cocatinate  (similar  to  INBY  on  the  CDC  6600) 
where  one  references  bits  within  the  word. 

1.  B6500  WORD  FORMAT 
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Exponent 

MANTISSA 

Floating  point 
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mj 
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m 

INTEGER 

Integer 

The  B6500  floating-point  number,  or  integer,  contains  a  39-bit  mantissa  in 
bits  0  to  39;  a  sign  bit  for  this  mantissa  is  contained  in  bit  46,  word  format 
use  bit  47.  For  the  floating-point  word  format  an  octal  (SIC)  exponent  is 
contained  in  bits  39  to  44  and  an  exponent  sign  bit  is  contained  in  bit  45. 

Note  that  this  is  not  the  normal  biased  exponent  that  is  used  in  many  other 
machines . 

In  the  integer  format  the  fields  used  by  the  exponent  and  its  sign  in  the 
floating-point  format  are  zero.  In  the  integer  format  the  largest  integer  that 
can  be  held  in  this  format  is  239-l.  Beyond  this  the  hardware  will  convert  the 
result  of  an  arithmetic  operation  to  floating-point  format.  Because  the  binary 
point  is  assumed  to  be  to  the  right  of  bit  zero,  this  conversion  from  integer 
to  real  formats  is  a  continuous  one— an  integer  appears  to  be  just  an  "unnormal¬ 
ized"  floating-point  number  with  zero  exponent.  The  value  of  a  floating  point 
A  is  calculated  from  the  mantissa  M  and  exponent  E  as 
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2.  ILLIAC  IV  WORD  FORMAT 
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The  ILLIAC  IV  has  a  64-bit  word  that  can  contain  either  48-bit  integers  or 
a  floating-point  word  with  a  48-bit  mantissa.  The  binary  point  for  the  floating¬ 
point  numbers  is  assumed  to  be  at  the  high  order  end  of  the  mantissa,  that  is, 
between  bits  15  and  16  in  the  above  figure.  A  sign  bit  for  either  the  mantissa 
or  for  integers  is  contained  in  the  high  order  bit,  bit  0.  For  integers,  bits 
1  to  15  are  zero.  For  floating-point  numbers,  this  field  contains  a  15-bit 
exponent.  This  exponent  is  biased  at  400008  and  this  alleviates  the  need  for 
an  exponent  sign  bit.  The  value  of  a  floating-point  number  A  is  calculated 
from  the  mantissa  M,  the  exponent  E,  and  the  sign  bit  S  as 

A  -  (-1)S  x  2<E-40000e> 

It  is  noted  that  a  32-bit  half-precision  mode  also  exists  in  the  ILLIAC. 
Because  GLYPNIR  does  not  support  this  mode  at  this  time,  the  word  format  for 
the  32-bit  half-precision  mode  will  not  be  described. 

3.  THE  CDC  6600  WORD  FORMAT 

The  format  of  the  CDC  6600  integer  and  floating-point  words  are  included 
here  because  we  feel  that  any  calculations  done  on  the  ILLIAC  will  be  compared 
at  some  time  with  results  arrived  at  on  a  CDC  6000  or  7000  series  machine. 
Fortunately,  the  mantissa  in  both  of  these  word  formats  is  the  same — 48  bits. 

The  difference  in  word  length  is  all  absorbed  in  the  exponent  field.  (The 
extra  bits  are  needed  on  the  ILLIAC  so  it  can  have  one  bit  per  PE  for  logical 
instructions. ) 
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Exponent  MANTISSA 

_ I _ — - 


LJ _ 1 - — 

Floating  Point 

_ l 

63  48  47 

0 

VXZZA  -  —1 

Integer 

The  CDC  6000  series  machines  are  ones-compliment  machines,  thus  there  is  no 
sign  bit  and  the  negative  of  a  number  is  its  ones-compliment.  Integers  must  be 
less  than  248  in  size  so  that  they  can  be  processed  with  the  same  hardware  that 
processes  the  48-bit  mantissa  of  a  floating-point  number.  Although  the  bits  48 
to  63  are  not  used  (i.e.,  are  zero)  for  positive  integers,  there  is  sign  exten¬ 
sion  through  this  region  for  negative  numbers,  and  thus  they  are  all  on  in  this 
case.  In  the  floating-point  format  the  binary  point  is  assumed  to  be  a  the  far 
right  of  the  mantissa,  and  the  exponent  in  bits  48  to  62  is  biased  at  20008. 

Because  the  mantissa  of  the  CDC  6600  word  and  that  of  the  ILLIAC  are  of 
the  same  length,  the  movement  of  algorithms  from  one  machine  to  the  other 
should  be  simplified,  as  should  the  movement  of  data  and  calculations.  The 
movement  of  data  for  later  calculation  between  machines  of  different  word 
(mantissa)  length  presents  a  number  of  problems  with  the  stability  of  the 
calculations  in  question. 
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SCODF  LIST  PPUfiA  SUMRY  SAVE  7IP 
PEGU 

*  THIS  T$  #2  ZONE  (OR  LESS)  WIPE  SHELL  IN  GLvPNIR. 

CINT  IMAY#JMAV».|PLK»J»J1,J2| 

CREAI  PRnB»rVCl.F»T,DT#DT|.A$T.WCA»NCP#NCC#WCD,wCF»Pl»FTH,STARl  FCT»| 
CU  RF  AL  RFLFRRPFEORE#  C  YCl.  ASTPR I  NT  »MA  XRFl  FPR*  G0*TP70NE  I 
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Y 
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FILF  LINE  (1  ROW)) 

Y 

* 
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t 

1 

SUPROUTINF  PUTCPR(CRFA[  y, PRFAL  out  Z#CImT  III  PFGIN 
IE  PFN.T  THFN  7*X  I 

FNPI 

* 
t 
« 

« 

SUPROUTINF  PUTCPKCINT  X#PRFAl  OUT  7»CINT  1)1  RFGIN 
CREAL  A I  <*XI  IE  PFN«I  THFN  Z*AI 

FNPI 

* 

* 

* 

X 

PREAI  SURROiiTINF  LN  AS  RGA  f PRFAL  X  AS  ROA)|  BFGTN 
S  INrlNS«10O#  ASWRATCH-1P 
s$tapf?»glypnir/gln  serTali 

FNPI 

Y 
* 

Y 

y 

PREAI  SUBROliT  I NF  ExP  AS  RG*  fPREAL  X  AS  RGA)I  BFGlN 
S  TNrlNS«100»  ASWRATCH«17 
$STAPE2«GLYPNIR/GEXP  SFRTALi 
FNPI 

y 

y 

y 

y 

PREAI  SUBROl'iTlNF  SORT  AS  RGA  (PREAL  X  AS  RG*)I  PEGIN 
S  I NC I NS»1 Ip#  ArkRATCH»14 
*STAPE2»GLYPNIR/r.S0RT  SERIAL  I 
FNPI 

y 
y 
y 
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« 

PREAl  SUBPOuTINF  ALOGlOtPREAl  X)l  BEGIN 
*LOG!0«-0.«S4?9<  *6l9o  32*51  «>«LN(X)I 

FNOJ 

X 

X 

X 

X 

PREAl  SUBROiiTINF  POwEROFlO  fPRFAL  N)l  BrdN 
*  POwrpnFio»io.**N, 
P0NFR0F10*FXP(2,302Se5o9?99A0*6«w)J 

fnoj 

X 

X 

X 

x 

SURROUTINF  pnwFRfPRFA'.  X.PRFAL  OUT  XX,PInT  OUT  NU 
XX«-A10g10(ARS(X>>| 

N*(xx-n.5*RTGN(xxm  innwT  lft  glypwir  pound. 

IF(N<Oi  TMFN  N«-N-1J 
XX*XKPnwrBnMO(-N)l 

FNB  J 


iUPROUTINE  NFGATT VEM*SS(e INT  J)l  BFGIN 
CTNT  J.'ll 
MODF«-iMonr» 

loop  j.i«-o»i#jmax  r»o  begin 

irf AMyt jjjso. jthen  Bfgin 

JiIMwRITFfLTRFrPAGE]»"NFG  MAS5"»"j  « 
SlMWRITrfLINF."APXf *]»»AMX[*J)I 
FPPOR(O) I 

Ewm 

END  I 

HOI  X  NEGATlVFMASS, 


SUBROUTINE  wFGATlVEENEPfeYKfNT  J}|  BFGlN 
Cl  NT  J.'ll 
HODE*TMOnn 

LOOP  J,‘i*Ot1#JMAX  00  BEGIN 

iFfAIytJJISO. )THEn  BEGIN 

SINWRITFO.  TNFrPAGFI*"NFG  FNFR"»"J  ■** 
RIMWRITFO  INF»"AlXl*1»,AlXt*m 
FPPORCO)! 

fnbj 

END  I 

MODE*TPUEJ 

SINWRlTFMAniSKl#7)l 

SlMNRlTFfLINF#«MASS  ANn  FNERGY  OK  AT  J  >  ",J1J 
FNDI  X  NEGATIVFFNERGY. 


RFAl  SUBROUTINE  SUMROW ( PRE Al  X)l  BEGIN 
ClNT  J.THBTOTHFJlBfinLEAN  OLONOBEI 


BFGIn 


»  J)l 


.  J)l 
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•am 


U1.111  ig'.wi— — mmmmm  mhmhbmhmr 


oldmode *mopf i 
moof*not  Mppri 
x*o.j 

M00F«-TBliri 

TWOTOTprjM  .1 

LOOP  J*0#1,5  DO  BfcnlN 

X*X*RTR! TWOTOTHFJ*mOPE*X)I 
TWPT0TMFJ*SHIFTI!1 .TWOTPTHEJU 

ENDI 

5UMR0N*XI 
M0DF*0i PMOnFI 
NO J  X  SIJMRPW 


SUBROUTINE  rPRORCCINT  J)»  RF6IN 

simwritfil tnfipagei* "Error  ■"» jn 
WCA*0.I 

WCAM.p  +  l'.O/WCAl  t ABORT • 

FNOI  (ERROR. 

* 

« 

« 

t 

SUBROUTINE  POT  I  REGIN 

PRFAL  wTNPJMFNSIONI 

CREAL  rfr,tnfxt#woa,wdp.wov»  CINT  LP#JJ  rOOLFAN  Ml »  MR | 

WDA*1.0R-?OI 

LOOP  J*0»  1  *  JMAX  DO  RFP.TN 
MOOE*TRUFl 
J1*J  nlV  A4I 
J?*J"A«KJ1I 
Nnv*r,RABONE{ovr 

mppfmmodei 

CS«-AMXf  Jl/f  WnVKTAllW 

EsfJJlXES  RETURNS  IIPPATED  P  ANP  PS  FOR  THF  ROW  ,V. 
MINCIMENSIPN^PX) 

IrfWPY<PX)  THEN  MTNPIMFNSION*WpYl 

HOB*  MAX{CS/MlNnlMFNSlON)l  «  TTMF  FPR  SOUND  TO  CROSS  MIN  CELL 
IF  !WPA<WDB>  THFN . WDA+WDR) 

WPP*  MAX{ABS(Ur Jll/Pxm  TIME  FOR  MATFRTAl  to  cross  peli.  tn  X 
f  DIRECTION 

IE  !WPA<WDB)  THFN  WPA*WDBI 

HOB*  MAXIABSfVr J])/WPY)j  t  TIMr  FOR  MATFRIAt  TO  CROSS  CfLl  IN  7 

*  PIRFCTION 

IF  (WOA<WDB)  THFN  wOA*WDBt 

FNOI 

MOOE*TRUEI 

0T*0.5*STAPIFCTR/W0AI 

DT*1 

T*T*DTi ASTI 

POWER! T#BfR#lP)l 

TNEXT*fT*OY)*POWER(lFlOf-LP)l 

Ml*  (TnFXT>ST0T)  ANn  BpoL F AN!  1  FFF7F7777'77600C000000f Bill 
M2*SHIrTL  f 1 #M1  )  J 

Ml *! Ml  PR  M?)  AND  WOT  !Ml  AND  M2)l 

IF!M1  tND  f  PFR<STOT  )  )  THFN  PT*!  STOT-BFR)kPOWFROF10!LP)I 
DTLAST*PTI 
CYCLE*PYCLF*1I 
FNOI  SCOT. 
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r 

t 

i 

t 

t 

t 

i 

* 

subroutine  pFznwr  <cint  jii  bfgin 
FRROPdOOIl 

FNOI 

t 

* 

* 

* 

SUBROUTINE  fS(CTNT  K)l 
BEGIN 

*  ES  USES  RpP  TIN  CS  ARRAY! »  AND  AI*!*!. 

*  IT  RETURNS  UPDATED  PfK1  AND  CS. 

*  MODIFICATION  OF  AIR  TO  FXTFNO  TO  2.F13  FRGS/GM.  II  MAY  69 

PRFAL  F  *  L  NFRPVPOZ. FP»  FPN»  FN»  HS»  bfta  »  FE»  6M0NF I 
PREAL  RMO  AS  CSI 
E*1 .0#-10*ATTtK!l 
LNEPOVpPZ*! N(773.39*?0A?5hRH0)I 

WSO  . 0 /< 1.0*FXPt (-8.5-0. 1550431 «KlNFROVPPZ*F)HEXPr -0.05*1  NEPOVNOZ 
♦o.OP53178o798!))» 

FO  ♦  NS  neXP(-0.2242i52a664*f)I 
FON*( 1^ O-NS  JNEXP(*0»150R2956259*F5f 
BETA*O.I 

IF(E>1.)TPFN  BETA*  ( 0 . 0069A87*MS*0 • Ol 38974 )*LN( E ! I 

WS*1.0>( 1.0*FXp((-f*EXp( 0.01 57*LNFRPVP07*3. *06662489!) *FXp(-0.065 
*(  NFBPVROZ-l. 3*629436)))! 

FN*WS*FXPf-0.0392l5666275*E)l 
BETA*BfTA*M .o-ns)+o.oas*ws» 

NPAOO*.  0*3.474  356*LNFRpvR0Zl 
IF  WPA  <6’.0  THEN  MPA*6.nl 
NPA^(160.0-F)/WPA» 

FE*1  .0>m'.0*FXP(nPA))I 

GMONEtf 0,1 ftl ♦0,255*FP*0.280*FON+0, 1 37xFN*0,P5*FE )*EXP( BETA 
*1 NEPOVROZ)! 

P(  K !*GmPNF*BHOxAIX t K 1 1 
CS*SCRT((1‘.0*GMONF)kP[K!/RHO)» 

FNOI  *ES, 

* 

* 

X 

t 

SUBROUTINE  INPUT!  BEGIN 
MODE»TPUEI 
SINRFAOf I«0I5K0»Z)! 

P*dB  ♦ftP»BONF(Z»OI! 

CY)CLE*f,PABONF(Z»l  I! 

CYCLASTPR!nT*CYCLFI 
T  ♦fiPABPNF ( Z»  2) I 
OT  ♦gRABPNF(Z»31I 
DTLAST*GRABONE(Z.«)» 

T*T-OTi AST| 

CYCLE*CVCLF-1! 

IMAX  ♦f,RABPNE(Z»5l! 

UMAX  *fiPABpNE(Z»6)l 
ETH*GRAB0NF(Z»7)I 
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STABLFrTP*fiRABONF(Z»R)l 
IF<STAPIFCTP<0,)  THFN  ST ABl FC TP* 1  . 1 
G0*GRAPPNFf7»9)J 
G0*0,0i 

T0ZPNE+GPARPNE(Z»10)I 
JBIK*JmAX  PTV  6AI 
TOPREFl  FCTTVF  ♦  FALSE! 

RIGHTRrFLFOTIVE  ♦  FALSE! 

BOTTOHPFFLfCTIVE  *  FALSE! 

!F(GRApPNFrZ*ll)PO,)  THFP  toprfflfctive*trufi 

IF(GRAPPNFf7»l?)PO,)  THFN  RIGHTPFFLFCTI  VF*TP(»E  I 

IF(GRApONEf7»13)pO,)  THFN  BOTTOMPFFl ECTIvE*TRUFI 

RElFRRPFrORE*GRABONF(Z.lA)J 

HAXRFLFPP*GRAB0NFfZ.I5il 

DX*X“RTPr1>MPOE»XlJ 

RC*o.5yf  x*ptr(  i»Mnor.xn> 

MODE *BnPLF ANT “0)1 
MODF*SHIFTR(IMAX*l.MODnl 
DX*RTRf?*»pX)J 
RC*RTRf l*»X)*0,5xp*» 

M0DE*BPPLFAN(“0)l 
DX*  RT|  n».nx)i 
RO-RTl  M».pr  )J 
MODE  *TPHF I 

IMCOE*wPT  PPOLEAN  (-0)  AND  PENSIMaXI 
WPA*Xxvi 

TAU*  PIxrWPA-BTR<!»MpnF,NPAJ)J 
LOOP  J*0»1 * JBLK  00  PEG! N 

ovr jJtRTUi»MOPE  .vrjDi 

iFfJPjBLK)  THFN  BrGIN 

TF  BOOLF*N(!)  THEN  PVtJJ*  RTLfl  »»Yr  J*M  >1 

FNPI 

DVt  JWOYt  J]-Y[  J!l 

END! 

HP  I  XINPUT'. 


SURROUTINF  PAT* I  BEGIN 
CODr 

USE  STOT»P!l  PFGTN 

hfpetami  SET  > 

OPG  PTI 

DATA  ».l«159?6535R9793l 
OPP  STDTI 

DATA  tirMirMi9»2iA«2^l<7i3>2<5«7it» 

?,6*3,0»3.?*S,A»3«6»3,fl»4,0»4.?»4,5»4,6»5.  P»5,3»5,6»fi’.  0* 
6*5»7«0»7,5»R,0»8*5*9, 0*9,5*10. 01 
OPG  HFREIAHI 
END  COnFI 
FNPI  XOATA. 


SUBROUTINE  pH  1 1  BEGIN 
CINT  PASS!  CBEAL  61 
PREAl  uA»VR»PA»PP»Bb»URPI 

mode*tb«>ei 

DTH*0. 1 
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URR*O.OI 

IT  RIGHTRrFLFCTIVF  THEN 

ERROR! i 5 

EISF 

BFGIN 

I  SFT  P  AT  AXIS  AMR  RIGHT 
MnOE«-R0OLEANt-0)l 

LOOP  j*0# 1 # JMAX  OP  Pt J)*RTL(1##PU))I 

MnOF*$HlFTRtTMAX+1 .HOOF) I 

LOOP  ,1*0#  1  #  JHAX  On  Pt  J)*RTR(?##Pt  J))l 

FNO) 

LOOP  PASS*) » 1 #2  on  begin 
MnPF*TRHEl 

If  BOTTOMRFFLECTlVF  THEN  BFGIn 
VR*0#l 
rf»PtO)i 

FNO 

El  RE  BEGIN 

VR*VtO)l 

WCA*  GRA80NE(0Yt0)#0)J 
WCB*  GRABONEtOYtO)#l )l 
WCC*1 ./fWCA*RCB)l 
PP*(Pt  l)wj<CA+PtO)i«WCR)»«WCCI 

OTH*DTH*0.5*TAHKWCA*VmHPTH(PrOI*Pm)xNrrl 

FNPI 

IF  RTgHTREFLFCTIVF  then  FPRPRt?) 

El  RE  BEGIN 

t  SET  II  AT  AXIS  AND  RIGHT. 

Hn0E«-B0niFANf»0)l 

I  OOP  J*0# ! # JHAX  00  Ut J)**PTL(1 ,#Ht Jl)l 
M0PE«-SHTFTR(IMAX*1»M0DF)I 
BPA*RTR(l##Rf )/RC) 

I  OOP  J*0 # 1 # JHA X  00  Ut JHWCAxRTPt 1 ##Ut J))l 

FnP  I 

MnOE*|MOOF  OR  BOOirANt-O)! 

WPA*r.O/(OX*RTL(l.#PX))l 
WPP*PcxRTl.(  1  #  #PX)xWPAJ 
BpC*BTL(1»»BC)kDXxWPAI 
HPPF  *TRUE I 

WPA*GRABONF(DYtO]#0)l 
HOPE*THOOF  t 

Loop  ii*o»!»jhax-i  no  begin 

TFt AHXt J)<0. )  THEN  FRR0R(3)I 
|(RR*Ut J)xwPB+RTLt !##Ut J))xHPC> 

M0PE*TRUE) 

,ll«-(J*l)  PIV  AAI 
.!?♦(  J+D-Mx.'ll  I 
WCR*GRARONF(OYfJl ) # J? ) I 
M0DE*IM0DF> 

WCC*1 .0/(WCA*WCB)I 
MCD*WC AnNCC) 

WCE*NCBxHCCI 

VA*Vf  J)i«WCF*vr  J*1  IxNCOI 

Aixt  J)*AlXt  J1*Pr  J)*OTxtO.«iX(VB-VA)xTAH*(RTRM»»llRR)-liRR) 
xP|xwCA)/AHXt J) I 
IF (PASS«1 )  THFN  PFGIN 

PA*PtJ)*WCF+Pt J+l)xWrOI 

MOOE*TRUFl 

G*6#3567BBI 

G*G/(G+GBABONF(Yt J1 ) # J? >-0 .5xNCA ) > 
g*gGxg*g> 
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moofmmopf  or  booleanc-oii 

PR*fP[J3«»RTL(l##0*)*BTLC1,#Prjn»0*)»cWP»l 

moofaimodfi 

Ut  J  ]♦  lit  JUT*UKHC*K(RTR(l#,PR)-PR)»rpT/fO*KAMXrj])  I 
V(  JUVC  Jl*TAU»«WCA>«(PR-PA)>«OT/(WCA>«AMXf  Jl)-GxOTI 
PB*P|| 

FNOI 
V R»VAI 
WCA*NC8l 
FwPi  t  J  LOOP 
J+JWAxl 

iFf AMxrJlSO.)  THFN  FRNONC«)l 
URN«-tlr  Jl*WPB*RTLC1  »»Ut  J])XWPCI 
IF  TPPREFIECTIVF  THEN  FRN0BC5) 

FI SF  PFGIN 

VA«-Vt JII 

t  P/UPBI  X  ALRFADY  OONF  ABOVE 
0TH«-0TH-0'.5KT*U»«HCBKVt  JjKf>T»«CPr  Jl+PCJ-njKWCCI 

Ffein » . 

ATXC  jl«-MXt  JI+PC  JlxDTK(0.5KrvB-VA)xTAU+(RTRf  1#  *UNN)»UNR) 
xPT*WCA)/AMXr.lIJ 
IF  PASS«1  THFN  PEfirN 

MOOF*MOOF  OR  ROOIFANC-O)! 

PRKPt  JIKRTLM  »#0X)*RTLC1»#PI  J]  )XOY)xWPAI 
MOPEMMOOFI 

Ilf  JJ*UC  J]milxWCAx(NTRCl##PN>«pR>xnT/COXXAMXr  J])J 

E*PI 

X  niTNF  WITH  WHOLE  GRID 
END!  X  PISS  LOOP. 

F TH^ETM* SORROW C OTH ) i 
FNOI  I  PHI 

SUBROUTINE  pw?l  BEGIN 

PF  RF A|  SPFCTFICFNERGY»OPNORTOPU#OONORTOP\/#OONORTOPSPFClFlCFNFRr.y# 

VlTTOPTNTFnF*CC:,OONORM4RSOlVlnFOBVOV#N*SRTHP|iTOP,MlSSTHRUBOTTOM. 

OONORRI GHT||»OONORRl ghtv*oonorrightspfciftcfnfrgy* 

DONORMI5SOTVlOEDBYDX#UlTRIGHTINTFRFICF#MISSTHRURTeHT,lHXNFw» 

ONEOVFRIRXNFW,UMOMFNTURTHRUTOP#VNOMFNTUMTHPUTOP,|iMOMFNT|inTHRURIGHT, 

VMONFNTIIRTHRURlGHT#FNERfiYTHRURlGHT#FNFRGYTHRUTOP.FNFP6VTHRUBPTTOM# 

UMOMFNUiMTHRUBOTTOM,vMnMFNTUMTHRUBOTTOMl 

RF2ONFTPP«-F*LSFJRF7ONFRTGHT*F«LSF|RF20NFP0TTOM»FILSp» 

0TH«-0.i 

modf^trufi 

WCI^GRiPO'-FfP  i'C0]»0)l 

*  SFT  IXIS  BOUNDARY  CONOlTlONS 

MODF«-BOOLFaNC-0)I 

HASSTHRURtfiHT«-O.I 

X  SFT  BOTTOM  BOUNOARV  fONpiTl  VS 

MOOF*TRUFI 

if  bottomPfflective  THFN  bfgin 

MARSTHRUBOTTOM^O.i 
ENFRftVTHRUBOTTOMAO. I 

umomfntumthruBottpm«-o.  I 

VMOMFMTUMf HRUBOTTOM*0»  I 

ENO 

FLSF  BFGIN 

MASSTHRUBOTTOM*vm*AMXCOJ*OT/WCAJ 
!FfMASSTHRUBOtTOM»0.)  THFN  MASSTHRiibOTVOM^O. | 

I Ff •MASSTHRUB0TT0m>T020NE xAMX(Ol)  then  REZONFBOTTOMaTRUF) 
ENFRGYTHR'JBOTTOM«-MASSTHRuBOTTOMxCAIx:OU0.3xCUrOI»«IIIOUVfOl* 
vrom.i 

OTH»FnFRGYTHRUROTTOM| 
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l'MPMFMTtlMTHRll&PTTnM*MASSTHR»lpnTT0Mx|it01l 
VMPRF*TtlMTHRllBnTTnM«-MASSTHRUflOTTOMKV[0]  I 

FNDI 

IF  RIGhTRFfIFCTIVF  THEN  FRRPRC6) 

FLSF  BFRIN 

X  SFT  RIGHT  BOUNDARY  CONDITION* 

MpPE«-SHIFTRCTMAX*1 .MODEM 
IrtPP  j«-0»l#JMAX  OP  PffilN 
liT  J1*RTR()#»IT  JIM 
tfcu[ji<o.)  tmfn  utjjx-ur jii 

EwPI 

FNDI 

MODE«-TMCIPFl 

LOOP  J+O.  1  •  JMAX  DP  BEGIN 

X  CPMpUTF  OUANTTTTfS  AT  TPP  INTERFACE 
SPFCTFlCENERGY*Alxr  J]KU[  Jl+vr  JlKVr  JIM 

DOPORTOPIX-Ur.ll  I 
OPNPPTPPVxVf J]l 

OPNPRTPPSPECTFICEEiFRGY+SPFCTFlfENFPGYl 
DPNPRp A  ESDI VT  DFOBYPY  +  AMXf  J ] /WC  A I 
IF  J<JM AX  THFN  BFBIN 
WPPE«-TRIIFI 
J1«-CJ*1  )  PIV  AAI 
J?*(J*n-6AKJ1l 
WfR*GRABPNEfPYrjll»J?M 
MOOEMMOOEI 

VATT0PINTFPFACF*0.5xfVr  JUVt  J*1 1M 

vATTOPINTEREACE«-VATTPPTnTFREACF/(1  .O  +  mj+1  1-VT  JD/NCA 
*PTM 

TF  (VATTPPINTFRFACE<0, )  THEN  pEGIN 

DONPRHASsOIVlDFOBYPY«-AMXr.l  +  n/NCBl 
DONOPTOPIM-Ut  J+1M 

donortopv«-v[  j+i  m 

DONPRTOPRPECIFICFNERGY«-ATX[  J*1  J*O.Sxfii[.l«-l  ]X|ir  J*1  !♦ 

vr  j*i  i«vr  j+ini 

fnpi 

MASSTHRIlTOP«-VATTOPINTFPFAfFxDT»DONORMASSDTVlPFDBYOYI 

ENP 

El  SE  BEGIN  t  DP  TOP  ROW  OTTERENTLY. 

TF  TOPREELECTT VE  THEN  FRRPR( 7 I 
Ft  SE  PEGTN 

VATTPPlNTERFACFxVr Jll 

MASSTHRIITPP*VATTOPINTEpF  ACExPTxOONORHASSPI vipfopydy I 
IF{MASSTHRUT0P<0.  )  THEN  MASSTHRUTPPxO',  I 
lFfMASSTHRllTPP>T070NFXAMXf  JM  THEN  RE70nETPP«-TRUEI 
DTHxOmMASSTHRHTnpxOONOPTOPSPEClFICFNERGYl 

FNPI 

FNPI 

X  COMPUTE  QUANTITIES  AT  RIGHT  TNTERfACF 
OPNpRp I GHTU*Ht  J]l 
OPNOPRlGHTVxVr J]l 

DnNOBRTGHTSPFCIFlCFMFRGY«-$PFClFlCENFRGYI 
OpNPPmARSDIVIOFDBYPY«-AMX[J]/(DXxRC)I 
UATRlf,HTINTERFACE«-0.5x(Ur  Jl+RTI  ( 1  »»H[ J1 ) M 

UATRlf,HTlNTFRFACE«-llATRlGHTlNTEPFACE/(I.O*(RTl  (1  »#lir  JI)»nr  JM 

/OXxDTM 

IF  OIATRIGHTINTEReACE<0.)  THEN  begin 

PPN0RMhS5DTVT0FDBYDX«-RTL(1  ##donormassoi  vs  pfdbydxm 

pnNORRlGHTtJXRTl  C  1  #  #OONPRRTGHTllil 
PPNORRlGHTV*RTl(I # .DONOR RIGHT V  > I 
PPNORRIGHTSPFCIFICENERGY* 
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RTLn»#PPNORRIGHTSPFCIFICFNFPGY)J 

Fwm 

MARSTHRURIGHT«-llATRTGHTlNTFRFACF>iOT>ipONPRMA3SPIVTDFnBYDXiiX» 

MPPF«-TRI'FJ 

Irf  GPAPONf (MASRTHRIlRIGHT»IHAX>>TnZPNE*GRABnNF(ANXf  J]#IMAX)) 

TUFN  rF7PN£RIGNT«-TRUEJ 

MpPEMMPDEJ 

AMXNFW«-AMXt  J 1-MASRTHRUYOP+M  A  SSYHPURPTTPM-M  A  SSTNPUPT  r.HT* 

RTP( 1 #»MASSTHBURlGHT) J 
PKFnvrPAMXNFR*! .O/AMXNENJ 

*  POMrNTUM  FIUXFS  AT  TOP  ANP  RIGHT 
UpPMFnTIIMTHBUTOPTMASSTHRUTOPhOPNOBTPPU I 

umphfntumthruRight«.pasrthrupight»cdpworpightiii 
VMPHFKiTIIMTHRUTOP*PAS5THR1ITOP>iDPNPRTPPVI 
VmPMFNTUMTHRURTGHT*MASRTHRURIGHT>«DOmORPIGHTVI 
FKirRGvTHRUTOP*MASRTHRllTPP>«OPNOPTPPRPFCTF  ICFwFRGYI 
Fk'TRGyTHRURIGHT^MARSTHRUPI  GHT  *PONOR  RIGHTS  pECTFirFKiFRGYl 
Ir  PFn«IMAX  THFN  PTH«-DTN-FNFRGYTHRIIBIGHTJ 
»  HPPATF  II 

UrJ]*(AMXtJl*UtJJ-llPOHFK'TllMTHRl'TOP*|iMOMrNTUMTHPIiPPTTPM- 

UMPMFNTUMTHRURTGHT*RTR( 1 »  »  UHOMFNTUmTHRIIRIGHT) jkpnfpvframxnfn j 

*  HPPATF  V 

Vr.l]«-(AMXtJl*VU]-VMnMFNTllMTHRUTOP«-vHPMFNTUPTHPIlROTTPM- 

VmPPFnTUMTMRUBTGHT*PTR(1»»VMOMFNTUMTHRI'BIGHT))hPNFPVFRAMXA'EWJ 

*  l'PP*TF  A I X 

AtytJU«AMXtJlw8PFClF!CFNFRfiY-FNFRGvTHP|lTOP*FNFPGYTMRUROTTPM- 
FAirPGYTHRURlGHT4-RTP(1»»ENFRP,YTHRIlRIf!HT>«0.RwAMXWF*lw(lirjlwMrJ]* 
VtJ)*VtJm*PNFOVFRAHXHFWl 
«  HPPATF  AMX 
AMXI JUAMXNFHI 

%  RTPrF  BCITTPM  OUANTITIFS  FOR  NFxT  POW 

UMPMFNTUMTHRIIBnTTPM«-UMOPFNTUMTHRllTPP| 

VMPMFNTUMTHP|lROTTnM*VHPMFRTllMTHRI)TnP| 

fnfrgythrubpttpm^fnfrGythrutopi 

MARSTHRHBOTTPM«.HARRTHRtlTnPI 

cndi  *  j  i pop 

FTH*FTM+SUMRPR(DTH)I 
FNOI  *  PH2 
8 
% 
t 
« 

SURRPUTINF  PUTPHTI  RFGlN 

CU  RFAI  TPTALFNFRGYHnW*RFLFRRNOR»FNGYCHrrK#NliHBFPPFrvci  rSl 

MODFMMPDFl 

NPA4-O.OI 

LOOP  ji.O,1,jMAX  OP  wPA«.wPA*AMXtJl*(ATX[JI*0.$>«Mir  J]*"t.lWV[  J]*Vf  Jl) 

)i 

TOTALFNFPGYNOR^SUMRnHtWPA )l 
RELFRRNPW*fTPTALCNFRGYPPW*FTH)/ETHI 
NuHBFRpFCVrLF5*CYCLr»CYCl  ASTPRINTl 
IF  (NUmBFUpfCYCLES-O’.O)  THFN  NUMRERPFC YCI  ESM .01 
ENGYCHFCK«-fRFLERRNtlW-RFl  FRRPFFORF  )/NUMBFP0FCYCLFS» 
RELFRRpFFnpF*RELERRNORi 
CYCLASTPRTNT*CYCLFI 
HOOF ♦TRUF I 
PUTCPRf CYCl F  #  Z# 1 )  I 
PUTCPRfT  »Z#?)I 
PUTCPRfPT  #Z#3)I 
PUTCPRfPTLAST»Z#«)l 
PUTCPRfFTH,Z#7)l 
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kh  > 


PUTCPRr<T«PirCTR,7»«)J 

PllTfPRrfiC.7.9)» 

PUTCPPfTP7pNF  *7*10)1 

PllTCPPf  P*  7  •  1  1  )IP|iTCPP(ft.?,l?)|PHTCPR(0»7, 1  3)  I 
IP  TPPRFHFCTIVC  THFN  PllTCPPC  1  *  Z#  1  !  3  * 

IF  PIGkTRFfI FrTIVF  THE>  PdTPPRf 1 ,7# 1 ?) » 

IF  PrTTPPRFFl FCTIVF  THrN  PUTfPP M . 7, 1 3 ) | 


PUTCPRf PFI FRPBFFPRF*7»14)j 
MM»PItF(  TAPISK1»7)| 

»IF (RFl  FRRnPW>MAXPF|  FRp)  THFN  FPPPRr?00>i 
*m»TFiiT . 


PATV  PROGRAM  STARTS  MFPF. 


SHFLI  t  PPFNf TA0|SK0#0# 1)  I 

npFKf'r«niRKi,o,n» 

Ik  RUT | 

ArM»'t  CDT! 

Output  i 
Phi  I 

RFGATT VFM« SSf  3)  JNFGATlVf FHFRRYf  3 ) t 
PMRI 

MFGATt  VFMASSf  «  1  »Nrr,ATlVFFMFPr,Yf  a  1 1 
»  AT  THIS  PPTMT  SHOULD  GO 
t  IF TPF70NFTPP1  THFH  PF7nvFfl)l 
t  Iff  pr  70NFPTGHT  )  TMfN  RF70WF(?)| 

*  TF rpr70NFROTTPPl  THFN  PF7PNFf3)J 
GP  TP  AGAIN! 

* 

STPPIFIWIFND. 
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APPENDIX  X 
SHELLN 


(CODE  LIST  DBUGA  SUMKY  SAVt 
BEGIN 

I  THIS  IS  IMAxjMAX  ZUNE  SHELL. 

BOOLEAN  MODEU 

Cl  NT  lMAX.jMX.jbLK, J, J1  ,  J2,  IF)Ll' .  IBLKP  1 » l,Ji 

CKEAL  PROb.CYCLE.T.DT.DTLAST*WCA»wCb.'WC(  .WCO.WCC.PI.ETH.STABLFCTR) 
CO  HEAL  RELERRBEF ORE->  C  YCL  ASTPR  INT  ,MA  XHEtERR.  GO.  TOZONE) 

PHEAL  CS.WPA.WRB.WPC,  OTH) 

PKEAL  VECTOR  0Y[0 J  .OX. RC . TAul 21 j 

BOOLEAN  BOTTOHREf LECT1VE,  KIGHTREF LECTIVF.  T OPRE FlEC T I VE) 
bOOLEAN  RtZONETOP. RE ZUNE RIGHT. RE  ZONE  BOTTOM) 

label  shell. again. stop. fin) 

PhEAl.  STDT)*  THIS  CARO  IN  OUTER  BlK  FOR  STD  TIME  COT. 

*  THE  FOLLOWING  CAROS  DEFINE  THE  COMMON, 

PKEAL  Li 

PKFAL  VtCTOR  U0MYlO],Ot2O].V0HYE03»Vt2OI»AMXDMYI0]»AMXC2O]» 
AlXOMYtOI. AlXtPO J»FDMY(0].Pl?0) ) 

PKEAL  VECTOR  XI2I.YI01# 

% 

FILE  14L1SN0“*,SHELL,,/',ST01"  (106  ROWS  FUI.  L  >  i 
FILE  14LISM«,,SHFLL"/"IT0B"  (106  ROWS  FOLD! 

FILE  LINE  (1  ROW ) i 

I 

% 

%  FOLLOWING  are  all  SUPROUT 1NES, 

1 

1 

sobrootine  hutcprccreal  x.preal  out  z.cint  d; 

BEGIN  IF  F’EN-I  THEN  Z*X)  END J 

t 

% 

X 

% 

SUBROUTINE  fUlCPKCIM  X.Ph'E  AL  OUT  Z.CINT  I>) 

BEGIN  CREaL  A i  A*X»  IF  PEN^I  THEN  Z*A)  ENOF 
1 

I 

1 

A 

PREAL  SUBROUTINE  LN  AS  RGA  (PR1AL  X  AS  RGA)) 

BEGIN 

S  INC1NS»100.  ASKRATCHbIB 
S*TAPE2=GLYPMR/GLN  SERIAL) 

END) 

* 

X 

X 

X 

PREAl  SUBROUTINE  E XP  AS  RGA  ( PREAL  X  AS  RGA ) ) 

begin 

X  INClNSelOO.  ASKRATCHslZ 
»»TAFE2«GLYPMR/GEXP  SERIAL) 

END) 

%' 

% 

% 

X 

PREAl  SUBROUTINE  sort  AS  RGA  (PREAL  X  AS  RGA ) 1 
BEGIN 

t  INCINS-1 1C.  ASKRATCH.H 
SSTAHE2BGLYPNIR/GS0RT  SERIAL) 
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END; 

% 

A 

A 

A 

PREAL  Su&KOlTINE  ALOGlOCPHEAl  X)F 
bLGU 

AL0GlC«-C.<t34?*Mfll9OJ*51tt*LMX)J 

EM'J 

A 

t 

% 

* 

PhEAL  SlbhdL'TlNE  PUrERCjFIO  (PHEAL  N )  J 
bLGIN 

1  PCikEhCF  1U=1U.**N. 

PUhl  P  OE  10*EXP  (2.302*6  509  29  ¥A0«6*N); 

END) 

A 

* 

t 

A 

SUBROUTINE  POMRCPKEaL  XxPKFAL  OUT  XX#P1N7  OUT  MU 
BEGIN 

XX«-ALbGlOCAESCX)  j; 

N*(xx-o.s*siGN(xxn;  adont  let  giypmr  round. 

IE ( N<1)  }  Then  N*N-U 

XA*XxPUkERLF 1U<-NJ; 

EXOJ 

s 

% 

) 

A 

SUEtprUTUE  ERRCRCClNT  J)i 
BEGIN 

S1PXH1TE  Cl 1NEIPAGE  J»«ERRUR  r"»J)> 
hlA*0» J 

hC A«- 1  *0  +  1 ,0/kCAI  AAbORT. 

EMDJ  AERRUK, 

* 

A 

A 

A 

SUBROUTINE  KE  GATIVEPASSCC1NT  J)» 

P  E  G 1 N 

C 1  NT  1I»J0J 

LOOP  JJ*0»1,JHAX  UO  bIGIN 
LOOP  11*0,1, IbLK  UU  bEGJN 
K*Jj*IbLKPl4l  l; 
lECAHMMSO.)  THEN  BEGIN 
S1MRRITE<IInEIPAGE],"NEG  HASS"»"J 
SIMkhlTE(LlKE»"APXl*l",AMXl*]); 

ERRUR(O); 

ENDIENUJENU; 

END;  *  NEGATIVEHASS, 

A 

A 

A 

A 

SUBROUTINE  NF GAT  I VEENERGY ( L 1  NT  U); 

BEGIN 
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CiM  1I#JoJ 

LuUH  oJ*0»l>JHAx  00  bLblN 
LOOK  1I*0,J»IBLK  UU  bt  G  I  N 
K* Jj* IHLKb 1 illj 
lKAlXEMSO.)  THE  N  BEGIN 

SIM|»RlTt(LlNE  l  A  A  G  E  J  *  •*  N  E  G  LNf  K"#  *  J  *"#0)) 
SlhM»ITl (L IkE  »"AlX(*]"*AIXE*J)J 
t  KPCjH  l  0  )  A 
EM'JE  NDJENOJ 

SlMimMTHL  InE#hMASS  AND  ENEAgY  DM1  J  =  "#0)J 
t  ND J  *  NEgATIVEEKEHGY. 

A 

A 

A 


A 

CKEAl  SlbhOlTiNE  SUMKL W < Hh t A L  >)J 
btf'IK 

CINT  wMbOTfJTHEoiBnuLLAN  C'LLMJllt; 
0LIAL0E*K1DE  J 

mjifmlt  muoe; 
x  4-  c .  i 

MILE  *  TRIE  ) 

Twl  TLTM  j«-|  j 
LULA  jM.l.b  UO 
BEGIN 

XM+RTAUMHOTRE  J#MODE  *>  >J 
HLIOTHE  J*ShlFTL(l»TACTC'THE  JU 

EM  ; 

SUE  RCA* Xi 


E  UL  E  ♦  t'L E  EiOtlF  l 
EMj  i  SOhApA 
* 

* 


* 


% 

CAE  AL  SLBAlHTiNE  CM  I Nl C  PE  AL  *.  lAlAl  Y  5  A 

bE  GIN 

CHIN*XJ 

IE  Y*X  lE’tN  CHIN*YJ 
ENDJ 
% 


% 

% 

% 

SLBPLL'T  I  Nt  INEUTJ 
BtGlN 

SIE AEAOl IflOlSKOf /)! 

PAGE)  «■  G  A  A  b  C  N  E  U*0); 
CYCLE«-GAAbCNElZ#n; 

CYCLASTE  AINTUYCLE  J 
7  4-GAAbLNE  <Z,2  )J 
01  4-GAAblM  l Z#  3) E 
DTLAS^GRAbtM  (Z*«>J 
IMAX  «-GAAbCiNf  EZ#b)J 
JMAX  «-GAAbONf  <Z#6)-1I 
E7H«-GAAbONE(ZM)J 
57  ABLE  CTR«-GRABOM (Z#6) ) 

I E ( STABLE C  7  PSO  t  )  THFN  STABLFCTRM  , ) 

GO«-GAAblNL(Z*V)j 

IMCCeO.)  THEN  G04980.) 

TOZlM  «-GRAbrNt(Z#10)J 
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IblM(lMXM)  D I  V  6<t; 

lbLKF  1  ♦  1  b L I*  ♦  )  j 

JblKtviFAX  OJV  6 6i 

UmFllLTllt  ♦  MtSti 
PlGFtlhFFLlcUVF  ♦  FALSE; 

(■umMUHtcim  ♦  false; 

IF  (GKAbLhE(2»ll)*U.)  THFN  1  TPPF  F  L[  CU  VF  ♦  TKUt  l 
IF  (GRA6LNE(7»12)FU.)  THEN  h  1  i,HT  PE  F  LFLT  l  VF  *T  KUE  i 
IKGPAfelNEC**  13)60.  )  IhFN  bt'T  T  UP  PF  F  L  E  f  T  1  VF  ♦  I  HUE  ; 

RE  l  E  Rhbt  F  0«f  *-GRAbl'NF  lit  14)1 
NAXRE LFNN*GRAbONF  ll,  15); 

ut-ouast; 

OCLE  ♦tVCLF -H 
LbOF  I*t»l»JPLK  DU 
ELGIN 

RF'A*hTR(  lrhliC’LrXtl  J)f 

IF(IXO)  THEN  IF  blJUltfcN(-O)  fcp  A  *■  R  7  K  C  1 »  #  >  f  1  - 1  J  )  j 
Kid  ]*0.5x(X(  1  1  +  NRA  )i 
uxn  j*x  1 1  j-xfa; 

TAL'L  1  ]»F1«()(  I  IxXC  1  J  -WPAXNF  A); 

end; 

hliDE  ♦tOLLtAF  (-O)i 

Dxto]*-Mui##Dxton; 
hcto 3 ♦ -h  T  n 1 ##PCtO] )1 
LOOP  .JUrlxjHlN  1)1) 
bl&lN 

IVlkil^HTl  CI#MUUF  »YTjJ); 
lFCoFublh)  TbFN 
E-  E  G 1  N 

IF  KULFAN(l)  THEN  0Y[JJ«  RTl  ( 1».  Yt  J+l  ] )  J 

end; 

m j]*iyui-yi j] ; 

end; 

END;  XlNblT, 

* 

% 

* 
t 

SOHRCl'T  I  Nt  E  5  (  C  I  M  I»L1M  J); 

*  MUDlF  It  AT  J  f!  N  OF  A1F<  TO  EXTENT  TD  2.F1J  EKGS/GP.  11  MAY  69 
BEGIN 

HNEAL  ErPF’U,FUVRO/!»lMRtlVKOx»PrwEN»FU»FUN»FN.rS»E  1 »  E  2#  riF  T  A  ,  F  E  »  G  h  0  N  F  i 

bUDLEAN  OLOMiLEI 

K  «■  Jx  (  1  bL  K ♦  1  )*1  ) 

t*l.(P-10xAl»lKJi 

J1*J  CIV  b<t; 

J2*J-64xjl ) 

oldpdde*mude ; 

MODE  +  TR0F ; 

*CA*1  AUxf.KApnNECDY  [  jl  ],  J?)J 

NUDEtULDMLDE i 
hFilUAPXl  J]/pt  a; 

ROVR02W73.3VS?xPF)O; 

LNFK:VRL2*LN(K0VPU2); 

PURlP*0.4j«?9A4blV032b2xLNEP0VPEZi *  PONFR* ALEI& 1 C( RLVR02 ) I 
El  *1 .Obbfcflixce.b-E); 

NS  ♦1.0/( 1,0+LXP( -(6.5-0, 3b7xbOWLR-t)/(U.97SxEXP(0.05xLNEPUVRU2)))); 

FD  ♦  hS  xEXP(-0.22*t2152xf  j; 

F0N*(1.0-NS  )xEXP(-0, 15C6295XE  ); 
b  E  T  A  <0 . ; 

IFCl>l.)TblN  BETA*  (0.04Bxhs  +  0 . 032x (  1 , -WS ) )xALOGl  U(  t )  ; 


Reproduced  from 
best  available 


copy. 


£i 
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Ek  ♦C.3J3J333j333J33ox(E-«0, ) J 

^ S«- 1  .0/(  1 .0*EXP<-0.?b*<f  -A5.0*EXP(0«0157xi  NEPOVROZ))* 

E  XF<-O.OBSxLNEROVROZ)>>l 
f  N*RS*EXf  l “0 • 0392 156*E ) I 
RF A« 30.O-fc.C*P0REHf 
If  RFA  >6.0  THEN  RPA*6,0J 
RPA*(E-160.C )/WPAf 
bETA*bETA*<  1 . 0-RS ) *0. 045*RS J 
E  E  *0 . 1 

If  RFA2-t>.0  THEN  f  E*1 .0/<  1  .0*f  XPC-RPA))  J 

GMONE  *<0. 161*0. 25b*Fti  *0.260*FON*0. 1  37xf  N*0,  ObOxFE  ) *EXP( RET A*l NEROVROZ  )  I 
P  l  K  J.GMLNtxRbDxAl  X  C  K  If 
CS*SfcRT( l 1 .O+oMONE  )xF[K]/KHfi)J 
END!  XEb. 

X 

» 

* 

* 

SlJHliCOTIf.E  O^IAJ 
BEGIN 

CODE 

USE  STDT.PIJ 

bEGIN 
HERE  IAN  I  GET  i 

ORG  FI  I 

LATA  3.1415V2653SB9793I 
OKG  STDTf 

DATA  l.l»1.2»1.3»1.4M.5»1.6»1.7»1.6»l.W2.0»2.1»2.3#?.b»2.6#?.6» 
3.0»3.2»3.«r4,6#3.e#A.0»A.2*A.S»A.tt»!>.0»5.3»5.6»6,O»6.b»/.()» 

7. 5»o. o» ». 5» 9. 0# y. 5* 10. Oj 
ORG  HE  PE  I  AM j 

FNr  code  J 

EnOJ  XOAIa. 

» 

» 

* 

X 

SODROCT INI  COlJ 
regie 

PhEAL  mindimfnsion; 

CREAL  BER.TnE  XT , RU A» RUB» *DY j  ClNT  LP»JJ  RDOLEAN  M1*M2j 
RD A* l . 0620| 

L00F  J*0»1#JMAX  DO 
bEGIN 

J1*J  D I  V  64F 
«2*J«6AXoi; 

RDT*GRAbONF(UTI Jl]#j2)f 
LOOP  I *0» 1 » I  PLK  00 
BEGIN 

N«-JxIBLK*1  l 
NODE 1*PEN*6A*ISIHAXJ 

If  (  I  PLK*0  )  THEN  MODE  1  *HODE  1  AND  NOT  bU0LEAN(-0)J 
HOOE*PODElj 

ESC  I » J)J*ES  RETURNS  UPOAlfD  P  ANO  CS  FOR  THE  ROE  K. 

X  I NO  I  KENS  I  ON* OX 1 1  ]J 

IF (RDY<MINDIMEN$ION)  THEN  HINDI  ME NSlON*ROYf 

NOB*  MlN<MINOIMENSION/CS)  J  *  TIRE  FOP  SOUND  TO  CROSS  MIN  CELL  DImN. 
NO  A*  C  Ml N( HO A»  ROB  )  J 

ROB*  MlM  OX ( I J/Ut  K ) ) IS  TIME  FOR  MATERIAL  TO  CROSS  CELL  IN  X  DIRECT  . 
FDA*  CMiMHOA»ROB); 

fDB*  MlN( RDY/Vt  K  3  > I  X  TIME  FOR  MATERIAL  TO  CROSS  CELL  IN  Z  DIRECTION 
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EDA*  LRlN(RDA#M'.b)l 


EM) 

END 

DT*0.5xSUblFCTF!xRDA) 

T*T*DIAST) 

MUi>.(7»bER*lP>) 

TMX7«-(1+DT}xRQ*EK0F  10(«LP)) 

Ml*  (TNEM>STET)  AND  EOOLEAM 1777777777776000COOOOU(»))) 


M  1  *  {  M  1  LK  M?  )  AND  Ntn  (Ml  AM)  P2>) 
II  (Ml  ANO  (BFh<STUT))  THf  N  J,  T  *(  ST  C  T 
UTLAsT*L7) 

CYCLE  *CKlE*1  J 


-hi  RJxRDM  ROF  10a.P)l 


LhOl  tCllf 


l 

» 

J 

» 


SUHKUTINE  PM) 

BEGIN 

C IP  T  MASS)  CREAL  G) 

PRfAL  VA»Vb»PA#PR»Pb#UHH*UHll 
D7HM.  ) 

IM  RlGHlRLMLKTlVt  THEN 
ERKOR(l) 

USE 

BEGIN 

X  SET  M  A  f  AXIS  AM)  HIGHT 
MULE  *EDGLE  AN(-0  )  I 
LOOP  J*G»  1 » JMAX  DO 
P[*)J*NTLII»»P[  JD) 

M00E*SH1F  TM(  IMAX+1  MiflPE  )J 
Ll)t'M  J*0»1»JMAX  DU  M(  J]*HTK(2»»M(  JJ)I 
1  nDJ 

LOOM  MASSM*!*?  DU 
BEGIN 

M  Q (> E  *TRUE ) 

IE  BOTTOMREE LECTIVf  THIN 
BEGIN 

Vb*0.) 

Pb*P  C  G  ) ) 

EM 

else  he&in 

VB*Y ( G 1 ) 

MCA*  GHAbpNE (OYlO)»G)) 

MB*  GRAbONE(DYlO)»l)l 

WCC*1  ./(hCA^HCB)l 

P6*(PI  1  Jxv.CA*P(O]xWC0)xwCC) 

DTH*OTH*Ot5xTAUXhCAxVtO]xoTx(P[C3  +  Ptl))>»LtC) 

EM) 

IE  RIGHT HEELFLTlVt  THEN  ERHCR(2) 

ELSE  BEGIN 

1  SET  U  AT  AXIS  AND  RIGHT. 

Hf L‘*bODLEAN(-0)) 

UHjF  j*D#1#JMaX  DO  UIJ]»NTL(1a*|i[J))I 
nor  E*SHIF  TP( IMAX+1 f NODE ) ) 

HCA*RYR(1»»NC>/KC) 

LOOP  J+U> 1 >  JH  AX  DO  U[  J)*NCA*RTR(  1»M»I  J])J 
END 

HOOE*NUT  bOOUAN(-l)  AND  (IMAXfcFENJF 
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NFaM.O/(DX  +  RTLC  W»DX))I 
WPL*Rt*KTl  ( l##DX]xhPA J 
WPC«-hTUl»/RC3*UX»«WPAF 
WCA*GKAbONE  (OYlUlfO)) 

LOOP  JMAX-1  I/O 

BEGIN 

IK  ArtXljISO,  ]  THEN  FRRUROIJ 
URR*L'l  J3*RPE+RTLl  l##Ul  J]  )*WFC  J 
PRKPt  J3*KTLU»»DX)*RTL(1»#H  jl  :  *0X  )*NPA  I 
Jl«(J+i>  Ol V  64 J 
j2*(  J  +  l  )-64*JlJ 
WCb«-CRARONE(UYt  Jl]»  J2>J 
WtCM.0/(NCA  +  «CB31 
►CD*HtA*WCCf 
WCE*wtBxNtCj 
VA*VU3*NLE*VC  J*1  JxWCUl 

AIXlj]*AIXl  J]*P[  JJ*DT*(0.5*(VR-VA3»«TAli*(BTR(l#»URP>*yRR>*PI*RCA3/AMXl  JJ) 
IflPASS.l)  THEN 
BLGIN 

PA*Pl  j3*WCE*P[JM  1*bCnJ 
6*6. 3567P8I 

G*G/(G*^RABONF(Y,J1]#J2]”0«5i«WCA]I 

G*G0»G*GJ 

Ulo3U'lo3*TAU><WtA><lHTR-<KAPK3-PR3*DT/lUX>cAHX[j3>J 
VI  J]*VlJ]*TAUxFtA*(PB-PA  3 *UT / ( »C A* AHX [ J 3 ) -GxUT J 
Ft  *P  A  J 
ENDj 
Vb*VAl 
NCa^wCRI 
t NO!  i  J  LOOP 
J*JR AX J 

IF  (AFX  t  v)  3  <  0  •  3  THEN  ERPOHUII 
URR«-UlJ]*WPb  +  RTLCl#>U[J3)>«wPCJ 

PK*(P[  J3"ML(1##OX)*RTL(1»aFI  j3 3*OX >*NPAI 

IF  10FREFLECTIVE  THEM  ERR0R(5] 
else  BEGIN 
VA«V[v>3! 

%  PA*PBJ  *  ALHFAOY  oone  abovf 

0TH*0TH-0.5*TAU*WCb»«Vl  J3*CT*lPl  J3*PCJ-1 3  3*WCCI 
AIXl J3*AI XI J3+PI Jl*OTxl0.5*<VR“VA3*TAU*(RfRll##URR)*URR3>«PI*NtA3/AMX[ J3J 

If  F  A  S  S  *  1  THEN  Ul  J  3 «- U 1 0  J*T  AU*WC  A* l  HTKC 1  *  •  PH  >*PR  3  xOT /t  OXxAHX  C  J 3 1 1 
»  DONE  WITH  WHOLE  uRID 
END!  %  PASS  LOOP. 

ETF'*,E1H*SUHNOw(UTH3J 

ENLFXHHI 


l 

t 

% 

l 

SUBROUTINE  PH2F 

Bt6INpE  RE A l  SPEC IF'ICENERGYaOONORTOPUaOONORTOFVaOONORTOPSPECIFICENEHGYa 
VAl TOP  INTERFACE aDONORHASSOI VI OEOB YD Y# HA SSTHRUTOPaMASSTHRUBOTTOMa 
IONORRIGHTUaDONORRIGHT  VaDONORRIGHTSPECIF ICENERGY, 

D0R0RHAS50I VI  Of UbYDXfUATRIGHTl NTFRF  ACEaMASSTHRURI GHT» AMXNEHa 

ONE  OVERAHXNEWaUHOMFNTUHTHRUTOPaVMOHENTUMTHRUTOPaUMOMENTUPTHRURIGHT t 

vhlhentohthruright,energythrurightafnergythrutop.fnergythrubottoha 

UHOHENTUHTHRUbOTTOM^VHCMENTUHTHRUBUTTOH; 

REiONETOP*F  ALSEJREZONERIGHT*f  AL5EF REZONEBOT TOH*F ALSF l 
DTH«-0.  j 
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NCA«-&hAbDNE<UYlO]»0)l 
1  it  1  AXIS  bUUNUxRY  CUfUlTP'NS 
K CL  E*btiOLEAN<’0)J 
►  ASS1HRI'K1GHT*0.) 

%  HI  bOlTClM  BOUNDARY  CONDITIONS 
F  OL  E* 1  RUE ) 

If  bonOMF.Ef  LECT1VE  ThEN 
bEGIN 

HASSlhMbOTTUM*U.I 
f  MRGYTHhUBD1TON*O.I 
lMLNENTUMThRUbOTTOM*O.I 
VMt  MENTUMTHPUBOllGMMtJ 

EM 

L  L  Si  bEC,IN 

PASSTHRUbUlUiM^VlOJxAPXtUJxnT/HCA) 

If  (NAsSYHRUbUYTliM>Ot)  THEN  MASSTHRUbOT  T0«*0 . 1 
If  <«HASSTHRuB01TUn>T020fE*AMXlOJ)  THEN  REZONF  UOTTPMURUE) 
tM  PGYTHRUbGTTUM*>  ASSTHRlJHGTIUMxt  A I  X  [  0  HO,  S*  (  U 1 0  ]  *0 1 0  ]♦  V  t  C  3* 
V[OmJ 

DTH*ENEKf.YTHRllb01  TOM) 

UHUHENTUMTHRUHOTTPM^HASSTHHUHUTTOHxUEO]  ) 

VPGME  NTUMTHNUBOT  TPM«-MASSTRKUBUTTOM*V  10 1 ) 


If 


E  NU 

RIGMREH.F  CTI  Vi  THtN  ERRf  H  (  6  ) 
fLSE  bEtIN 

X  SET  RIGHT  BOUNDARY  CONDITIONS 
FULE*SHlFTR<lMAXM»MGUF)) 

LOOP  j«-0»  1 »  JM A  X  DO 
bEGIN 

UI  J]*RTN(l#,Ut  JUA 
If  (UCJ]<0.)  THfN  li( JJ*-UU>I» 

f  NU 
ENU 

PODE*NOT  HOOLFAN(-i)  AND  <IMAX*PFN)I 


LOLP  J«-0#i#JMAX  DO 


bEUN 

i  COMPUTE  OUANTITIES  AT  TOP  INTERFACE 

SpECIFIUENERGY«-AI>l  Jl^O.SxtUt  JIxUl  Jl  +  Vl  JJxVl  J])» 

UUNURTOHUMI  J]f 
DnNJRTOHV»Vl J] ) 

OUNORtOP SPECIF  ICE NERGY*SPEU I FICENi ROY J 
DONORMASSDlVinEOBYDY«-AMX[  Jl/NCAJ 

IF  J< jMAX  THEN 


bf  GIN 

D I V  64) 

J2« ( J+ 1 J-64XJ1) 
hCB^GRAUONE  (OY  U'l  j*J2)i 

VATT0P1NTERF  ACF  «-0«5x(V[J]*V  l  1  1  )) 
If  (VATTtlPINTEHF  ACr<C.)  THEN 


OEUIN 

DONnRMASSl’IVIDEOBYDY«-AMX(J]/NCBj 

UONORTGPUMlJ+lJ) 

DONORTUP Y* V C  J+ 1 1 )  .  ,  ,  , 

DUNOHI  OPSPEC  If  ICENEKGY*AIXU*n*0.5x(Ut  J*l3xu[  J*1 
V l J*i JXVEJ*  1 1 )  T 


]♦ 


MASSTHRUTOPnATTOPlNTERFACExDTxDONORMASSOIVIDEDBYDY) 


IF 


end 

ELSE  BEGIN  X  00  TOP  RON 
TOPREFLECTIVE  THEN  ERRLHC7) 

else  begin 


OIF  FEKENTL Y • 
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VATTOPINTLRE  ACf  *V[j]J 

MASSTHKUTOP«-VATTOPlNURFACE*UTxOONQRMASSDlVIOECBYDYl 
IF  (MASSTHRUTCIP<0.  )  THEN  MASS'  HHUT OP^O*  I 
IF (MA$STHRUTOP>TOZONExAMX[ JJ)  1  HEN  REZONl  TOP«-l  RUEl 
OTH«-DTH*MASSTHRUTOPxDUNOKTOPSPECIFICLNFRGYI 

END' 

end; 

J  COMPUTE  QUANT  I  TIES  AT  RIGHT  INTERFACE 
I UNORRlGHTU*Ut JJi 
OONORRIGHTV*VE JII 

OONORRIGHTSHECIF ICENERGY*SPECIF I  CINE  RGYi 
iifNORNASSOI  VIDEDBYUX«-AMXt  J  J/C  L'XxRC  )  J 
UAThlGHl  INTLRFACL*0.5x(UIJ1*RU(1#»U[  JJJJJ 
IF  CUATHIGHTlNTEHf ACf <0. )  THEN 
BEGIN 

DONURMASSOI V  I  [if  06  YOX*HTL  ( 1 # * DUNORMASSDI V I OFDBYDX ) J 
OONURHlGHTU*RII  ( 1 »  » CONOR K I GHTU ) I 
OONORRIf,HTV*RTL  ( 1  #  #  DONPRK  I GH1  V  )  I 
DflNUPRIGHTSPECIFICENEMGY* 

rtl(  j*»upncihright  Specif  icfnergyij 

E  NO  f 

MASSTHRUHIGHT4-UA1RIGHTINTEHF ACFmDTxOUNORMASSOI VlOEUbYUXxXi 
I;  (GRABUNF(MSSTHRURIGHT»IMAX)>TUZ0NE*GRABUNECAMX[  JI»IMAX)) 
THEN  REZONERIGHUTRUFJ 

APXNEN«-AMx[  Jl-MASSTHRUTOP  +  MASSTHRUHOTTOM-MASSTHRORIGHU 
RTPt 1»»MASSTHRUHIGHT )l 
UfEOVERAMXNEW*! .U/ANXNf Wl 
*  MOMfNTUM  FLUXES  AT  TUP  AND  RIGHT 
OHn«ENTUMTHhUTUP*MASSTHRUTUPxOONORTUPUI 
UMPMENTUMTHRURI GHT^NASSTHRURI UHTxUOnORRIGHTuJ 
VMOMFNTUMTHRUTUH«-MASSTHRUTOPxUONOHTOPV> 
VM0mEnTUMTHHURIGHT«-NA5STHRURIGHT*D0N0RRIGHTVJ 
ENf hGYTHRUTUP^MASSTHRUTOPxUONURTUPSPEClF ICEnERGYJ 
ENEHGYTHRiiRlGHT*MASSThRURlGHTxpQNUNRIGHTSPECIF ICENFRGYl 
If  PEM1MAX  THEN  OTR*UTh-ENERGYTHKURIGHTJ 
%  ikoate  u 

UfJ]*(AMX[JJxUtJj-uROMENTUMTHKUTUP+|iMOMENTUHTHRUBOlTOM- 
OmPHENTUHTHRURIGHT^RTRC 1 #>UMOMFNTUHTHRURlGHT) )xONEUVEHAMXNrR> 

%  UPOATE  V 

V[JJ«-(AMXfJ]xV[jJ-VMGMFNTUM*HNUTOPPVMOMENTUMTHRUBOTTOM- 
VMDHENTUMTHHURIGHT+RTR< 1»» VPOHFNTUMTHHURIGHT) IxONEOVERAMXNERi 
X  UPOATE  A I  X 

AlXtJ]«-(AMAtJJxSPFCIFICENENGY-ENERGYTHRUTOP*ENERGYTRHUBOTTUM- 
ENFRGYTHHUHIGHT+HTR( 1»#ENFHGYTHRUHIGHT)-0,SxAMXNENx(UI JJxUl JJ+ 
V[ JJxVt JI nxONEOVERAMXNEW) 
x  update  amx 
APX [ J ] ♦ AMXNEWI 

X  SI  ORE  BOTTOM  QUANTITIES  FOR  NEXT  RON 

umcimentumthruruttom*umomentumihruiopi 

VM0MEnTUMTHHUR0TT0M*VM0MENIUM1HRUT0PJ 
FNFRGYTHRUBOTTOMx-FnFRGYTHRUTOPI 
MASSTHRUHOITOMi-MASSTHRUTOPI 
ENui  %  J  LOOP 
ETH«-LTH*SUMRUN(01H)I 
ENDl  *  FH2 

s 

i 

% 

% 

SOBROUT  1NL  RF.ZONE  <  C  I  NT  J)J 

blgin 
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ERROH(IOO); 


END) 

> 

X 

X 

X 

SObROUTiNE  OUTPUT  f 
BEGIN 

CO  REAL  IoTAIF-NE  RGYNOrt# PF LEBHROHjf  NGYCHtCtuNORBEROF  CYCLES! 

HUDt«-NU»  bOCLEAN(-O)  AND  PtRSlMAXj 
RKA«-G«vl 

LOUR  J«-0/>l>jHAX  DO  HP  A«-HPA  +  AMX  t  J  1*  (  A  !  X  (  J  )*0 . 5x  ( U[  J]  *D[  J]4  V  t  J]  x  V  [  J]  ) )  ; 
TOTALED  RGYnCR«-SOMHOW(*»PA)I 
RtLERRNUN«-(T01  ALtNERGYNOH-ETHT/ETRl 
NuMbf  RUKYCLES4-CYCLE-CYCLASTRRIKT) 

If  CNUMbtRUFCYCLESsO.O)  TREE  NUR6ERUF  CYCl  E  S «- 1  ,0! 
ENGYCRELK«-(bELERRNUR-htLf  HNBEFIIKE  I/NOMBERfiFCYCLLSI 
RELERRBtFORf  «-hELE  KRNOhl 
CYCLaSUR1NT«-CYCLE  J 
MllDEONOE; 

POTCPRCLYCLf  tit  1 ) ! 

PDTtpRCT  tl> 2)f 
PUTCFRCOT 

PDTtPR(LTLAsT>Z,4); 

POTCPR(ETh,7,7); 

PUTCPR(STAblM.TR»Z»8); 

P0TCFR(i>0,a,9); 

PDTCFR(10<UNF*Z,1U); 

P0TCPR(C,t»iniPUUPR(0»Z»l?)IPUTCPR(0,Z,13)J 
If  HPRtFLECTm  THEN  PUTCPRT  1  »  ZM  1)  J 
If  RlGMRtfLF  CTOVF  THEN  PUTCPPt  1  >  Z  >  1  'i )  I 
If  BOTTORbff IECTIVE  IREN  PUTCPPt 1»Z»13)J 
R0TCFR(RELF  RRREFOREW.  1A)J 
SlMhRlTEr  U'ISKW); 

*Jf(RHtRRNOR>MAXRELEKP)  TRFN  EPR[lH(  200 )  J 

END;  xootput. 

X 

> 

X 

X 

X  RAIN  PROGRAM  STARTS  HFRE, 

X 

X 

SRELLI  OPENt  HD1SK0,0,  1)1 
OPEN(  I  4D  I  SK  1  *  0*  n; 

inpi.i1; 

NEGATIYEMASST  niNtGATIvEENERGYd); 

AGAIN!  COT) 

ooyfut; 

negativlmass(<:);nlgativeenergy(2); 

PR1; 

NEGATIVE  MASS ( 3)|NEGAT1VEENERGY(3)| 

RR2; 

negativemass(«);negativeenergy(«); 

X  AT  Tills  POINT  SHOULD  GO 
X  lF(REZONETUP)  THEN  HEZONE(l); 

X  If  CHI ZONERIGHT)  THEN  HFZ0N£(2)| 
t  IF(REZUNEBUTTOR)  THEN  REZONE(J)) 

GO  TO  again; 

X 

STOPif  IMEND. 
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APPENDIX  XI 
SHELL/OF/THE/FUTURE 


DBUGA  LIST  ZIP 


SHELL /OF/ THE/ FUTURE. 
THIS  IS  IMAX*uMAX  SHELL 
R.W.  CLEMENS  14  MAY  71, 


WRITTEN 

REVISED 


FOR  OUT  OF 
24  AUG  71. 


CORF  JOBS. 


CURRENTLY  DIMENSIONED  FOR  A  MAX 
50  ROW  BLOCKS. 

NB  A  GIVEN  I  KAY  (WILL)  BE  MORE 


OF  638*1280  IN  A  MAXIMUM  OF 
THAN  ONE  OF  THESE  ROWS. 


I  AND  J  REFER  TO  THE  UPPER  RIGHT  CORNER  OF  A  CELL. 

THERE  IS  NO  OUMMY  ROW  AT  THE  TOP  OR  BOTTOM  OF  THE  GRIO. 
THERE  IS  A  DUMMY  COLUMN  AT  THE  LEFT  ANO  RIGHT  OF  THE  GRIO. 


SSUMRY 
BEGIN 
X 
X 
7. 

/. 

X 

/. 

7. 

7. 

7. 

7. 

X 
X 

C IN  T  IK AX, J MAX ,I,J,K,J1,J2,I3LK, IBLKP1 , IROWS , IRCWSM1, JS TRIPS »JBOT , 
JMAXPERSTRIP, JPL, JPN, JB, JC, JJ, IROWS2, 1  ROWS 3, IR0WS4, JSTRIPSP1 , 

CREAL  PROB, CYCLE, T,DT,HCA,WC8,WCC,WCD,WCE,WCF, PI, ETH.STABLFCTR, 

RELERR8EF0RE,C YCLASTPRINT , MAXRELERR,GO , TO ZONE ,TOUMP ,CYCLESTOP , 
WCAjTOTALENERGYNOWJ 
PINT  IR,IR2,IL; 

RREAL  Z,CS,WPA,WPB,WPC,DTH,STDT ,PA,VA,ZSALUTEJ 
BOOLEAN  0OTTOMREFLECTI VE ,RI GHT REFLECT  I VE, T  OPREFLECTIVE,MOOEl , 

F EZONETOP, REZONERIGHT,REZONEBOT TOM, STOP THISCYCLE? 

BOOLEAN  OLDMO DE ?  7. V. 7. X 7 X 7 7. 7. V. 7. 7. 7. 7. 7. 7. 7. X 7 X X 7. 7. 7. 7. 7. 7. V. V. V. V. X 7. V. 7. 7. X 7. X 7. 7. 7. 

BOOLEAN  VECTOR  SWI6JJ 
LABEL  SHELL, AGAIN, STOP, fin; 

y.  U,  V,  AMX,  AIX,  AND  P  MUST  BE  DIMENSIONED  TO  AT  LEAST 

X  2* (BLOCK  SIZE>-1. 

7. 

PREAL  VECTOR  U, V , AMX, A  IX ,P( 99) ? 

7. 

7. 

X 
7. 

PREAL  VECTOR  L AST , NEXT  I  249 1 ; 

X 
X 

ym 

PREAL  VECTOR  X , DX , RC, TAU, VB,PR, PB, URR, WP1 , WP 2, WP3, 

MASSTHRUBOT  TOM , ENERGY THRUBOT TOM, UMOMENTUMT HRUBOT TOM, ENERGYTHRUTOP, 
LKOMENTUMTHRUTOP  ,VMOMENTUMTHRUT OP  » UMOMENTUMT HRURIGHT , 
VM0MENTUMTHRURIGHT,ENERGYTHRURIGHT,30N0RRIGHTU,D0N0RRIGHTV, 
DCNORRIGHTSPECIFICENERGY , MASSTHRURIGHT , MASSTHRUTOP , 
CCNORMASSOIVIOEOBYOX ,VMOMENTUM THRUBOT TOM, SPECIFIC ENERGY! 9) i 
BOOLEAN  VECTOR  IMO DE , IMODELI 9) 5 

y*  the  FOLLOWING  ITEMS  MUST  BE  OIMENSIONEO  AT  LEAST  JMAX  (SIC). 

7. 

CREAL  VECTOR  Y,OY[1280i; 

X 

X 

I4IN0UT  (BUFFERS=0,ACCESS=RANOOM,PHYSICALRECOROSIZE*250  RONS, 
TITLE=*SHELL/I4IN0UT.*)J 

I4SCR  (BUFFERS*!) , ACC ESS* RANOOM, PHYSIC ALRECOROSIZE* 250  ROMS, 

t  it  le=*  Shell/ i4S(Tft.  * ) ; 

I4SALUTE  (BUFFERS*0,ACCESS=RANOOM,PHYSICALRECOROSIZE=1  ROW, 
TITLE=+SHELL/I4S ALUTE.*) J 


THE 


FOLLOWING 
5* (BLOCK 


ITEMS  MUST 
SIZE)-1. 


BE  OIMENSIONEO  TO  AT  LEAST 


THE  FOLLOWING  ITEMS  MUST  BE  DIMENSIONED  AT  LEAST  (IMAX«-1)/E4. 


XI 

FILE 

XI 

XI 

FILE 

XI 

XI 

FILE 

XI 

11111 

11111 

11111 

11111 

11111 

11111 
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m 


StHESi 


I 


FILE  l4JISKO-*SHELL/I4CLAf1/'OUT6QBY15+  (250  ROMS  FULL)  5 
FILE  I4DISK7=  +  SHELL.  /  I4SCR*  (250  ROMS  FULL)? 

FILE  LINE  (1  ROM) ?  . 

X 

X  FOLLOMING  ARE  ALL  SUBROUTINES. 

x 

y. 

PKEAL  SUBROUTINE  ALOG  AS  RGA  (PREAL  X  AS  RGA) ? 

BEGIN 

1  INUINS=100,  ASKRATCH= 18 
SSTAPE2=GLYPNIR/GLN  SERIAL? 

END?  XLN. 

V. 

V. 

V. 

X 

PREAL  SUBROUTINE  EXP  AS  RGA  (PREAL  X  AS  RGA)? 

BEGIN 

$  INGINS=103,  ASKRATCH=17 

81tape2=glypnir/gexp  serial; 

END?  XEXP. 

7. 

7. 

7. 

7. 

PREAL  SUBROUTINE  SORT  AS  RGA  (PREAL  X  AS  RGA)? 

BEGIN 

t  INCI ’S=ll 3  »  ASKRATCH=14 
SiTAPE2=GLYRNIR/GSQRT  SERIAL? 

ENC?  XSGRT. 

X 

7. 

7. 

7. 

PREAL  SUBROUTINE  ALOGIO (PREAL  X)? 

A LOGlu=0. 43429 44819 Q32518*ALOG (X) ? 

7. 

7. 

7. 

7. 

PREAL  SUBROUTINE  POWERCF10  (PREAL  N)? 

POMtR0F10=EXP(2.302585092994046*N)  ?  7.  POMEROF10=10.**N. 
7. 

7.  ...  ... 

7. 

7. 

SUBROUTINE  POMER  (PREAL  X, PREAL  OUT  XX, PINT  OUT  N)  ? 
BEGIN 

N=TRUNCATE(ALOG10 (ABS(X) ) ) ? 

IF (N<0)  THEN  N=N-1? 

XX=X*POMEROF10 (-N) ? 

ENB ?  7. POMER. 

7. 

V. 

X 

7. 

SUBROUTINE  ERROR (CINT  J) ? 

BEGIN 

XI  EXIT ( J) ? 


11111 

11111 

11111 


11111 
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S IHWRITEi LINE! PAGE 1,+ERROR  =  +  , J)} 

WCA=0. J 

WCA=1.0+1.0/WCA;  '/.ABORT. 

ENOi  XESROR. 

'/. 

'/. 

7. 

SUBROUTINE  DATA; 

BEGIN 

CODE 

USE  STDT,PI,IR,IR2,IL; 

BEGIN 

HEREIAM  SET  ; 

ORG  pi; 

DATA  3.141592653589793, 

ORG  stdt; 

DATA  1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2. 0,2. 1,2. 3, 2. 5, 2. 6, 
2. 8, 3. 9, 3. 2,3. 4, 3. 5,3. 8, 4. 0,4. 2, 4. 5, 4. 8, 5. 0,5. 3, 5. 6, £.0, 
6. 5, 7. 3, 7. 5, 8. 0,8. 5, 9.0,9.5,10.0; 

org  ir; 

DATA  (0)63, -i; 

ORG  IR2J 

DATA  (0)62, (-1)2? 

ORG  11} 

DATA  1,(0)63; 

ORG  hereiam; 
end  code; 
end;  xoata. 

X 

7. 

7. 

7. 

SUBROUTINE  ES (CINT  K) $ 

BEGIN 

7. 

X  THIS  IS  11  MAY  69  VERSION  OF  AIR. 

7.  IT  WILL  GIVE  SLIGHTLY  DIFFERENT  ANSWERS  THAN  THAT  CODE 

X  AS  HERE  WE  COMPUTE  THE  EXPONENTIAL  SWITCHS  IN  ALL  CASES. 

X 

X  ES  USES  RHO  [IN  CS  ARRAY),  AND  AIX(K). 

7.  IT  RETURNS  UPDATED  PI  K 3  AND  CS. 

7. 

P REAL  e,lnerovroz,fo,fon,fn,ws,beta,fe,gmone; 

FREAL  RHO  AS  CS; 
c=1.0+-10*AIX[Kl ; 

LNEROVROZ»ALOG(773.39520495*RHO>; 

WS=t.0/(1.0+EXP< (-8. 5-0. 155043 14*LNER0VR0Z+E) *EXP (-0 . 05*LNEROVROZ 
+0.02531780798) )) ; 

FC  =WS  *EXP(-0 , 22421 52 4664*E) * 

;-'ON=(1.0-WS  )*EXP(-0.15082956259*E) ; 

E£TA=0 . ? 

IF ( E>1 .) T  HEN  BET  A*  <0 . 0069487*WS+0 . 0 138974) *ALOG <E) ; 
WS=1.0/!1.0+EXP( (-E+EXP(0. 0157*LNEROVROZ+3. 806662 489)) *EXP(-0.C85 
•LNEROVROZ-1. 38629436) ) ) ! 

FN=WS*EXP (-0. 0  39215686275*E) ; 
e£TA=BETA*(1.0-WS)+0.045*WS; 

WPA=3  0 .0+3.47  4356*LNER0VR0Z; 

IF  WPA  <6.0  THEN  wpa=6.o; 
wpa=(160.o-E)/wpa; 


11111 

mu 

urn 
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KNKXMO  ri/iUXX 


r 


FE31.0/(1.0*EXP(HPA)); 

C-MONE3  (0  „161+0.255*FO+0.280*FON+0.137*FN+0.05*FE)*EXP(8ETA 
•LNEROVROZ)  5 
P  [K]=Gi10NE*RHO*AIXt  K1 » 

CS“SQRT ( ( 1 • 0  +GM0NE ) *P[ KI /RHO) } 

END?  X E S . 


UBROLTINE  COT » 


GIN 

preal  hindimension;  creal  hob; 

LCCP  I=0,1,IBLK  30  BEGIN 


*=i; 

LOOP  JJ=1,1» JMAXPERSTRIP  DO  BEGIN 
J=JJ+J30T* 


M03E=IM0DEI I I » 

CS  =  AMXtKJ/(0Yl  JI*TAU[in; 
ES(K) ; 

MINOIHENSION=Oxm  J 
IF  3YI J XMIN3IMENSI0N  THEN 
hob*hax (Cs/minoimension) ; 
IF  W3B>H3A  THEN  W3A  =  W3B; 
WDB  =  HAX(ABS(UTKI)/OXm)  5 
IF  WDB>HDA  THEN  W3A  3  WDBJ 
hdb=hax (abs(vck])/oy[J]> ; 
IF  H3B»H3A  THEN  H3A  3  W3BJ 
K=K+IBLKP15 
end;  XJLOOP. 
end;  xiloop. 

03E=true; 

imwrite(line,*cot-hoa*,woa); 


MINJIMENSIONOYl  J] $ 

X  INV  TIME  FOR  SOUND  TO  CROSS 
X  MINIMUM  CELL  3IMENSION. 

X  TIME  FOR  MATERIAL  TO  CROSS 
•/.  CELL  IN  X  DIRECTION. 

X  TIME  FOR  MATERIAL  TO  CROSS 
CELL  IN  Y  DIRECTION. 


NJ5  X  CDT • 


*/ 

SUBROLTINE  CHOOSEOT; 

BEGIN 

CREAL  REG.TNEXT;  CINT  LP;  BOOLEAN  Ml,M2,EVERYSTOTIME5 

V. 

X  THE  BOOLEAN  EVERYSTDTIME  CHOOSES  ONE  OF  TNO  OPTIONS  FOR 
X  COMPUTING  STANDARD  TIMES.  THE  TNO  OPTIONS  WILL  DIFFER  ONLY 
DURING  THE  FIRST  FEH  CYCLES  OF  A  CALCULATION  WHEN  IT  IS 
X  POSSIBLE  FOR  T+DT  TO  VAULT  ACROSS  SEVERAL  STANDARD  TIMES. 

X 

tVERYSTDT I ME3 TRUE  t 
jt=o.5*stablfctr/hda; 

*iJA  =  1.0t-20J 

t=t*dt; 

POWER ( T  »REG»LP ) ; 

TNEXT=(T+OT) *POWEROF10 (-LP)  5 
IF  EVERY ST3TI ME  THEN  BEGIN 

X  THEN  CLAUSE  FINOS  EVERY  STANDARD  TIME. 

Ml3 (REGSSTDT)  OR  BOOLEAN (1/77777777(0)); 

M2=SHIFTR<1,M1) ; 

Ml3 (Ml  OR  M2)  ANO  NOT  (Ml  ANO  M2) 5  _ 

IF  Ml  AND  (TNEXT>ST3t)  THEN  BEGIN 
OT=(STDT>REG)*POHEROF10(LP)  *, 


modestrue; 
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mmrmmtimm 


swww 


mxmimmmmm 


SIMHKITE(LINE,*OT+,OT>  ? 

TOTALENERGYNOH=0; 

T0UMP=CYCLE+2? 


7.7.7. 
Z.’/.7. 
7.7.7. 
7.7.7. 
7. 7. 7. 
7.7.7. 
7. 7. 7. 
7. 7. 7. 
7. 7. 7. 


7.7.7. 
7.7.7 
7.7.7. 
7. 7. 7. 
7.7.7. 
7.7.7. 
7.7.7. 
7.7.7. 
7.7.7. 


end; 

7. 

7. 

7. 


HERE  PULSE  OSK  TO  LET  IT  KNOW  THAT  A  DUMP  IS  COMING  AT  THE  END 
OF  THIS  CYCLE.  THIS  SHOULD  GIVE  IT  TIME  TO  LOAD  THE  CORRECT 
STRIP  TO  LASER  STORE. 

FORMAT  OF  OSK  STATEMENT  IS  OF  YET  UNDEFINED. 


end; 

ENO 

ELSE  BEGIN  ...  ..  _ 

V.  ELSE  CLAUSE  FINDS  THE  LARGEST  STANDARD  TIME  LESS  THAN  TOT. 
Ml=  (TNEXTiSTDT)  AND  B00LEAN(1777777777776Qr' 1000000 (8) ) 5 
M2=SHIFTL  Cl, Ml) } 

Ml= (Ml  OR  M2)  AND  NOT  (Ml  AND  M2)J 
IF  Ml  ANO  «REG<STDT)  THEN  BEGIN 

DT=(ST JT-REG)*POHEROFlQ(LP>  i 

totalenergynowso; 

TDUMP=CYCLE+25 


HERE  PULSE  OSK  TO  LET  IT  KNOW  THAT  A  DUMP  IS  COMING  AT  THE  END 
OF  THIS  CYCLE.  THIS  SHOULO  GIVE  IT  TIME  TO  LOAD  THE  CORRECT 
STRIP  TO  LASER  STORE. 

FORMAT  OF  OSK  STATEMENT  IS  OF  YET  UNDEFINED. 


end; 

end; 

CYCLE=CYCLE+i; 

stopthiscycle=false; 

IF  CYClEJCYCLESTOP  OR  SHI  4 1  THEN  STOPTHISCYCLE=TRUE5 
3CCHOOSEJT. 


■/. 

SUBROUTINE  SETDXDYETC5 
BEGIN 

MODE=NOT  BOOLEANC-O)? 
LOOP  1-0,1, IRLK  DO  BEGilN 
IMODElI)=MODE5 
MOOE=TRUE{ 

IMO DEL  II l*MODE5 


end; 

IMODEIIBLKJ*IMOOEI IBLK] j»ND  IMAX>64*IBLK*PENt 
I  MODEL  1 1 BLK 1= I  MODEL  1 1  BLK  1  AND  IMAX>64*IBLK«-PENf 
LOOP  I=0,1-IBLK  DO  BEGIN 
MODE=iMODEI 1)5 
WPA=RTR (1 , , XI I+IR1 ) 5 
Rcm=o.5*(xm»wPA)  ? 
oxt i )=xi i j-wpa; 

TAUt  I)=Pf*"<xril*xTIl-¥PH*"wPA) ; 


end; 

MODE-BOOLEAN  (-0) ? 
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end; 

x 

x 

v. 


cxm=RTL<itiOxm> ; 

«C[a]*-RTL(l»,RCtO)>} 
M0Cc=REVR(IMAX*1-£4MBLK,M00E)  } 
JXtIBLKl=RTR<2,,3XC  IBLKMR2]) 5 
RC  t IBLK ) =RTR ( 1 »  »X  1 1 BLK*IR I ) *0 .5*DXt I^LKl » 
LOOP  1=0,1, IBLK  DC  BEGIN 
MODE-IMODELCIlf 
WPA=RTL (1 < »OXt I+1L1 ) 5 

hpiiii*i.o/oxiii*mpa>; 

HP2C  Il=HPA*RCt  H*HPlt  I)  f 

WP3[I3=RTL<l,,RCU*ILl)*3XIIl*WPim; 


enc; 

M0JE=TRUE; 

LCCP  J=1 » 1 , JM AX  00  DYl Jl  =  Yt J]-Yt J-ll f 
•/.SET  JX  3YETC  • 


*/.  OYtOl  IS  NOT  USED  IN  CODE. 


input; 


THE  CURRENTLY  ASSIGNED  LOCATIONS  IN  THE  2-BLOCK  ARE 


2. 

3. 

4. 

5. 

6. 


10. 

11. 

12. 

13. 

14. 

15. 

16. 
17. 


problem  no. 

CYCLE. 

TIME. 

JT. 


jM AX  MUST  BE  A  MULTIPLE  OF  JMAXPERSTRIP. 


SUBROUTINE 
BEGIN 
/. 

7. 

X 

■/. 

Z 
Z 
Z 
Z 
Z 
Z 
Z 
Z 
Z 
Z 
Z 
X. 

Z 
Z 
Z 
Z 
Z 
Z 
Z 
Z 
Z 

mode=true; 

SIMNRITE(LINE, *ENTER  INPUT*); 

ZI  READ(I4INOUTCOJ,NEXT) ; 

SINRE A J ( I4DISK0I0) , NEXT) |_Z  READS  Z,Xf 

z=next [ o i ; 

frgb  sgrabone cz,o ) ; 

cycle=grabone(Zi1); 

cyclastprint=cycle; 

t=grabone ( z  *  2  > ; 

3T=GRA30NE(Z*3) ; 

IMA  X  =GRAB0NE(Z*5) ; 

JMAX  =GRAB0NE<Z,6) ; 
eth=grabone<z,7) ; 


IMAX. 

JMAX. 

ETH. 

STABLFCTR,  SHOULD  BE  0.5. 

GO.  GRAVITY  AT  THE  EARTHS  SURFACE. 

TOZONE,  PERCENT  MASS  LEAVING  A  BOUNDARY  CELL  TO  TRIGGER 
A  REZONE. 

TOP  REFLECTIVE,  NONZERO  FOR  REFLECTIVE. 

RIGHT  REFLECTIVE,  NONZERO  FOR  REFLECTIVE. 

BOTTOM  REFLECTIVE,  NONZERO  FOR  REFLECT  IV-. . 

RELATIVE  ERROR  AT  LAST  PRINT. 

MAXIMUM  RELATIVE  ERROR  ALLOWED. 

CYCLE  STOP. 

JMAX  PER  HORIZONTAL  STRIP  PROCESSED. 


AN 3  Y.  CLOSE  PACKED. 


11111 

11111 

11111 

mu 
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STA3LFC1'R=  GRAB  ONE  (Z,8)  ? 

if{sta8lfctr<o.)  then  stablfctr=i. ; 

GO=GRABONE( Z,9> ? 

T  0 ZONE =GRA BONE (Z, 10 )  $ 

TOPREFLECTIVE  =  FALSE? 
n IGHTREFLECTI Vt  =  FALSE? 

80TTOMREFLECTI V£  =  FALSE? 

IF(GRABON£(Z,11>*0.)  THEN  TOPREFLECT IVE=TRUE? 

IF (GRA80NE ( Z , 12) X  0 . )  THEN  RIGHTREFLECTI VE=TRUE ? 

I F  (GRAUONE ( Z, 13) *  0 • )  THEN  BOTT OMR£FLECTIV£=TRUE» 
R£L£RR3EF0RE=GRA90N£ ( Z» 14 ‘ ? 

MAXR£LERR=bRA80NE (Z,15) ? 

C VCL ESTOP =GRA BONE <Z,16) ? 

JNAXPEK3TRIP=  GRABONE <Z,17) ? 

JSTRIPS= JMAX  JIV  JMAXPERSTRIP  $ 
JSTRIPSP1=JSTRIPS*1? 

X 0 LK= ( IMAX+1)  JIV  64? 


I  bLKPi  =1 8LK*1  ? 

I R0HS=IBLKP1* JMAXPERSTRIP? 

1kowsmi=irohs-i; 

IRCHS2=IROWS*IROHS? 

IROHS3=IROWS2*IROHS? 

IROWS4=IROWS3*IROHS? 

LCCP  1  =  0,1,  IBLK  DO  X  II)=NEXT  [  I  *1  ]  *, 

LOOP  J=0 , 1 , JMAX  TO  BEGIN 
J1=J  JIV  64? 

J2=J-64*J1? 

Yt J]=GRA0ONE (NEXTI Jl+IBLK+2) ,J2) ? 

END? 

T=T-OT? 

cycle=cycle-i? 

3IMWRTTE<LIN£,*TO  SET*) ? 

SETJXJYETC? 

SIMWRITEILINE, *BACK  FRCM  SET*)? 

XI  REAJ(I4IN0UT(  1]  ,LAST)  ? 

XI  R6AD(I4IN0UTI2J,NEXT) ? 

XI  RAIT(I4IN0UTtll) ? 

SIMREAD(14DISK0tl ],L AST) ? 

SIMREAJ(I4JISK0t2],NEXT) ? 

LOOP  JC=J,1,IR0HSM1  DO  BEGIN 
U  ( JC]=LASTtJC] ? 

V  t JC]=LAST[ JC+IROWS]? 

AMXt JC)=L ASTI JC* IRON S2]? 

AIX(JC]=LAST(JC*IROHS3] ? 

P  t JC]=LASTt JC+IR0WS41 ? 

cNJ? 

JPL=0? 

JPN=2? 

W0A=1.3t-20? 

JEOT  =  0  ? 

LOOP  J3=1,1,JSTRIPS  00  BEGIN 

SIMHRITE(LINE,*JB  ETC* , JB, JBOT , JPL ,JPN , AI X I  0 J , AIXti ) , AIXt 2 j , AIX 1 3 ) ) ? 
COT? 

XI  WAIT (I4SCRI JPLI) ? 

XI  WAIT ( 1 4l NOUT  t JPN1 ) ? 

JPL= JPL+1 $ 

JPN=JPN*1? 

IF  JPL  =JSTRl PSP1  THEN  JPL  =  1? 

IF  JPN=JSTRIPSP1  THEN  JPN=1? 

LOOP  JC=0,1, IROWSM1  00  BEGIN 


11111 

11111 

11111 

11111 

11111 


11111 

11111 


X  X 


LAST ( JC  ]=u  IJC15 

L AST [ JC  ♦ IROWS  1=V  [JC]; 

LAST[JC+IR0WS2J=AMX[ JCJ  5 
LAST  I JC +  IR0WS3 ) = AI X t JC ] 5 
L AST [ JC+IR0WS4) =P  [JC]? 
7. 

U  ( JC1=N£XT[JCJ ; 

\l  I  JC]  sNEXTIJC  +  IROWS  H 
AMX[JCl=NtXT[JC+IROHS2]} 
AIXC JC]=NEXT[ JC+IROWS3JJ 
P  [JCJ  =NEXTt JC+IR0WS4) ? 
EN3J  7 JC  LOOP. 

/.I  WRITEII4SCRI  JPL)  .LAST)  ; 

'/.I  REA3<I4IN0UTIJPN1,NEXT) 

SIM WRITE  (I40ISK7C JPL] .LAST) } 
SIMREA3 ( I43ISK0IJPN] .NEXT) J 
Jt)OT  =  JBOT*JMAXPERSTRIPJ 
tNJJ  X JB  LOOP. 

end;  xim>ut. 

7. 

7. 


11111 

11111 

11111 

11111 


SUBROUTINE  PHI? 

BEGIN 

cint  pass.jm;  creal  g; 

HQCE=TRUE  » 

SIMWRITEUINE.+INTO  PHI*)? 

oth=o; 

URRI0]=0. ; 

/.  SET  P  AT  AXIS. 

ROGE=BOOLcAN(-0) ; 

LOOP  K=O.IBLKP1,IROWSM1  30  PI K ]  =  RTL ( 1 , ,P[ K  ) )  ; 

IF  RIGHT  REFLECTIVE  THEN  ERRORU) 

ELSE  BEGIN 

7.  SET  P  AT  RIGHT. 

M0DE=REVR(IMAX+1-64*IBLK.M0DE) ; 

LOOP  K  =  I0LK.IBLKP1. 1  ROW SMI  JO  P(K]=RTR(2»  .PIK+IR2) ) J 

cno; 

LOOP  PASS=1 .1,2  00  BEGIN 

mooe=true; 

IF  JB=  0  THEN  IF  BOTTOMREFLECTIVE  THEN  LOOP  I  =  0,1,IBLK  30  EEGIN 

vb( i j-o . ; 

PBiiJspm; 

END 

ELSE  begin 

wca=  jv ( i ] ; 

wcaovm; 

WCC=1.1/(WCA+WCB); 

LOOP  I=0,1,IBLK  30  BEGIN 
MOCE=IMOOE[ II ; 

vb[ i i=vt i ) ; 

PB [ 1 1= (Pt I+IBLKP1]*WCA+Ptl ]*WCB)*WCC5 
OTH=OTH+0.5*HCA*WCC*OT*TAUm*V[I)*(PII) 
♦Ptl+IBLKPl)) ; 

end; 

enj; 

7.  SET  U  AT  AXIS. 

MOO E=  BOOLEAN  1-0)  *, 

LOOP  K=0 , IBLKP1, IR0WSM1  30  U[K)=-RTL<1, .UCK+IL] ) } 


MOOE=REVR ( IM AX  +  1-64* IBLK , MODE) { 

X  SET  U  AT  RIGHT. 

IF  RIGHTREFLECTIVE  THEN  E RROR (2 ) 

ELSE  BEGIN 

WCA=RTR  <l,,RCt IBLK+IHIT/RCt I BLK)  5 

LOOP  K= IBLK, IBL  KP1 , I R OH SMI  00  Ut<] =MCA*RTR < 1 , , Ut K+ IR ] ) J 

eno; 

V. 

X  HAVE  SET  ALL  BOUNDARY  CONDITIONS  BOTTOM, LFFT, RIGHT . 

X  NOW  00  SOME  PHYSICS. 

X 

WCA=0V[JB0T+11 ; 

k=o  ; 

IF  J8= JSTRIPS  THEN  JM=JMAXP£RST R IP-i  ELSE  JM=JMAXPERSTRIP ? 
LOOP  JJ=1,1,JM  DO  BEGIN 
J=JJ+JBCT? 

WC8=0Y[J+1] ; 

WCC=1.0/(HCA+WCB)  ; 

wco=wca*hcc; 

wce=wcb*wcc; 

WCF=YtJ]J 

LOOP  1  =  0,1, IBLK  00  BEGIN 
MOOE=IMOOE  t I] ? 

I F  ( AMX [ K] 5  0  • )  THEN  ERROR  (3 ) T 

URRtI  ]sU[K'J*HP2[  IJ+RTL  (1, , UCK+IL1)  *HP3TI1  ; 

MOCE=IMOOEL III; 

PR[Il=(PtK]*RTL(l,,OXII+IL])+RTL(l ,,P[K+IL 1 ) *0X111) 
♦HPitn; 

k=k+i; 

eno;  ziloopi. 

K=  K-IBL  KP1 * 

LOOP  1=0,1, IBLK  00  BEGIN 
HOOE=IMOOE  til ; 

VA  =  V[K]*WCE+VU+IBLKP1J*NC0; 
AIX[K]=AIXtK]+P[KJ*DTM0.5*  (VBt I]-V*)*TAUC 1 1 

+<RTR<1, ,URRCI+IR) ) -URR[ I ) ) *P I*WCA) /AMX ( K) 5 
IF  <PASS=1)  THEN  BEGIN 

PA=Pt  K J*HCE*P£  K+IBL  KP1  ]*NCOJ 
G=6.3567t8; 

G=G/ <  G+WCF-O. 5*MCA ) J 

g=go*g*g; 

UtK]  =  Ut<]  +  TAUm*WCA*  < RTR (1, , PR[ I+XR ] ) -PRt I ) ) 
*JT/(OXm*AMXTK]»  ; 

V[KJ=VCK]+TAU£I]*  (PBm-PA)*OT/AMX[K]-G*OT{ 
pb[I]=pa; 

enc; 

vbci j=va; 
k=k+i; 

END;  ZIL00P2. 

hca=hcb; 
end;  •/.  j  loop 

IF  J 8= JSTRIPS  THEN 'bYGIN  ZTOP  ROW  COOE. 

LOOP  1=0,1, IBLK  00  BEGIN 
MODE=IMOOEIIi; 

IF  (AMXtKlSO.)  THEN  ERROR  <i»)  ? 

URR[IJ=U[K]4WP2CIJ  +RTL ( 1 , , U [K  +  IL ] ) *WP3[  I ] ; 
MOOE=IHOOELIIJ5 

PR [ I ]=  <P[ K 1 *RTL< l, , OX  1 1 +ILJ)+RTLC 1 ,,Pt  K*I L ) ) *  Ox  III ) 
♦WP1IIJ5 
K--K+1? 


END;  I L00P3 • 

K=K-I8LKPi; 

IF  TOPREFLECTIVE  THEN  £RR0R(5) 

ELSE  BEGIN 

LOOP  1=0) i, IBLK  JO  BEGIN 
M03E=IM0DEIIi; 

DTH=OTH-Q . 5#T  AUt I ]*HCB*V£  K J*UT* (P! K] 
-PtK-IBLKPil )*HCC; 

k=k+i; 

END,* 

K  =  K- IBLKP1 » 

end; 

LOOP  I=0,1,IBLK  DO  BEGIN 

MooE=iHootm ; 

VA=VtKi; 

AlXt K]  =  AIX!KOP!K]*DT*  <0.5* < Vet  I ] - VA> *T AUtll 

MRTR(1,  ,URRt  I  +  IR])-URRm)*PI*HCA)/AMXCKl  ; 
IF  <PASS=1)  THEN  UtKl=UIK]+TAUII]*FCA 

MRTR(1,,PR!I*IR]>  -PRm)*DT/OXm*AMX[K]>  J 
k=k+i; 

end; 

end;  /!TOP  RON  COOF. 

END?  •/.  PASS  LOOP 
mcde=true; 

C  TH=ETH+ROWSUM  (DO)  ; 

end;  ZFHl. 

•/. 

7. 

7 

7. 

SUB-HiL  T  INE  PH<J  J 
BEGIN 

CIST  K15 

FRLAL  VATTOPINTERFACE.DONORMASSDIVIOEOBYOY , 

DONORTOPU,OONORTOPV,OONORTOPSPECIFIC ENERGY, 
ONEOVERAMXNEH.UATRIGHT INTERFACE,  AMXNEH,*' 

rezonetop=false;rezoneright=false;rezonebottoh=false; 

mode=true; 

SiMRRITE(LINE,*IN‘!0  PH2+>; 

DTH=Q  . ; 

HCA=OY[  JBOOl]  } 

V.  SET  AXIS  BOUNDARY  CONDITIONS. 

MASSTHRURIGHTt  J]=0. ; 

Y.  SET  BOUNDARY  CONDITIONS  AT  RIGHT. 

IF  RIGHT REFLECTIVE  THEN  ERROR (6) 

ELSE  BEGIN 

MODE=REVR(  IMAX  +  1-6<#*IBLK»  BOOLEAN  (-0)  >  J 
LOOP  K=I3LK,IBLKP1, IROWSMi  DO  BEGIN 
UIK]=RTR(i, ,U( K+IR] ) } 

IF  (U(<J<0.>  THEN  U(Kl=-UtKi; 

end; 

cnd; 

Y.  SET  BOTTOM  BOUNDARY  CONDITION. 

IF  JB=0  THEN  IF  BOTTOMREFLECTI VE  THEN  LOOP  1=0,1, IBLK  DO  BEGIN 
MODE=TRUEJ 

MASS THRU BOTTOM 1 1 1=0. J 
ENERGYTHRUBOTT  OMt 1 1  =  0 . ? 

UMOMENTUMTHRUBOTTOMt I )=0, 8 
V MOMENTUM  THRU BOTTOM! I J=0 . J 

END 
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ELSE  LOOP  I=0,1,IBLK  00  BEGIN 
M00c=IM00££IJ} 

Mf SSTHRU30TTCM£I1=V£I]*AMX£I )*JT/WCAJ 
IFJMASSTHRUBOTTOMII J>0.)  THEN  MASSTHRUBOTTOM! IJ  =  0. 5 
IF(- MASS THRU BO TTOM£I]>TOZONE*AMXIIJ)  THEN  REZONEBOTTOM=TRLEJ 
ENERGYTHRUBOTTONII  ] =M ASS THRU BOTTOM!  I1MAIXIU+0.5*  (UC  I )  *UC  I ) 

♦vm*vm>> ; 

oth=oth+ energy  thrubot  tomi  lit 

UMO MEN TUM THRU BOTTOM! I ]  =  MA  SS THRU BOTTOM! I ) *UC 1 1 » 
VHOMENTUMTHRLBOTTOM! I 3  =  MA SSTHRUBOTTOMIIJ * V!  I] 5 

enj; 

k=c; 

LOOP  JJ= 1,1,JKAXPERSTRIP  JO  BEGIN 
J=  JJ+JBOTJ 

7.  COMPUTE  QUANTITIES  AT  TOP  INTERFACE  OF  CELL. 

IF  J<JMAX  THEN  WCB=DY! J+i J { 

LOOP  1=0,1,  IBLK  JO  BEGIN 
MOOE= IMOOEI 11} 

SPECIFICENERGY!  I]=AIX£KJ+0.54(U£K]*IJ£K]+V£KJ*VIKJ)! 
JONORRIGHTSPECIFICENERGYII J-SPECIFICENERGYIII  ,* 
UONORMASSOIVIOEOBYOX!I]=AMXIK)  /  (DX  1 1  ]*P.CI  I )  > ' 

jonorrightuiij=uik); 

OONORRIGHTV£I]=VIKJ? 

K=K+15 

ENO  !  7. 1  LOOP  1 . 

K-K- IBLKP1 ! 

LOOP  1=0,1, IELK  00  BEGIN 
MOOE=IMOOEII]{ 

J0N0RT0PU=UIKJ5 
OONORTOFV=V£K] J 

JONORTOFSPECIFICENERGY=SPECIFICENERGYII  ]; 
OONORM4SSDIVIOECBYDY=AMXEKJ/WCA; 

IF  J< JMAX  THEN  BEGIN 
K1=K+IBLKP15 

VATTOPINTERFACE  =  0.5MV£K)«-V£K1] ) 5 
VATTOPINTERFACE=VATTOPINTERFACE/(1.0  +  (V£K1  J-VIK1) 
•OT/WCAI?  7.  THIS  GIVES  ENC  OF  SLUG  VELOCITY. 

IF  <VATTOPINTERFACE<0.)  THEN  BEGIN 

00N0RMASS0IVI0E0BY0Y=AMXIK1 J/WCBJ 
J0N0RT0PU=U£K1JJ 
00N0RT0PV=V£K1 JJ 
J0N0RT0PSPECIFICENERGY=AIX£K1 J 

♦  0.5MU£K1)*U£K1]+V£K1)*V£K1J)  $ 

enc; 

MASSTHRUTOPtI J = VAT TOPINTERF AC E* JT 
•OONORMASS  UIVIDEDBYDY } 

ENO 

ELSE  BEGIN  7.  00  TOP  ROW  DIFFERENTLY. 

IF  TOPREFLECTIVE  THEN  ERROR (7  > 

ELSE  BEGIN 

VATTOPINTERFACE=VIKJ } 

MASSTHRUTOP! I J=VATTOPINTERFACE*JT 
•DONORMASSOI VIDEO BY OV  J 

IF(MASSTHRUTOPII )<0.)  THEN  MASSTHRUTOP£IJ=0. 5 
IF  MASSTHRUTOPt I ] >TOZONE*AMX  £  K  J  THEN  REZONETOP 
=  TRUE  ? 

JTH=3TH+MASSTHRUT0P£I ]* JONOR TOPSPEC IF IC ENERGY? 

End; 

end; 

V.  COMPUTE  QUANTITIES  AT  RIGHT  INTERFACE  OF  CELLS. 
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UATRIGHT  INTERFACED.  5*  (U(KJ  +  RTL<1,  .UlK+ILl) )  } 
UATRIGHTINTERFACE=UATRIGHT INTERFACE/  (1. 0+(RTL (1.  ,UlK  +  IL] > 
-U£K])*JT/JXI  I]  >*,  X  THIS  GIVES  REAR  OF  SLUG  VELOCITY 
IF  ( UATR I GHt INTERFACE <0 • )  THEN  BEGIN 

OONORMASS  31 VI 3E  3BY JXt I ]  =  RTL  (1 . . 30NORMASS3I VI 3E3BY3X 
ti+iD) ; 

DONORRIGHTUII3=RTL (1, , OONORRIGHTUC I*IL  3 )  } 
30N0RRIGHTVCI]=RTL ( 1 , , DONORRIGHTVt I+IL ] ) 5 
OONORRIGHTSPECIFICENERGYII]* 

RTLU,  ,30N0RRI GHT  SPECIF ICENERGYl I+IL ] >  ! 

END} 

MASSTHRURIGHTt I ] =U AT  RIGHT  I  NT  ERF  ACE  * 3T 
*UONORMASSOIVIOEOBYOXm*xm  } 

X  MOHENTUM  FLUXES  AT  TOP  ANO  RIGHT 
UM0MENTUMTHRUT0PtI]  =  MASSTHRUT0Pm*30N0RT0Pu; 
UliOMENTUMTHRURIGHTt  I1=MASSTHRUR  I  GHT  1 1 1  *DONORRI  GHTUtll  *, 
VM0MENTUMTHRUT0PII]  =  MASSTHRUT0Pm*30N0RT0Pv; 
VMONfcNTUMTHRURIGHTl I) =HA SS THRUR IGHT 1 1 1  * DONORRI GHT V (I! 5 
ENERGYTFRUTOPI I ]=MASSTHRUTOP [ I]* 30NCRT0FSPECIFICENERCY ; 
ENERGYTHRURIGHTII1=MASSTHRURIGHT[I] 
•30N0RRIGHTSPECIFICENERGYC I)} 

k=k*i; 

END  *  ZIL00P2. 

K=K-I0LKPi; 

MODE-TRUE} 

J1=IKAX  OIV  64? 

J2=IMAX-64*J1} 

IF  (GRABONE (MASS THRU  RIGHT [ J1 ] , J2 ) > 

TOZONE*GRA0ONE<AMX(K+J1],J2) >  THEN  REZONERI GHT  =  TRUE  » 

IF  (PEN=J2)  THEN  PTH=OTH-ENERGYTHRURIGHTt Jl) 5 

LOOP  I~0 . 1 . I0LK  DO  BEGIN 
M03E=IMC3El I J ! 

AMXNEH=AMX(K]-MASSTHRUTOP(I J+MASSTHRUBOTTOMII] 

-MASSTHRURIGHTt I ] +RTR ( 1 »  »MA  SST HRURIGHTI I+lR] ) } 

oneoveramxnem=i.o/amxneh; 

•/.  UPDATE  U 

U[  K  ]  -  IAMXIK]*U(K3-UM0NENTUM  THRU TOP (13 

♦UMOMENTUMTHRUBOTTOMl I] -UMOMENTUMTHRURIGHT 1 1 ] 

♦RTR ( 1 » » UMOMENTUMT HRURIGHT 1 1 ♦ IR  1 ) ) ♦ONEOVERAMXNEWJ 
X  UPOATE  V 

V[K]  =  (AMX(K]*V(K ]- VMOMENTUM  THRUTOPl 1 1 

♦V MO MENTUMTHRU80TT0MII1- VMOMENTUM I HRURIGHT (  I ] 

♦  RTR (1 1  »V MOMENT UMT  HRURI GHT lI+IR])>*ONEOVERAMXNEK} 

X  UPDATE  A IX 

A IX ( K  ]  = ( AMX  CK! ’SPECIFIC ENERGY! I J  ^ENERGY  THRU TOP t 1 1 
♦ENERGYTHRUBOttONtI ]-ENERGYTHRURIGHTt  I ) 

♦RTR ( 1, .ENERGY T HRURIGHT tl+IR] ) -0.5*AMXNEW 
•IUI  K]*UIK]+VIK]*VIK]) ) ♦ONEOVERAMXKEWJ 
X  UPDATE  AMX 

AMXt  <j=amxnew; 

•/.  STORE  BOTTOM  QUANTITIES  FOR  NEXT  ROM 
UMCMENTUMTHRUBOTTOMt Il=UMOMENTUMTHRLTOPd]  } 
VMOMENTUMTHRUBO  TTOMf I ]= VMOMENTUM THRU TOP t 11} 
ENERGYTHRUBOTTOMmsENERGYTHRUTOPtl]  J 
MASS  THRU BOTTOMl I ]=MASSTHRUTOPt I )  } 

k*k+i; 

END;  XILOOP3. 

end;  xjloop. 

*ooe=true; 

£  TMETH+ROWSUN  <DT HI  ; 
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END;  V.  PF2. 

X 

X 

X 

/• 

SUBROUTINE  SUiiTOTflL ENERGY? 

BEGIN 

t'OJEsTRUE  ? 

wpfl=o • ; 

LOOP  I=0,1,IBLK  30  BEGIN 
HOOE=IMOOEtI  1? 
k = i ; 

LOOP  J=1,1,JMAX  30  BEGIN 

WPA=NPA*ANXlKlMAIXtKl*0.5MUIK]*UCK]*V[K  J*V[K1 )) 5 
K=  K* I BLKP1 ? 
end;  •/.  j  loop 
end;  x  i  loop 
moje=true; 

TOT  ALENERGYNOW=T0T  ALENERGYNOW+ROWSUM  CHPA) } 

end;  /.sumtotalenergy, 

X 

X 

X 

suewotTiNE  output; 

BEGIN 

creal  relerrnow,engycheck,numberofcycles; 
moge=true; 

SIMWRITE(LINE,*INTO  OUTPUT*)? 

XL  READ(I4SALUTEIQ],Z SALUTE)  ? 

XI  WAIT<I4SALUTE[01) ? 

LOOP  1=1, 1,6  DO  IF  GRABONt (ZSALUTE  » 1-1 ) =0  THEN 
SW[ I ]  =  F  ALSE  ELSE  SWII5=TRUE? 

IF  PEN=6  THEN  ZSALUTE=CYCLEJ 
IF  PEN=7  THEN  ZSALUTE=T? 

IF  PEN=8  THEN  ZSALUTE=DT? 

7.1  HRITE(  14 SA LUTE!  1)  , ZSALUTE)  ? 

XXX 

y.x’x 

XXX  HERE  PULSE  OSK  TO  LET  IT  EXECUTE  SALYUT  IN  DURING  CLAUSE  OF 

XXX  ICL.  TRANSFER  SALUTE  INFORMATION  FROM  AN  ILLIAC  IV  FILE  TC  THE 

XXX  OUTSIDE  WORLD  FOR  EXAMINATION, 

XXX 

XXX 

X 

X 

X 

IF  CYCLE=  TOUMP  OR  STOPTHI SCYCL E  THEN  BEGIN 
RELERRNOW=(TOTALENERGYNOW-ETH)/ETH? 
NUMBEROFCYCLES=CYCLE-CYCLASTPRINT  J 
IF  <NUMBEROFCYCLES=0.0)  THEN  NUMBER0FCYCLES  =  1 . 0  ? 
tNG  YCHECK= (RELERRNOW-REL  ERRBEFORE) /NUMBEROFC YCLES? 
R£LERRBEFORE=RELERRNOW; 

CYCLASTPRINT=CYCLE? 

IF  PEN=1  THEN  7=CYCLE? 

IF  PEN=2  THEN  Z=TJ 
IF  PEN=3  THEN  Z=DT? 

IF  PEN=7  THEN  Z=ETHJ 

IF  PEN=8  THEN  z=stablFCTr; 

IF  PEN=14  THEN  Z=REL ERRBEFORE ? 

XI  WAIT (I4IN0UTI JPL1) ? 
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7.  1 


11111 


WAIT  II4SCRIJPL])  ? 

LAST  C  Q ]  =  Z ; 

LOOP  1=0,1, IBLK  DO  LAST! I*1J=XI II' 

LOOP  j=o,i,Jmax  DO  BEGIN 
J1=J  OIV  64? 

J2=J-64*J1? 

IF  J2=PEN  THEN  LAST! Jl* IBLK+2 1 =DY I J] ? 

end; 

7.1  WRITEII4INOUTIO],Z>;  Hill 

SIMWRITEI I4DISK0I0] ,Z> ?  11111 

7.7.7.  INTERUPT  HERE  FOR  OSK...  WE  HAVE  OATA  IN  FILE  I4IN0UT  TO  BE 

77.7.  COPIED  EITHER  TO  THE  LASER  STORE  OR  TO  TAPE.  A  PROGRAM  IN  THE 

7.77.  DURING  CLAUSE  OF  ICL  WILL  DO  THIS  COPY.  ***MUST***  LOCK  OUT 

7.77.  FURTHER  WRITES  BY  GLYPNIR  TO  I4INOUT  TILL  THIS  COPY  IS  COHPLETE 

7.7.7.  HENCE  MUST  HAVE  A  RESULT  DESCRIPTOR  FROM  OSK. 

7.77. 

7.7.7  THIS  IS  NOT  YET  OEFINEO.  ♦••NOTE***  PROEAELY  BETTER  TO  FIRST 

7.7.7  PULSE  THE  COPY  OF  THE  FILE  OUT,  RETURN  THIS  DESCRIPTOR,  ANC 

77.7  THEN  PULSE  AN  EOIT  OF  THE  COPIEO  FILE  FOR  PRINTER  OUTPUT. 

7.7.7  THIS  USE  OF  TWO  TASKS  IN  THE  DURING  CLAUSE  OF  ICL  WOULO  ALLOW 

7.7.7  MORE  OVERLAP,  ANO  WE  NEED  ALL  WE  CAN  GET . 

7.7.7 

IF  RELERRNOW»MAXRELERR  THEN  ERROR(200i? 

END? 

end;  '/.OUTPUT. 

•/. 

7. 

7 

7. 

SUBROUTINE  REZONE; 

7.  NB.  REZONE  WILL  READ  ALL  OF  I4SCR  AND  REWRITE  IT  AFTER 

/.  THE  REZONE.  IF  CYCLE=TOUMP  OR  STOPTHISC YCLE  IS  TRUE  THEN 

7.  THE  FILE  I4IN0UT  MUST  BE  REWRITTEN  ALSO. 

ERROR  (100)  ; 

7 

7. 

7. 

7. 

7.  MAIN  PROGRAM  STARfs  HERE. 

7. 

7. 

SHELL 


7.1  OPEN  ( 141 NOUT )  ;  11111 

7.1  0PEN(I4SCR)  J  .  11111 

7.1  OPEN ( I4SALUTE)  ;  11111 

INPUT?  */  READ  Z-BLOCK,X,  AND  Y.  COPY  HYORO  TO  I4SCR. 

7.1  READII4SCRI1  ] , L AST)  J  Xii.ll 

7.1  READ(I4SCRI21,NEXT)  ;  11111 

71  wait;  11111 

SIMREADII4DISK7I 11,LAST) 5  _  11111 

SIMREAD(i40ISK7I2],NEXT) }  11111 


LOOP  JC=0,1,  IROWSM1  DO  BEGIN 
U  I JC 1 =L ASTI JC ] $ 

V  l JC ] =L ASTI JC+IROWS  1? 
AMXl JC1=LASTIJC*IR0WS2] J 
AIXI JC1=LASTI JC+IR0WS3]; 
P  I JC1=L ASTI JC+IR0WS4] J 
7 

U  I JC*1R0WS1=NEXTI JCi; 
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V  [JC+IROWS]=NEXT[ JC+IROWS  ] 
AHX( JC+ IRONS] =NE XT [ JC+IR0WS2] 

$ 

AIXt JC+IROWS] =NEXT( JC+IR0WS3 ] 
P  [ JC+ IRO WS1= NEXT [ JC+IROWS 4] 

end;  7.JC  LOOP. 

7.1 

READ(I4SCRt31,NEXT)  *, 
SIrtREAO(l4DISK7t31.NEXT)  *, 

AuAlN 

choosedt; 

jpl=o; 

JPN  =  3{ 

JBOT=OJ 

LOOP  JB=1 ,1, JSTRIPS  DC  BEGIN 

phi; 
ph2  ; 
cot; 

mode=true; 

IF  CYCLE=TDUMP  THEN  SUMTOT ALENERGY 5 
7.  I  WAIT  ( 14  SCR t JPL  ] )  ?  X  LAST  WRITE. 

7.1  WAIT  (I4SCRI JPN1 )  t  •/.  NEXT  READ. 

jpl=jpl*i; 

JPN=JPN+1J 

IF  JPL  = JSTRIPSP1  THEN  JPL=15 
IF  JPN=JSTRIPSP1  THEN  JPN=1? 

LOOP  JC-0  1 1. IRO WSH1  00  BEGIN 
lastijc  ]=u  c  jc i ; 

LAST  I JC+IROWS  I=V  [JCi; 

LAST! JC*IR0WS21=AMX[ JCI 5 
LAST! JC+IR0WS11=AIXI JCI 5 
LAST! JC+IR0WS41=P  tJCi; 

7. 

U  [JC]=U  IJC+IROWSi; 

V  [JC1=V  ( JC+IROWS ) ? 

ArtXt JC1=AMXI JC+IROWS] ? 
AIX(JC1=AIXIJC ♦ I ROWS] ; 

p  tjci=p  ijc+irowsi; 

x 

U  tJC+IROWS]=NEXTtJC] » 

V  [JC+IROWS]=NEXTl JC+IROWS  )5 
AH XI JC+IROWS] =  N£XTt  JC+ IR0WS2J  } 
AIXt  JC+IROWS] = NEXT t  JC  +  IROWS 3) } 

P  { JC+IROWS]=NEXTt JC+IR0WS4) ; 

END  $  7.JC  LOOP. 

IF  CYCLE=TOUMP  OR  STOPTHISCYCLE  THEN 


7.1 

WRITE<I4IN0UT[ JPLJ.LAST) ? 

SI MWRITE (I  4DISK0I JPL] » L AST) » 

XI 

WRITE CI4SCRI JPLl, LAST); 

7.1 

READ (I4SCRt  JPN1 » NEXT ) { 
SIMWRITE(I40ISK7t JPL], LAST)  *, 

1 

S IMREA D  <I4DISK7 l JPN 1 .NEXT)  ? 

jbot=jbot+jmaxperstrip; 
end;  7.JB  LOOP. 

uui kui  » 

IF  REZONETOP  OR  REZONERIGHT  OR  REZONEBOTTOM  THEN  REZONE  $ 
IF  NOT  STOPTHISCYCLE  THEN  GO  TO  AGAIN  $ 

TCP  FIN  END. 
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11111 

11111 

11111 

liill 

11111 
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APPENDIX  XII 

CODE  REQUIRED  TO  TRANSFER  A  B6500  CLAM  OUTPUT  FILE  TO  AN  ILLIAC  IV  INPUT  FILE 


SL  1ST 

FILE 

FILE 


SINGLE 

4=CLAM/0UT,UNIT  =  DISK, RECORDS  100  0 
1=CLAM/CUTG,UNIT=OISK,SAVE=1,AREA=6Q,R£CORD=540 


PROGRAM  I4CLAM  FOR  SHELL /OF/THE/FUTURE. 

THIS  DECK  WILL  CONVERT  ANY  SIZE  CLAM  WRITTEN  1000  WDS/RECORD. 


REAL  Z(150),U(200),V(20Q) , AMX ( 20 0) , A I X (200) ,P (20 0 ) »X ( 1J03  )  , 
C  Y (15  CO) *  ZZ(  16000  ) 

EQUIVALENCE  (IHAX,Z ( 33) ) ,  ( JMAX , Z ( 35) ) ,  (KHAXA  ,Z ( 3ft) ) 

INTEGER  STRIPS 
LOGICAL  FIRST 
DATA  FIRST/. TRUE./ 


STRIPS^  5 
IROWSH=5C 
KLH=  5*64* IROWSH 

READ (4)  WS  ,  WSA 
READ (4)  Z 

JRO  WS=  JMAX-1 
JPS=JROWS/STRIPS 
IF (JPS*STRIPS.EQ.JROWS)  GO  TO  6 
WRITE (6,5) 

FORMAT (15H  BAD  CLAM  DATA.) 

STOP 

K 1=  1 
10  KL  =  KMAX  A 

IF(KL.GT.K1+199)  KL=K1+199 

READ  (4)  OMY 

K1=KL*1 

IF  (KL.LT . KMAXA)  GO  TO  10 

c 

READ  (4)  Xn,<xm,DMY,I=l,lMAX) 
READ  (4 )  YO,  (Y(J)  ,J=1,JMAX) 

C 

DO  20  K=1,KLH 
20  ZZ(K)=0 

ZZ(1)=Z(1> 

ZZ (2) =Z (2 ) 

ZZ ( 3)  =  Z  (84) 

ZZ(4)=Z(3) 

ZZ (5)=Z (26) 

ZZ  (6)  =  I(*AX 
ZZ (7)= JROWS 
ZZ  (8)  =  Z (13) 

ZZ(9)=Z(139) 

ZZ(10) =  Z  ( 108) 

ZZ  ( 1 1)  =  Z (75) 

ZZ  (12) =0 
ZZ (13) =0 
ZZ (14) =1 
ZZ ( 15)  =  Q 
ZZ  ( 16)  =1 
ZZ  <  1 7)  =  1 . 

ZZ (18) =  JPS 
C 
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00  40  I =1 » IN AX 
40  ZZ(H-65>=X<I) 

C 

IBLKP1= (IMAX+65) /64 
IW0S1=64*I8LKP1 
IROWS=I ELKPl* JPS 
IW0S=64*IR0WS 
iwnsr=iwosi+64 
C 

KU  =  0 

KV=KU+IH0S 

KM=KV+INOS 

KI=KM*IMOS 

KP=KI+IWJS 

C 

DO  60  J  =  1 t JMAX 
60  ZZ(J+IWCSY»=Y(J> 

C 

CALL  HRTI4D (ZZ» KLH) 

C 

REMIND  4 
READ (4)  OMY 
READ (4 )  DMY 
DO  50  Kal.KLH 
50  ZZ(K)~0 

C 

KN=i 

K=2Q0 

DO  100  J=?t JMAX 
KL=KN 

JO  95  1  =  1)  IMAX 
70  IF(K.LT .200)  GO  TO  80 

READ  (4)  (U(KK» , V(KK) , AMX ( KK ) , AIX ( KK> , P ( KK) ,KK  =  1,20  0) 

K=K-20Q 

IF  (FIRST )  K=IMAX*1 
FIRST=. FALSE. 

GO  TO  70 
C 

80  KsKtl 

KN=KN+1 
ZZ  (KNtKU) =U( K) 

ZZ (KN+KV) =V ( KJ 
ZZ  (KN+KM) = AMX (K ) 

ZZ (KN+KI) =AIX (K) 

ZZ(KN*<P)=P<K> 

95  CONTINUE 
KN=KL*IWDS1 

IF ( ( ( J- 1 > / JPS) * JPS.NE. J-l>  GO  TO  100 
CALL  WRTI4D(ZZ«KLH) 

DO  96  1=1, KLH 

96  ZZ(I)=0 
KN=  1 

100  CONTINUE 
LOCK1 
STOP 
END 

SUBROUTINE  WRTI4D(Z,N) 

DIMENSION  1(540)  ,ZTN> 

WRITE<6,99> 

99  FORMAT ( 1H1) 
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11  =  0 
NN=N 
5  K=NN 

IF (< .GT .256)  K=  256 

00  11  1=1, < 

L=LL+I 

IFCZtU.EQ.O)  GO  TO  10 
WRIT  £ (6 ,98)  L,Z(L),Z(L) 

98  FORMAT (I1G,F20«8,E2Q,8) 

CALL  B5I4F(Z(L) ,A(2*I-1) ,A (2*I>) 
GO  TO  11 

10  A  <2*1-1 )=  0, 

A  (2*1)  =0  • 

11  CONTINUE 

HRITt(l)  (A(I) ,1=1,540) 

NN=NN-K 

IF(NN.EC.O)  RETURN 

LL=LL+2E6 

GO  TO  5 

END 

SUBROUTINE  85I4F(A,B,C) 
EQUIVALENCE  (XX, II) 

C 

XE=CONC AT (0, A, 6, 45, 7) 

IF(Xt.EO.O)  GO  TO  10 
CALL  BI4F(A,B,C) 

RETURN 

C 

10  XX= A 
X  =  II 
X=X+0.1 
X=X-0.1 

CALL  BI4F (X, B,C) 

RETURN 

END 

SUBROUTINE  BI4F(A,B,C> 
EQUIVALENCE  (XE,IE8) 

DATA  I  BIAS/0400  0  0/ 

IF(A.NE.O)  GO  TO  10 

B=0 

C=0 

RETURN 

C 

10  XS=CONCAT (0, A, 0,45,1) 

XE=CONC AT  (OjAj.5,44,6) 

IF(XS.NE.O)  IE8=-IE8 
IE8=-IE8 

B9=C0NCAT(0, A, 0,38,1) 

B10=C0NC AT (0, A, 0,37,1) 

IE8E=2 

IF(B10.NE.0LIE8E=1 
IF(B9.NE.0>  IE8E=0 
IE=3*IE8+IE8E 
IE=IBIAS-IE+39 
BB=CONC AT (0, A, 31, 46,1) 

BB=CONCAT (BB, IE ,30 ,14,15) 

B=CONC AT (BB,A,15,38-IE8E,16) 
C=CONC AT  ( 0  ,  A,31 .22-IE8E  ,23-IEA'. ) 
RETURN 
END 


1(W 


SUBROUTINE  BI4HA,B,C) 
BB=CONCAT(Q,Af 31,46,1> 
B=CONCA T  < 8B *  A , 6, 38 , 7) 
C-CONCAT  < P  » A , 3i  »  31 . 32) 
RETURN 
END 


o  o  o  o  o 
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APPENDIX  XIII 

CODE  REQUIRED  TO  EDIT  AN  OUTPUT  DUMP  FROM  SHELL/OF/THE/FUTURE 


FILE  4=SH£LL/I4CLAM/0UT6QBY15,UNIT=DlSKtREC0RD=54Q 

EDIT  FOR  SHELL/OF/THE/FUTURE. 

COMMON  2(16000) 

DIMENSION  U(10QQ)tV(10QQ)»AMX(lQQQ)»AIX(100Q)*P(10QO)* 

C  X(IOOO) ,  3X(  1000) ,TAU( 1000)  ,Y  (1500) >JY (1500) 


IN0ST=5*64*5Q 
C 

205  CALL  REJI4J (Z» IWDST ) 

CALL  PRINTZ 
C 

PROB-Z ( 1 ) 

CYCLE=Z  12} 

T=Z ( 3) 

DT=Z  (4) 

IMAX=Z(6) 

JMAX=Z(7) 

ETH=  Z (8 ) 

RELERR=Z ( 15) 

JPS*Z (18) 

STRIPS3 JMAX  /JPS 

IBLKPl=(IMAX*65)/64 

IWOSl=64*iBLKPl 

IR0HS=IBLKP1*JPS 

IWJ$=64*IR0MS 

IH0SY=IH0Sl«-64 

KU=0 

KV=KU+IWOS 
KM3KV*INJS 
KI=KM*IWOS  • 

KP=XH-IHOS 

C 

JO  10  I=1,IMAX 

10  X(I)3Z(I*65>  - 

C 

WSA»0 

JO  15  I31 » IHAX 
OX (I)=X (I)-MSA 
15  MS A  3  X ( I ) 

C  . .  .  -  — 

JO  20  J=1 i JMAX 
20  Y ( J) -Z ( J+IM3SY ) 

C 

WSA=0 

JO  25  J=l* JMAX 

OY(J)  =  Y(J)-MSA  _  -  —  - 

25  MSA*Y ( J) 

C 

WSA  =  0 

JO  30  1=1 t IMAX 
HSB=X(I)**2 

TAU(I)= 3. 14159265358929* (WSB-MSA) 
30  MSA»MSB 

C 

KK=0 
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30  600  JB=1, STRIPS 
CALL  RE  31 43 (Z» I  HOST ) 

KL  =  0  1 _  _ 

00  610  J=1,JPS 
30  620  I=1*IMAX 
KK=KK+1 
K=KL4l*l 
U (KK)-Z (K+KU) 

V(KK)=Z(K+KV) 

AHX(KK)=Z(K+KH) 

AIX ( KK) =Z ( K+KI ) 

P(KK)=Z (K+KP! 

620  CONTINUE 

KL=KL*-IH0S1 
610  CONTINUE 
600  CONTINUE 
C 

C  SHORT  PRINT 

C 

HR"Tt<6,470>  CYCLE 
00  200  1=1,4 
200  PR(I)=0 

K=0 

00  221  J=1,JHAX 
DO  220  I~1  » I  MAX 
K=K4  1 

IF ( ANX ( K) 1300,220,210 
210  WS0*O.5MU(K>»*2+V(K)**2> 

PR ( 1 )=PR ( 1 )♦  ANX (K)*AIX(K) 

PR(2)=PR(2)*WSB*AMX(K) 

PR(4)=PR!4)+ANX(K) 

220  CONTINUE 

221  CONTINUE 

PR( 3)=  PR ( 1) +PR ( 2! 

HRITE  <6 ,3601  PRO0, CYCLE, T,DT 
WRITE  <6,3701  (PR ( I ) , 1*1 , 4) 

WRIT E < 6 , 380 >  ETH ,REL  ERR 
C 

C  LONG  PRINT  C03E. 

C 

HRITE (6 ,360)  PROB, CYCLE DT 

00  200  1= 1 , IMAX 

HRITE (6 ,420)  I,X(I),3X(I) 

K=IMAX* JHAX+I 

J=JNAX4l 

00  270  L  =  1„JNAX 

J=J-1 

K=K-I1AX 

IF ( AhX ( K ) • EQ.O • )  GO  TO  270 
RHO=AMX (K)/(TAU(I)*OV(J) 1 

HRITE  (6, 430)  J,U(K>  »  V  (K )  ,P (K)_,3HY  ,ANX  (K)  jAIX(K)  ,  RHO  ,  Y(  J) 
270  CONI INUE 
283  HRITE( 6,600) 

GO  TO  205 


C 

C 

C 

300 


NEGATIVE  NASS. 


C 


hRITE(6,490) 

STOP 


107 


o  o 


360 
370 
38  0 
423 

430 

470 

490 

900 


C 

99 


20 

10 


C 


100 


10 

11 


C 

200 

201 


C 

C 


c 

5 

C 

10 

C 


FORMATS. 


FORMAT  (  8HiPR0BLEN6X5HjSYJCLE9X4HTIME13X2HDT/F9.4  »F11. 0  ,lP2X2r.  1.6 . 7 

FORMAT (////17X2HAI16 X2HAK14X5HAI* A K15X 2HAM/1P7X4E1 8.7) 

FORMAT  (///16X3HTHE12X9HREL  ERROR/1P7X2E 18.  7//// » 

FORMAT ( 4H  I  =I3,6X6HX<I)  =  F 12 . 3 , 6X7H0X (I>  =F12.3//3H  J8X1HU 

X13X 1HV1 3X 1HP 14X4X1 1X3HAMX 11 X3HAI XI 2X3HRH011X1HV/) 
F0RMAT(1P1XI3,1X8E14.6) 

FORMAT (///21H  TAPE  4  OUMP  ON _CTCLEF6 .0////> 

FORMAT (27H1NEGATIVE  MASS  ENCOUNTERED.  > 

FORMAT  (/////> 


ENO 

SUBROUTINE  PRINTZ 
COMMON  Z (64) 


WRITE (6  »99) 
FORMAT  (1H1) 
OO  23  1=1.32 


RETURN 

END 

SUBROUTINE  RE3I4D<Z,N) 
REAL  Z(N> ,ZZZ(540) 
LOGICAL  FIRST 
OATA  FIRST/. TRUE./ 


IF (FIRST)  REHINO  4 
FIRST=. FALSE. 

N 1=  0 

READ (4,END=2C0)  (ZZZ  (I) » 1=1 »540> 

00  11  K=1 .256 
I  — N 1 lK 

CALL  I4BF(ZZZ(2*K-1) ,ZZZ(2*K) *Z  €I>  > 

IF(I.GE.N)  RETURN 

CONTINUE 

Nl=Nl+256 

GO  TO  100  .  - 


WRITE(6 .201)  „„„ 

FORMAT (24H1EN0  OF  FILE  ON  14  OISK.) 

STOP 

ENO 

SUBROUTINE  I4J BFCA.B.C)  - 

I NT EGER 01 00 

OATA  IBIAS, 0100/040000, 0100/ 
EQUIVALENCEfXE.IE) 


I F ( A . N E •  0 ._0_. OR, B, N E ._Q. 0 )  _GO_TO  5 - 

C*0 

RETURN 

XE=!'ONCAT<0,A,14,30,15> 

IE=I8I AS+39-IE 

IE8  IS  OOHERJ3F  8  _  .  _ 

IIE8EEIS  RT  SHIFT  OF  MANTISSA  AT  IE8 
IE8E=IE-3*IE8 
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IF(IESE.GE.O)  GO  TO  20 

IE8E=IE8E+3 

I E  8 = I E  8 -1  _  _ 

20  IE8='-IE8 

IF (IE8.LT >0)  IE8=0100-IE8 
C  PICK  UF  SIGN 

OUT  =CONCAT (0, A, 46,31,1) 

C  PICK  UP  EXPONENT 

OUT  =CONCAT (OUT , IE8 ,45  ,6 , 7 ) 

C  PICK  UP  TOP  OF  MANTISSA 

Ot;?-CONCAT(OUT,A,38-IE8E,15,i6> 

C  PICK  UP  BOTTOM  OF  MANTISSA 

C=CONCAT(OUT,B,22-IE«E,31,23-IE8E) 

RETURN 

END 

SUBROUTINE  I4BI(A,B,C> 
IF(A.NE.O.O.OR.B.NE.O.O)  GO  TO  5 
C*0 

RETURN 

C 

5  CC=CONCATtO, A, 46,31,1) 

CC*CONCAT (CC , A , 38 ,6 , 7 ) 
C=CONCAT(CC,B,31,31,32> 

RETURN 

END 
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